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SECTION  1 
INTRODUCTION 


Timely  design  and  analysis  of  missiles  requires  the  use  of  rapid  and 
accurate  analytical  procedures  to  determine  their  aerodynamic  stability,  con¬ 
trol  and  drag  characteristics.  The  most  complete  compendium  of  methods 
available  is  the  USAF  Stability  and  Control  Datcom,  Reference  1.  Since  it 
emphasizes  aircraft  configurations,  it  does  not  address  the  range  of  geometric 
and  flight  condition  parameters  unique  to  missiles,  such  as  high  angle  of 
attack  and  bank  angle.  Although  many  aerodynamic  methods  are  available  in 
the  literature  for  missile  configurations,  there  is  no  collection  of  these 
techniques  in  a  form  suitable  for  efficient  missile  design. 

This  study  did  determine  the  feasibility  of  developing  a  Missile  Datcom. 
Specifically,  the  four  objectives  of  this  study  were  as  follows: 

(a)  Determine  the  range  of  geometric  and  flight  conditions  for  which 
methodology  should  be  specified. 

(b)  Determine  a  structure  of  a  handbook  and/or  computer  program  which 
allows  rapid  and  accurate  use  of  available  methodology  in  the 
missile  design  environment. 

(c)  Survey  the  literature  for  applicable  methodology  and  determine 
those  areas  lacking  in  appropriate  methods  and  needing  further 
development. 

(d)  Assess  the  feasibility,  or  probability  of  success,  of  developing 
a  Missile  Datcom  and  recormend  a  handbook  and/or  computer  program 
format. 

This  report  addresses  each  of  these  objectives.  Throughout  the  report 
reference  is  made  to  the  conceptual,  preliminary  and  point  design  phases  of 
missile  development.  The  definitions  of  these  phases  are  as  follows: 

Conceptual  Design  -  The  design  process  during  which  the  proof  of  concept 
is  demonstrated.  This  task  demonstrates  the  feasibility  of  designing  a 
missile  system  which  performs  the  tasks  required,  such  as  payload  carry 
capability,  range,  speed  and  altitude. 

Preliminary  Design  -  This  stage  of  the  design  process  explores  several 
variations  of  a  configuration.  Trade  studies  are  identified,  individual 
components,  such  as  nose  shape  and  wing  type,  are  analyzed  and  the  impact 
on  system  design  due  to  configuration  aerodynamics  is  determined. 
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Point  Design  -  This  is  the  last  major  step  in  the  design  process  prior 
to  hardware  development.  A  baseline  configuration  is  selected  and  per¬ 
formance  characteri sties  determined.  Perturbations  to  the  design  are 
often  performed  to  accommodate  changes  in  subsystem  design  or  to  provide 
an  extension  of  capability. 

To  be  useful.  Missile  Datcom  should  apply  to  all  three  design  phases,  but 
emphasize  conceptual  and  preliminary  design. 

Since  the  missile  design  process  is  often  of  shorter  duration  than  that 
for  aircraft,  quick  and  accurate  methods  are  needed.  Automation  of  the 
techniques  fulfills  the  quickness  requirement.  However,  low  analysis  costs 
are  also  important.  These  costs  include  the  man-hours  required  to  set-up, 
execute  and  interpret  the  results  of  the  program.  The  advantage  of  sophisti¬ 
cated  accurate  theoretical  methodology  must  be  weighed  against  the  costs  re¬ 
quired  to  obtain  the  results.  Determining  the  proper  choice  between  accuracy, 
efficiency  and  cost  was  a  goal  of  this  study. 

An  automated  program  does  not  provide  the  aerodynamicist  with  the  neces¬ 
sary  background  or  methods  choices  available.  A  thorough  user's  manual,  or 
a  handbook  similar  to  the  USAF  Stability  and  Control  Datcom,  will  supply 
this  additional,  yet  essential  information.  Since  methods  more  suitable  to 
accurate  design  may  be  difficult  to  employ  in  a  handbook  technique,  the  choice 
between  alternate  methodology  and  generation  of  computer-based  design  charts 
was  also  explored. 

The  methodology  collected  was  assessed  based  upon  the  criteria  given  in 
Table  1.  Ideally,  the  methods  selected  should  be  theoretically  based,  have 
minimal  parameter  inputs,  be  accurate,  and  cover  a  wide  range  of  configuration 
variables  and  flight  conditions.  Since  no  single  method  will  meet  each  of 
these  goals,  the  eight  criteria  of  Table  1  will  help  identify  those  techniques 
that  should  be  retained  for  the  quantitative  assessment  phase  of  Missile  Datcom 
development . 
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TABLE  1  ASSESSMENT  CRITERIA 


1.  METHOD  APPROACH 

—  Theoretical 

—  Semi-Empirical 

—  Empirical 

2.  EFFICIENCY  (HANDBOOK,  COMPUTER) 

—  Number  Computations 

—  Complexity  of  Logic 

—  Number,  Type  of  Inputs 

—  Iterative 

—  Detail  of  Geometry  Required 

3.  ACCURACY 

—  Existing  Validation 

—  Compatibility  with  Accuracy  Requirements 

—  Sensitivity  of  Output  to  Input  Accuracy 

—  Derivation  Assumptions  —  Theoretical 

—  Range  of  Data  —  Empirical 

—  Geometric  Model 

4.  STATUS 

—  Current  Use  in  Industry 

—  Handbook  Method  Available 

—  Method  Coded  and  Used  Locally 

—  Does  it  Need  Modification 

—  Is  Modification  State-of-the-Art 

5.  RANGE  OF  APPLICABILITY 

—  Flight  Conditions  (Mach,  a,  <p) 

—  Geometry 

6.  UTILITY  OF  OUTPUT  PARAMETERS 

—  Compatibility  with  Other  Methods 

—  Thoroughness 

7.  GENERAL  UTILITY 

—  Understandability 

—  T  raceability 

—  Modifiability 

8.  VALIDATION  STATUS 

—  Existing  Validation 

—  Data  Base  Available 

to  Complete  Validation 
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SECTION  2 


REQUIREMENTS  AND  MISSILE  DATCOM  ARCHITECTURE 


2.1  Requirements 

To  be  a  useful  design  tool,  Missile  Datcom  must  address  those  geometric 
and  flight  condition  requirements  of  interest  to  the  missile  design  community. 
As  illustrated  in  Figure  1,  configuration  shape  is  often  described  by  factors 
beyond  the  control  of  the  aerodynamicist,  such  as  seeker,  warhead  size, 
launching  platform  or  cost.  In  addition,  missile  designs  are  strongly  driven 
by  the  mission  to  be  performed.  As  illustrated  in  Figure  2,  the  type  of 
mission  influences  the  choice  of  components  which  make  up  the  configuration 
design.  The  matrix  of  configurations  which  satisfy  the  mission  requirements 
are  comprised  of  both  mission  and  design  requirements.  The  vehicle  class 
(cruise  versus  intercept)  and  configuration  type  (glide,  boost-glide  or  boost 
sustained)  are  choices  which  depend  upon  the  mission  to  be  accomplished,  and 
are  usually  range  or  speed  dependent.  On  the  other  hand,  the  selection  of 
individual  configuration  components,  such  as  wing  size,  body  shape,  or  tail 
arrangement,  arise  from  both  physical  (launcher  constriants,  steering  mode 
or  payload  size  and  weight)  and  mission  dependent  (propulsion  type)  require¬ 
ments.  Hence,  a  range  of  configurations  can  be  specified  which  satisfy  the 
mission  goals;  such  a  wide  range  of  design  variables  must  be  addressed. 

Selection  of  the  physical  properties  and  flight  conditions  most  useful 
to  the  missile  design  community  was  performed  through  analysis  of  the 
world's  missile  systems.  Their  characteristics  were  extracted  from  "Jane's 
Weapon  Systems"  (reference  2),  and  the  experimental  data  summaries  of  the 
"Aeromechanics  Survey  and  Evaluation"  report  (Reference  3).  Design  experience 
obtained  from  numerous  in-house  missile  programs  have  also  been  considered. 

The  most  useful  results  were  obtained  through  categorizing  the  data  of 
the  Aeromechanics  Survey.  The  types  of  configurations  tested  do  reflect 
current  trends  in  missile  designs.  The  classification  of  results  are  shown 
in  Figure  3,  and  were  obtained  from  1824  separate  test  summaries.  These 
results  indicate  that  the  typical  missile  design  is  a  body-tail  configura¬ 
tion  of  overall  fineness  ratio  10.  It  is  tested  across  the  Mach  spectrum  to 
approximately  30  degrees  angle  of  attack  and  the  lifting  surfaces  are  straight- 
tapered  with  an  aspect  ratio  between  one  and  four.  A  sharp  tangent-ogive 
is  the  most  common  nose  shape.  Since  this  typical  missile  is  too  specific 
for  design  purposes,  two  sets  of  requirements  evolved  from  these  results. 
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1)  Priority  1  capability,  the  range  of  parameters  which  Missile  Datcom  must 
address,  and  2)  Priority  2  capability,  the  range  of  parameters  which  extend  the 
utility  of  the  method  collection,  but  are  not  detrimental  to  the  development 
of  Missile  Datcom  or  its  usefulness.  Coverage  of  the  priority  1  range  of 
conditions  was  used  to  determine  feasibility  of  Missile  Datcom. 

The  priority  1  range  encompasses  75%  of  the  conditions  identified  in 
the  Aeromechanics  Survey  plus  additional  conditions  based  upon  in-house 
design  experience.  This  results  in  the  range  of  conditions  given  in  Table  2. 

One  crucial  parameter  is  Mach  number.  Mach  number  limitations  vary  with 
altitude,  and  are  primarily  driven  by  structural  and  thermodynamic  considera¬ 
tions.  The  Mach-altitude  boundaries  selected  are  shown  in  Figure  4.  The 
boundaries  are  based  upon  structural  requirements,  thermodynamic  properties 
of  materials,  speed  required  to  maintain  flight,  and  effective  aerodynamic 
control.  Included  on  this  figure  are  the  requirements  for  the  Reynold's 
number  per  foot  of  length  from  the  1962  U.S.  Standard  atmosphere. 

2.2  Missile  Datcom  Architecture 

The  most  common  missile  aerodynamic  analysis  technique  is  “component 
build-up".  A  configuration  is  analytically  modeled  as  a  combination  of 
components  such  as  body,  wings,  tails  and  inlets.  The  aerodynamic  coeffi¬ 
cients  of  each  of  the  components  are  estimated,  interference  effects  among 
them  determined  and  all  coefficients  summed  to  determine  full  configuration 
aerodynamics.  This  technique  is  advantageous  in  the  preliminary  or  con¬ 
ceptual  design  process  since  many  configuration  components  are  screened  to 
establish  the  best,  configuration. 

The  structure  of  Missile  Datcom  as  presented  in  Table  3  in  outline 
form  has  been  assembled  assuming  a  handbook  format  and  using  the  "component 
build-up"  approach.  A  handbook  version  of  Missile  Datcom  is  necessary  to 
supply  the  detailed  documentation  necessary  for  effective  use  of  the  methods 
selected.  There  are  nine  sections  in  the  handbook. 

Section  1  will  include  a  synopsis  of  the  Missile  Datcom  methods  contained 
in  the  entire  document.  The  synopsis  gives  the  user  a  quick  overview  of  the 
available  methodology  and  the  method  limitations. 

Section  2  will  1)  define  the  standardized  notation  used  throughout  the 
document,  2)  supply  equations  or  charts  for  calculation  of  geometric  characteris¬ 
tics,  such  as  nose  mold  line  contours,  wetted  area,  planform  area,  panel  local 
chord,  and  3)  provide  flight  condition  parameters,  such  as  Reynolds  number  and 
speed  of  sound  at  altitude. 
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The  remaining  sections  will  contain  the  method  descriptions  covering 
the  following  areas: 

(a)  a  concise  description  of  the  physical  phenomena  modeled, 

(b)  a  description  of  the  method  including  equations,  tables  or  charts 
necessary  for  its  use, 

(c)  tables  or  charts  which  illustrate  the  accuracy  of  the  method,  and 

(d)  a  bibliography  of  reference  documentation  or  related  material. 

Since  the  methods  will  be  automated,  approximate  equations  will  also  be 
presented  to  enable  quick  hand  calculations. 

Section  3  through  9  have  been  structured  in  a  heirarchy  of  component 
class,  speed  regime,  and  aerodynamic  parameter.  This  heirarchy  will 
minimize  redundant  explanations,  and  enables  computation  of  all  component 
forces  and  moments  within  one  section.  This  approach  also  serves  as  a 
guide  to  a  modular  computer  program  development. 

The  notation  used  in  missile  design  is  standardized  by  usage  convention. 
No  known  written  standard  is  in  use,  although  "A  Compilation  of  Aerodynamic 
Nomenclature  and  Axes  Systems",  NOLR  1241  (Reference  4)  was  published  some 
18  years  ago  for  use  in  the  U.S.  Navy.  However,  the  USAF  Stability  and 
Control  Datcom  symbols  and  nomenclature  are  in  wide  use  in  aircraft  design 
and  is  recommended  for  missiles.  With  few  exceptions,  this  set  of  nomen¬ 
clature  allows  little  ambiguity  and  is  applicable  to  missiles.  Some 
exceptions  are:  the  reference  length  should  be  defined  as  the  maximum  body 
diameter/width/height,  and  the  reference  area  as  the  area  of  an  equivalent 
circular  section  whose  diameter  is  the  reference  length. 

By  convention,  the  body  axis  system  is  in  primary  use  in  government 
and  industry.  This  convention  should  be  accepted  as  the  Missile  Datcom 
standard.  It  is  recognized  that  other  axes  systems  would  be  preferable 
for  a  specific  task,  and  conversion  capability  among  axes  systems  is 
required.  Static  and  dynamic  axes  system  conversion  equations  are  recom¬ 
mended  for  inclusion;  reference  documentation  is  given  in  Reference  4  to 
8.  It  is  also  convenient  to  express  the  relationship  between  pitch  and 
yaw  angle  of  attack  to  total  angle  of  attack  and  bank  angle  in  equation 
and  graphical  form.  A  proposed  form  of  this  information  is  given  in 
Figure  5. 

Since  missile  performance  is  a  result  of  aerodynamic  analysis,  it 
is  proposed  that  suitable  data  from  the  1962  U.S.  Standard  Atmosphere, 
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Standard  Day,  and  from  the  Military  Standard  210A/B,  non-standard  day, 
be  included  for  ready  reference.  A  typical  means  of  presenting  the  results 
is  shown  in  Table  4.  Automated  atmospheric  routines  are  readily  available 
for  computer  program  use,  with  an  average  error  of  less  than  0.5%.  The 
routine  in  Digital  Datcom  is  an  excellent  example  of  the  Standard  Day 
routine. 

Readily  available  equations  for  the  geometric  characteristics  of  a 
missile  configuration  will  facilitate  the  design  process.  Computation  of 
wetted  area,  nose  mold  line  contours,  and  panel  sweep  angle  are  examples 
of  those  fundamental  relationship  that  should  be  included. 

Thrust-drag  accounting  requirements  occur  when  the  aerodynamicist  and 
the  propulsion  engineer  derive  system  performance.  The  aerodynamicist  is 
responsible  for  estimating  drag.  The  propulsion  engineer,  must  determine 
installed  engine  performance.  To  avoid  a  complex  bookkeeping  task  and  avoid 
double  accounting  of  some  drag  components,  the  Thrust-Drag  Accounting 
Committee  of  the  JANNAF  (Joint  Army-Navy-NASA-Air  Force)  Airbreathing 
Propulsion  Working  Group  has  recommended  a  standard  procedure.  This 
procedure,  documented  in  the  "Airbreathing  Propulsion  Manual",  CPIA/M6, 
Reference  9,  is  recommended  for  inclusion. 

The  Datcom  provides  useful  methods  for  predicting  missile  mass  and 
inertial  characteristics.  This  information  is  recommended  for  inclusion 
because  of  its  utility  in  preliminary  design. 

The  remaining  sections  of  this  report  will  describe  Sections  3  through  9 
of  the  Missile  Datcom  Outline  and  present  the  method  recormiendations . 
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Figure  1.  Factors  Which  Affect  Missile  Design 
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Figure  2.  Matrix  of  Possible  Configurations 
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Figure  3.  Aeromechanics  Survey  Results 
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c)  Mach  Number 


d)  Angle  of  Attack 


Figure  3.  (Continued) 
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f)  Nose  Fineness  Ratio 


Figure  3.  (Continued) 
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g)  Wing  Planform 
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h)  Tail  Planform 

Figure  3.  (Continued) 
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i)  Wing  b/d 


j)  Tail  b/d 


Figure  3.  (Continued) 
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k)  Wing  Aspect  Ratio 


1 )  Tail  Aspect  Ratio 

Figure  3.  (Continued) 


16 


m)  Panel  Combination 


n) 


Control  Device  Type 


Figure  3.  (Continued) 
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TABLE  2  RANGE  OF  GEOMETRIC/FLIGHT  CONDITIONS 


PARAMETER 

SYMBOL 

PRIORITY  1 

PRIORITY  2 

SOURCE 

ANGLE  OF  ATTACK,  DEG. 

a 

-20<q<T0 

-180<a<180 

AEROMECHANICS  SURVEY  AND 

ANGLE  OF  YAW,  DEG. 

B 

-2O<0<2O 

-180<6_480 

WORLD’S  MISSILE  SYSTEMS 

AERODYNAMIC  ROLL,  DEG. 

0<|  4>|  <45 

(Ki^isiso 

i 

MACH  NUMBER 

M 

0<M<6 

0<M40 

j 

BODY  FINENESS  RATIO 

UAOb 

6<(V(1)q<2Q 

l<(2/d)B<30 

NOSE  FINENESS  RATIO 

(Vd)N 

.5<(£./d)^5 

0^(2.  /d )  ^<7 

FIN  EXPOSED  SPAN  TO  DIAMETER 

b/d 

l<b/d<6 

ft_b/d<10 

FIN  ASPECT  RATIO 

AR 

0.6<AR<4 

0 .  kAR<10 

FIN  PLANFORM 

TRIANGULAR 

TRAPEZOIDAL 

ALL 

WING/TAIL  ORIENTATION 

IN-LINE 

ALL 

CONTROL  METHOD 

ALL  MOVEABLE 

FIN 

ALL 

REYNOLDS  NUMBER/FT 

rN 

3xl05<Rm<2xl07 

103<Rm<3x107 

MACH-ALTITUDE  BOUNDARY 

FIN  DEFLECTION/INCIDENCE,  DEG. 

6 

0<6<30 

0«c<5  <60 

MISSILE  SYSTEM  ANALYSIS 

ROLL  RATE,  RAD/SEC. 

P 

0<|p|<l 

0<|p|<8 

l 

PITCH  RATE,  RAD/SEC. 

q 

0<|q|<1.5 

0<|q|<3 

YAW  RATE,  RAD/SEC. 

r 

0<|r|<1.5 

0<|r|<3 

FIN  DEFLECTION  RATE,  RAD/SEC. 

6 

0<|«|<10 

0<|  il< 28 
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Re/ft. 


MACH  NUMBER 


Figure  4. 


Mach-Altitude-Reynolds  Number  Requirements  for  the  1962 
Standard  Atmosphere 
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TABLE  3  MISSILE  DATCOM  OUTLINE 


1.  INTRODUCTION  AND  METHODS  SUNNARY 

2.  GENERAL  INFORMATION  AND  EQUATIONS 

2.1  NOTATION 

2.2  AXES  SYSTEMS  AND  TRANSFER  EQUATIONS 

2.3  BODY  PARAMETERS 

2.4  LIFTING  SURFACE  PARAMETERS 

2.5  ENGINE  AND  INLET  PARAMETERS 

2.6  SIMPLIFIED  EQUATIONS  OF  NOTION 


3.  CHARACTERISTICS  OF  BASIC  BODIES 

3.1  AXISYNMETRIC  BODIES 

3.1.1  SUBSONIC 

3. 1.1.1  AXIAL  FORCE 

3. 1.1. 2  NORMAL  FORCE 

3. 1.1. 3  PITCHING  MOMENT,  XAC,  XCP 

3. 1.1. 4  ASYMMETRIC  VORTEX  SHEDDING  AT  HIGH  ALPHA/BETA 

3.1.2  TRANSONIC 

3. 1.2.1  AXIAL  FORCE 

3. 1.2. 2  NORMAL  FORCE 

3. 1.2. 3  PITCHING  MOMENT,  XAC,  XCP 

3.1 .2.4  ASYMMETRIC  VORTEX  SHEDDING  AT  HIGH  ALPHA/BETA 

3.1.3  SUPERSONIC/HYPERSONIC 

3. 1.3.1  AXIAL  FORCE 

3. 1.3. 2  NORMAL  FORCE 

3. 1.3. 3  PITCHING  MOMENT,  XAC.  XCP 

3.2  TWO-AXIS  SYMMETRICAL  BODIES 

3.2.1  SUBSONIC 

3. 2. 1.1  AXIAL  FORCE 

3. 2. 1.2  NORMAL  FORCE  AND  SIDE  FORCE 

3. 2. 1.3  PITCHING  (XAC,  XCP),  YAWING  AND  ROLLING  MOMENTS 

3. 2. 1.4  ASYMMETRIC  VORTEX  SHEDDING  AT  HIGH  ALPHA/BETA 

3.2.2  TRANSONIC 

3.2.2. 1  AXIAL  FORCE 

3. 2. 2. 2  NORMAL  FORCE  AND  SIDE  FORCE 

3. 2. 2. 3  PITCHING  (XAC,  XCP),  YAWING  AND  ROLLING  MOMENTS 

3. 2. 2. 4  ASYMMETRIC  VORTEX  SHEDDING  AT  HIGH  ALPHA/IETA 

3.2.3  SUPERSONIC/HYPERSONIC 

3.2.3. 1  AXIAL  FORCE 

3. 2. 3. 2  NORMAL  FORCE  AND  SIDE  FORCE 

3. 2. 3. 3  PITCHING  (XAC,  XCP),  YAWING  AND  ROLLING  MOMENTS 

3.3  ARBITRARY  SHAPED  BODIES 

3.3.1  SUBSONIC 

3. 3. 1.1  AXIAL  FORCE 

3. 3. 1.2  NORMAL  FORCE  AND  SIDE  FORCE 

3. 3. 1.3  PITCHING  (XAC,  XCP),  YAWING  AND  ROLLING  MOMENTS 

3.3.2  TRANSONIC 

3.3.2. 1  AXIAL  FORCE 

3. 3.2.2  NORMAL  FORCE  AND  SIDE  FORCE 

3.3.2. 3  PITCHING  (XAC,  XCP),  YAWING  AND  ROLLING  MOMENTS 

3.3.3  SUPERSONIC/HYPERSONIC 

3. 3.3.1  AXIAL  FORCE 

3. 3.3.2  NORMAL  FORCE  AMD  SIDE  FORCE 

3. 3. 3.3  PITCHING  (XAC,  XCP),  YAWING  AND  ROLLING  MOMENTS 
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TABLE  3  (CONTINUED) 


3.4  EFFECT  OF  PROTUBERANCES 

3.4.1  SUBSONIC 

3. 4. 1.1  AXIAL  FORCE 

3*4. 1 •?  NORMAL  FORCE  AND  SIDE  FORCE 

3. 4. 1.3  PITCHING  (XAC,  XCP),  YAWING  AND  ROLLING  MOMENTS 

3. 4. 1.4  ASYMMETRIC  UORTEX  SHEDDING  AT  HIGH  ALPHA/ BETA 

3.4.2  TR ANSON IC 

3. 4.2.1  AXIAL  FORCE 

3.4.2.?  NORMAL  FORCE  AND  SIDE  FORCE 

3. 4. 2. 3  PITCHING  (XAC,  XCP).  YAUING  AND  ROLLING  MOMENTS 

3. 4. 2. 4  ASYMMETRIC  UORTEX  SHEDDING  AT  HIGH  ALPHAS BETA 

3.4.3  SUPERSONIC/HYPERSONIC 

3.4.3. 1  AXIAL  FORCE 

3. 4.3.2  NORMAL  FORCE  AND  SIDE  FORCE 

3. 4. 3. 3  PITCHING  ( XAC.  XCP>#  YAWING  AND  ROLLING  MOMENTS 

3.5  PROPULSION  SYSTEM  EFFECTS 

3.5.1  IASE/FXHAUST  PLUME 

3. 5. 1.1  SUBSONIC 

3.5. 1.2  TRANSONIC 

3. 5. 1.3  SUPERSONIC /HYPERSONIC 

3.5.2  AIRBREATHING  INLET  AERODYNAMIC  CONSIDERATIONS 

3.5.2. 1  SUBSONIC 

3.5.2.2  TRANSONIC 

3. 5. 2. 3  SUPERSONIC/HVPERSONIC 

3.5.3  INLET  COUERS 

3.5.3.1  SUBSONIC 

3. 5. 3. 2  TRANSONIC 

3. 5. 3. 3  SUPERSONIC/HYPERSONIC 

3.5.4  INLET  PLUGGED 

3.5.4.1  SUBSONIC 

3.5.4.2  TRANSONIC 

3.5.4.3  SUPERSONIC/HYPERSONIC 

4.  CHARACTERISTICS  OF  LIFTING  SURFACES 

4.1  SECTION  CHARACTERISTICS  AND  DESIGN 

4.1.1  FLAT  PLATE  SECTIONS 

4. 1.1.1  SUBSONIC 

4. 1.1. 2  TRANSONIC 

4. 1.1.3  SUPERSONIC/HYPERSONIC 

4.1.2  DOUBLE  WEDGE  SECTIONS 

4. 1.2.1  SUBSONIC 

4. 1.2. 2  TRANSONIC 

4. 1.2. 3  SUPERSONIC/HYPERSONIC 

4.1.3  HEXAGONAL  SECTIONS  (MODIFIED  DOUBLE  WEDGE) 

4.1.3. 1  SUBSONIC 

4. 1.3. 2  TRANSONIC 

4. 1.3. 3  SUPERSONIC/HYPFRSONIC 

4.1.4  CIRCULAR  ARC  SECTIONS 

4. 1.4.1  SUBSONIC 

4. 1.4.2  TRANSONIC 

4. 1.4.3  SUPERSONIC/HYPERSONIC 

4.1.5  NACA  AIRFOIL  SECTIONS 

4.1.5.1  SUBSONIC 

4.1.5.2  TRANSONIC 

4.1.5.3  SUPERSONIC/HYPERSONIC 

4.1.6  SUPERCRITICAL  AIRFOIL  SECTIONS 

4. 1.6.1  SUBSONIC 

4. 1.6. 2  TRANSONIC 
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TABLE  3  (CONTINUED) 


4. 1.6. 3  SUPERSONIC/HYPERSONIC 

4.2  THREE-DIRENSIONAL  EFFECTS 

4.2.1  STRAIGHT  TAPERED  SURFACES 

4. 2. 1.1  SUBSONIC 

4. 2. 1.2  TRANSONIC 

4. 2. 1.3  SUPERSONIC/HYPERSONIC 

4.2.2  CRANKED  SURFACES 

4.2.2. 1  SUBSONIC 

4. 2. 2. 2  TRANSONIC 

4. 2. 2. 3  SUPERSONIC/HYPCRSONIC 

4.2.3  SLEWED  SURFACES 

4.2.3. 1  SUBSONIC 

4. 2. 3. 2  TRANSONIC 

4. 2. 3. 3  SUPERSONIC-' HYPERSONIC 

4.2.4  CURVED  SURFACES 

4.2.4. 1  SUBSONIC 

4. 2. 4. 2  TRANSONIC 

4. 2. 4. 3  S'JPERSONIC/HYPERSONIC 

4.2.5  VERY  LOU  ASPECT  RATIO  SURFACES 

4.2.5.1  SUBSONIC 

4.2.5.2  TRANSONIC 

4.2.5.3  SUPERSONIC/HYPERSONIC 

S.  INTERFERENCE  EFFECTS 
S.t  CARRYOVER 

5.1.1  BODY  /  FIN.  FIN  /  BODY 

5. 1.1.1  SUBSONIC 

5.1.1. 2  TRANSONIC 

5. 1.1. 3  SUPERSONIC ''HYPERSONIC 

5.1.2  ADJACENT  FIN  EFFECT  (FIN  /  FIN) 

5. 1.2.1  SUBSONIC 

5.1.2. 2  TRANSONIC 

5.1.2. 3  SUPERSONIC/HYPERSONIC 

5.1.3  BODY  /  INLET,  INLET  /  FIN 

5. 1.3.1  SUBSONIC 

5. 1 .3.2  TRANSONIC 

5.1.3.3  SUPERSONIC/HYPERSONIC 

5.1.4  NULT1PLE  BODIES 

5. 1.4.1  SUBSONIC 

5. 1 .4.2  TRANSONIC 

5. 1 .4.3  SUPERSONIC/HYPERSONIC 

5.2  VORTICES 

5.2.1  VORTEX  STRENGTH  AND  TRACKING  —  BODY 

5.2. 1.1  SUBSONIC 

5.2. 1.2  TRANSONIC 

5.2. 1.3  SUPERSONIC/HYPERSONIC 

5.2.2  VORTEX  STRENGTH  AND  TRACKING  —  LIFTING  SURFACES 

5.2. 2.1  SUBSONIC 

5.2. 2. 2  TRANSONIC 

5.2. 2. 3  SUPERSONIC/HYPERSONIC 

5.2. 3  VORTEX  STRENGTH  AND  TRACKING  —  INLET 

5.2. 3.1  SUBSONIC 

5.2. 3. 2  TRANSONIC 

5.2. 3. 3  SUPERSONIC/HYPERSONIC 

5.2. 4  NULTIPLE  VORTEX  INTERFERENCE 

5.2.4. 1  SUBSONIC 

5.2. 4. 2  TRANSONIC 

5. 2. 4. 3  SUPERSONIC/HYPERSONIC 
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S.  CONFIGURATION  SYNTHESIS 

6.1  BASIC  BODY  AND  INLET/ENGINE 

6.1.1  SUBSONIC 

6.1. 1.1  AXIAL  FORCE 

6. 1.1. 2  NORMAL  FORCE  AND  SIDE  FORCE 

6. 1.1. 3  PITCHING  (XAC,  XCP).  YAWING  AND  ROLLING  NONENTS 

6. 1.1.4  ASYMMETRIC  UORTEX  SHEDDING  AT  HIGH  ALPHA/BETA 

6.1.2  TRANSONIC 

6. 1.2.1  AXIAL  FORCE 

6. 1.2.2  NORMAL  FORCE  AND  SIDE  FORCE 

6. 1.2.3  PITCHING  (XAC.  XCP),  YAWING  AMO  ROLLING  MOMENTS 

6. 1.2. 4  ASYMMETRIC  UORTEX  SHEDDING  AT  HIGH  ALPHA/BETA 

6.1.3  SUPERSONIC/HVPERSONIC 

6. 1.3.1  AXIAL  FORCE 

6. 1.3.2  NORMAL  FORCE  AND  SIDE  FORCE 

6. 1.3. 3  PITCHING  (XAC,  XCP),  YAWING  AND  ROLLING  MOMENTS 

6.2  BODY-TAIL  OR  BODY-WING 

6.2.1  SUBSONIC 

6. 2. 1.1  AXIAL  FORCE 

6.2. 1.2  NORMAL  FORCE  AND  SIDE  FORCE 

6.2. 1.3  PITCHING  (XAC,  XCP),  YAWING  AND  ROLLING  MOMENTS 

6.2. 1.4  ASYMMETRIC  UORTEX  SHEDDING  AT  HIGH  ALPHA/BETA 

6.2.2  TRANSONIC 

6.2.2. 1  AXIAL  FORCE 

6. 2. 2. 2  NORMAL  FORCE  AND  SIDE  FORCE 

6. 2. 2. 3  PITCHING  (XAC,  XCP).  YAWING  AND  ROLLING  MOMENTS 

6. 2.2. 4  ASYMMETRIC  UORTEX  SHEDDING  AT  HIGH  ALPHA/BETA 

6.2.3  SUPERSONIC/HVPERSONIC 

6.2.3. 1  AXIAL  FORCE 

6. 2. 3. 2  NORMAL  FORCE  AND  SIDE  FORCE 

6. 2. 3. 3  PITCHING  (XAC,  XCP).  YAWING  AND  ROLLING  MOMENTS 

6.3  BODY-WING-TAIL  OR  BODY-CANARD-UING-TAIL 

6.3.1  SUBSONIC 

6. 3. 1.1  AXIAL  FORCE 

6.3. 1.2  NORMAL  FORCE  AND  SIDE  FORCE 

6. 3. 1.3  PITCHING  (XAC.  XCP).  YAWING  AND  ROLLING  MOMENTS 

6. 3. 1.4  ASYMMETRIC  UORTEX  SHEDDING  AT  HIGH  ALPHA/BETA 

6.3.2  TRANSONIC 

6.3.2. 1  AXIAL  FORCE 

6. 3. 2. 2  NORMAL  FORCE  AND  SIDE  FORCE 

6. 3. 2. 3  PITCHING  (XAC,  XCP).  YAWINC  AND  ROLLING  MOMENTS 
6. 3. 2. A  ASYMMETRIC  UORTEX  SHEDDING  AT  HIGH  ALPHA/BETA 

6.3.3  SUPERSONIC/HVPERSONIC 

6. 3. 3.1  AXIAL  FORCE 

6. 3. 3. 2  NORMAL  FORCE  AND  SIDE  FORCE 

6. 3. 3. 3  PITCHING  (XAC,  XCP).  YAWING  AND  ROLLING  MOMENTS 

7.  CONTROL  DEUICES 

7.1  ALL  MOUABLE  LIFTING  SURFACES 

7.1.1  BASIC  AERODYNAMICS 

7.1. 1.1  SUBSONIC 

7. 1.1. 2  TRANSONIC 

7. 1.1. 3  SUPERSONIC/HYPERSONIC 

7.1.2  HINGE  MOMENTS 

7. 1.2.1  SUBSONIC 

7. 1.2. 2  TRANSONIC 

7. 1.2. 3  SUPERSONIC/HYPERSONIC 
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7.1.3  I ENDING  NONE NTS 

7.1.3. 1  SUBSONIC 

7. 1.3. 2  TRANSONIC 

7. 1.3. 3  SUPERSONIC /HYPERSONIC 
7. a  PLAIN  TRAILING  EDGE  PLAPS 

7.2.1  BASIC  AERODYNAMICS 

7.2.1. 1  SUBSONIC 

7. 2. 1.2  TRANSONIC 

7. 2. 1.3  SUPERSONIC/HYPERSONIC 

7.2.2  HINGE  NONE NTS 

7.2.2. 1  SUBSONIC 

7. 2. 2. 2  TRANSONIC 

7. 2. 2. 3  SUPERSONIC/HYPERSONIC 

7.2.3  BENDING  MOMENTS 

7.2.3. 1  SUBSONIC 

7. 2. 3. 2  TRANSONIC 

7. 2. 3. 3  SUPERSONIC /HYPERSONIC 

7.3  JET  INTERACTION/REACTION  JET  CONTROL  (JI/RJC) 

7.3.1  BASIC  AERODYNAMICS 

7.3.1. 1  SUBSONIC 

7. 3. 1.2  TRANSONIC 

7. 3. 1.3  SUPERSONIC/HYPERSONIC 

7.4  THRUST  VECTOR  CONTROL  (TUC) 

7.4.1  BASIC  AERODYNAMICS 

7. 4. 1.1  SUBSONIC 

7. 4. 1.2  TRANSONIC 

7. 4. 1.3  SUPERSONIC/HYPERSONIC 

7.5  SPECIALIZED  CONTROL  DEUICES 

7.5.1  BASIC  AERODYNAMICS 

7.5.1. 1  SUBSONIC 

7. 5. 1.2  TRANSONIC 

7.5.1. 3  SUPERSONIC/HYPERSONIC 

7.5.2  HINGE  MOMENTS 

7.5.2. 1  SUBSONIC 

7. 5. 2. 2  TRANSONIC 

7.5. 2. 3  SUPERSONIC/HYPERSOHIC 

7.5.3  BENDING  MOMENTS 

7. 5. 3.1  SUBSONIC 

7.5. 3.2  TRANSONIC 

7. 5. 3. 3  SUPERSONIC/HYPERSONIC 

7.6  CONFIGURATION  TRIM  AERODYNAMICS 

t.  DYNAMIC  DERIUATIUES 

8.1  BODY  OR  BODY-INLET-ENGINE 

8.1.1  PITCH  DERIUATIUES 

8. 1.1.1  SUBSONIC 

8. 1.1.2  TRANSONIC 

8. 1.1. 3  SUPERSONIC/HYPERSONIC 

8.1.2  YAU  DERIUATIUES 

8. 1.2.1  SUBSONIC 

8. 1.2. 2  TRANSONIC 

8. 1.2. 3  SUPERSONIC/HYPERSONIC 

8.1.3  ROLL  DERIVATIVES 

8. 1.3.1  SUBSONIC 

8. 1.3. 2  TRANSONIC 

8. 1.3.3  SUPERSONIC/HYPERSONIC 

8.2  FIN  OR  FIN  COMBINATIONS 

8.2.1  PITCH  DERIUATIUES 
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8. a.  1.1  suisomc 

8. 2. 1.2  TRANSONIC 

8. 3. 1.3  SUPERSONIC/HYPERSONIC 
8. a. a  yau  derivatives 

8. a. a. l  subsonic 
8. a. a. a  transonic 

8. a. a. 3  SUPERSONIC/HYPERSONIC 

8.2.3  ROLL  DERIVATIVES 

8. 2. 3.1  SUBSONIC 
8. a. 3. a  TRANSONIC 

8.  a. 3. 3  SUPERSONIC/HYPERSONIC 

8.3  BODY-WING  OR  BODY-TAIL 

8.3.1  PITCH  DERIVATIVES 

8. 3. 1.1  SUBSONIC 

8. 3.1. 2  TRANSONIC 

8.3. 1.3  SUPERSONIC/HYPERSONIC 

8.3.2  YAU  DERIVATIVES 

8.3.2. 1  SUBSONIC 
8. 3. 2. 8  TRANSONIC 

8. 3.  a. 3  SUPERSONIC/HYPERSONIC 

8.3.3  ROLL  DERIVATIVES 

8. 3. 3.1  SUBSONIC 

8. 3. 3.  E  TRANSONIC 

8. 3. 3. 3  SUPERSONIC/HYPERSONIC 

8.4  BODY-WING-TAIL  OR  BOOY-CAHARD-U1NC-TAIL 

8.4.1  PITCH  DERIVATIVES 

8. 4. 1.1  SUBSONIC 

8. 4.1. 2  TRANSONIC 

8. 4. 1.3  SUPERSONIC/HYPERSONIC 

8. 4. 2  YAU  DERIVATIVES 

8. 4.2.1  SUBSONIC 

8. 4.  a. a  TRANSONIC 

8. 4.  a. 3  SUPERSONIC/HYPERSONIC 

8.4.3  ROLL  DERIVATIVES 

8.4.3. 1  SUBSONIC 

8. 4. 3. 2  TRANSONIC 

8. 4. 3. 3  SUPERSONIC/HYPERSONIC 

9.  SPECIALIZED  PROBLEMS 

9.1  TUNBLlNG  NOTION 

9.1.1  PLATES 

9. 1.1.1  SUBSONIC 

9. 1.1. 2  TRANSONIC 

9. 1.1. 3  SUPERSONIC/HYPERSONIC 

9.1.2  CYLINDERS 

9.1.2.1  SUBSONIC 

9.1.  a. a  TRANSONIC 

9.1.  a. 3  SUPERSONIC/HYPERSONIC 

9.1.3  CONES 

9. 1.3.1  SUBSONIC 

9. 1.3. 2  TRANSOHIC 

9. 1.3.3  SUPERSONIC/HYPERSONIC 

9.2  SPINNING  NOTION 

9.2.1  NAGNUS  EFFECT 

9. 2. 1.1  SUBSONIC 
9-a.l.a  TRANSONIC 

9.  a. 1.3  SUPERSONIC/HYPERSONIC 
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YAW  ANGLE  OF  ATTACK.  -0to«o) 


RELATIONS  BETWEEN  BODY  AXES 
ANGLES,  a  ANO  2.  AND  TOTAL 
ANGLE  OF  ATTACK,  a': 

0  «  -  ARCTAN  (TAN  a' SIN  0} 
a  *  ARCTAN  (TAN  a'  COS  0) 

a'  -  ARCTAN 

a+TAN1# 


PITCH  ANGLE  OF  ATTACK,  a  (deg) 


Figure  5. 


Conversion  Between  Pitch  and  Yaw  Angle 
Angle  of  Attack 


of  Attack  to  Total 
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TABLE  4  PRESENTATION  OF  STANDARD  ATMOSPHERE  DATA 


Altitude  «t) 

Preuure  (ib/vq  ft) 

Temperature  (deg  R) 

Speed  of  Sound  (ft/sec) 

0 

2.1 162003 

618.G7 

111645 

6000 

1.7609E  *03 

500.84 

1097.09 

10000 

1.455GE  *03 

483.03 

1077.40 

15000 

1.1948E  »03 

465.22 

1057.35 

20000 

9  7327E  *02 

447.42 

1036.93 

25000 

7.8633E  *02 

429  62 

1016.10 

30000 

6.29G6E  *02 

411.84 

994,85 

35000 

4.99 34E  *02 

394.06 

973.14 

40000 

3.9313E*02 

389.97 

968.07 

45000 

3  0945E  *02 

389  97 

968.07 

60000 

2.4361E  *02 

389.97 

968.07 

65000 

1.9180E*02 

389.97 

968  07 

60000 

1.5103E*02 

389  97 

968.07 

65000 

1.1893E+02 

389.97 

968,07 

70000 

9.3727E*01 

392.25 

970  89 

75000 

7.3990E+01 

394.97 

974.26 

80000 

5.851 1  £>01 

397.69 

977.61 

85000 

4.6350E  +  01 

400.42 

980  95 

90000 

3  6778E*01 

403.14 

984.28 

95000 

2.9232E +01 

405.85 

987.59 

100000 

2.3272E  *01 

408.57 

990.89 

105000 

1.8557E+01 

411.29 

994.18 

110000 

1.4837E+01 

418  38 

1002.72 

115000 

1.1912E+01 

425  98 

1011.79 

120000 

9.6013E400 

433.58 

1020.77 

125000 

7.7688E+00 

441.17 

1029.66 

130000 

6.3094E  *00 

448.76 

1038.48 

135000 

5.1426E400 

456.34 

1047  22 

140000 

4.206  IE  +00 

463  92 

1055.88 

145000 

3.4517E  400 

471  50 

1064.47 

150000 

2.84 18E  400 

4  79  07 

1072.99 

165000 

2.3471E  400 

486.64 

1081.43 

160000 

1 X419E  400 

487.17 

1082.01 

165000 

1.6068E  400 

487.17 

170000 

1.3297E  400 

4«»- 

175000 

1.1000E400 

180000 

9.0836E-01 

185000 

7.485'*r 

BEST 


r  Pi  pH  p  pppy 


Oenuiv  (vlugv/cu  ft) 

2.37G9E  03 
2  0182E  03 
1.7556E  03 
1.49G2E-03 
1.2G73E  03 
1.0CG3E-03 
8.90G9E-04 
7.3320E-04 
5.8728E04 
4.6227E  04 
3.6392E-04 
2.8652E04 
2.2561  E-04 
1.7767E-04 
1.3920E-04 
1.09t3E*04 
8.5710E-05 
6.74  34  E -05 
5.3147E-05 
4.1959E-05 
3.3182E  05 
2.6285E-05 
2.0659E-05 
1.6290E-05 
1.2900E-05 
1.0259E-05 
8.1907E-06 
6.5G50E-06 
5.28I8E-06 
4. 264  8  E -06 
3.4557  E -06 
2.8097 E -06 
2.322  IE -06 
t.92t5E-06 
1.5901  E -06 
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SECTION  3 


BODY  METHODOLOGY 


3.1  INTRODUCTION 

This  section  summarizes  those  methods  available  to  compute  the  aerodynamic 
characteristics  of  bodies.  Those  body  alone  methods  chosen  for  Missile  Datcom, 
surmarized  in  Table  5,  are  the  result  of  the  qualitative  assessment  criteria 
given  in  Table  1.  Rationale  used  in  the  selection  of  these  recommended  techniques, 
as  well  as  available  alternate  methodology  are  discussed  in  the  following  para¬ 
graphs. 

The  most  complete  compendium  Of  methodology  available  is  for  axi sym¬ 
metric  bodies.  Theoretical  analysis  is  greatly  simplified  over  more  complex 
shapes.  The  large  quantity  of  experimental  results  also  allows  for  more 
rapid  and  accurate  analysis.  Since  Max  Munk  first  derived  the  potential  flow 
over  pointed,  slender  bodies  of  revolution,  other  authors  have  developed  more 
sophisticated  techniques  which  allow  accurate  analysis  for  more  realistic 
configuration  shapes.  With  high  speed  computers  being  almost  common-place, 
the  mathematically  complex  methods  are  being  used  more  often  for  design. 

The  capability  exists  to  determine  the  inviscid  normal  force  and  pitching 
moment  characteristics  of  axi symmetric  bodies  across  the  Mach  spectrum.  At 
subsonic  speeds  the  Neumann  Potential  Flow  process  of  simulating  a  body  by  a 
series  of  source-sink  pairs  and  applying  the  condition  of  zero  flow  velocity 
normal  to  the  body  surface  is  well  known.  At  supersonic  speeds,  the  work  of 
Van  Dyke,  Syvertson  and  Dennis,  and  Lavender  in  solving  the  inviscid  Euler 
equations  are  recognized  as  great  strides  in  analytical  prediction  capability. 

The  transonic  speed  regime  has  been  more  elusive,  having  been  solved  just 
recently  by  Klopfer  and  Chaussee  in  a  manner  suitable  for  engineering  design 
purposes  through  numerical  solution  of  the  unsteady  Euler  equations.  Al¬ 
though  theoretical  methods  can  be  identified  at  all  speed  regimes,  use  of  all 
these  methods  is  not  desirable  for  Missile  Datcom.  Neumann  Potential  Flow 
is  a  large  computer  program  which  is  costly,  and  requires  significant  geometry 
detail.  The  solution  of  the  unsteady  Euler  equations  in  the  transonic  Mach 
regime  is  an  iterative  solution,  shown  to  utilize  3900  seconds  computer  time 
per  case  on  a  CDC  7600.  Although  these  methods  are  attractive,  they  are  not 
the  quick  and  inexpensive  methods  desirable  for  Missile  Datcom.  Generation 
of  design  charts  is,  however,  an  acceptable  way  of  using  these  techniques. 
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Available  methods  were  assembled  through  extensive  literature  searches 
of  the  Defense  Technical  Information  Center  (DTIC),  the  National  Aeronautics 
and  Space  Administration  (NASA)  and  the  McDonnell  Douglas  Corporation  libraries 
and  the  archive  journals  and  meeting  papers  of  the  American  Institute  of 
Aeronautics  and  Astronautics  (AIAA).  The  reports  selected  were  categorized  by 
subject  and  screened  according  to  the  required  range  of  geometry  and  flight 
conditions.  This  process  identifies  most  of  the  state-of-the-art  prediction 
methods . 

Of  primary  consideration  is  the  method  type:  theoretical,  semi -empirical , 
or  empirically  derived.  Theoretical  or  semi -empirical  techniques  are  desirable 
because  they  provide  greater  capability  than  empirically  derived  methods.  Since 
method  extrapolation  is  often  required,  the  techniques  with  a  theoretical  base 
can  often  be  extrapolated  without  great  loss  in  accuracy.  Empirical  methods 
using  cruvefits  of  test  results  generally  do  not  follow  any  rationale  physical 
phenomena,  but  are  chosen  for  convenience.  Hence,  theoretical  or  semi -empirical 
techniques  are  preferred  if  the  accuracy,  range  of  applicability,  and  ease  of 
use  of  the  methods  are  acceptable.  Empirical  techniques  are  not  completely 
excluded.  A  proper  mix  of  the  method  types  will  allow  analysis  of  a  large 
range  of  the  selected  requirements. 

3.2  AXIAL  FORCE 

Configuration  stability  and  control  is  of  importance  in  configuration 
design  but  vehicle  drag  (axial  force)  is  the  determining  factor  in  assessing 
the  ability  of  a  configuration  to  deliver  its  payload  to  the  target.  Since 
the  body  is  a  major  component,  its  axial  force  must  be  predicted  accurately. 

The  body  axial  force  is  comprised  of  friction  drag,  nose  pressure/wave 
drag,  interference,  boattail  pressure/wave  drag,  base  drag  and  protuberance 
drag.  The  Mach  number  variation  of  the  major  contributors  are  shown  in 
Figure  6.  The  following  paragraphs  describe  the  recommended  and  available 
methods  for  each  component. 

Skin  Friction  Drag  -  Van  Driest  Method  II  is  recommended  for  skin 
friction  at  all  Mach  numbers.  The  method  requirements  for  skin  friction  are 
based  upon  the  expected  range  of  Reynolds  number  and  Mach  number.  These  ranges 
are  set  by  the  vehicle  flight  profiles,  as  well  as  a  typical  wind  tunnel 
model  size  and  test  conditions. 


30 


The  classic  means  of  evaluating  skin  friction  is  to  use  flat  plate  results 
and  apply  correction  factors  for  three-dimensional  flow  or  compressibility. 

Of  the  many  skin  friction  methods  available,  only  the  Blasius  laminar  flow 
theory  is  considered  theoretical;  the  remaining  methods  are  equations 
which  model  experimental  results.  Many  of  the  smooth  skin  friction  methods 
used  are  summarized  in  Table  6.  The  Mach-Reynolds  number  capability  of  each 
of  these  methods  are  shown  in  Figures  7  and  8.  The  methods  in  Figure  8 
are  presented  assumming  no  compressibility  correction.  The  characteristic 
length  of  one  foot  is  representative  of  a  fin  panel  mean  geometric  chord 
for  a  full  scale  vehicle;  the  ten  foot  characteristic  length  is  representative 
of  a  body.  The  turbulent  Sommer  and  Short  T',  Van  Driest  methods  I  and  II, 
and  Spalding  and  Chi  techniques  are  shown  to  be  candidate  for  further 
review;  Van  Driest  has  been  shown  to  be  superior  to  others  in  its  class  and 
is  the  recommended  automated  approach.  Since  the  Van  Driest  method  is  an 
iterative  technique,  the  Blasius  and  Schoenher  methods  are  chosen  for 
handbook  use.  These  methods  are  given  in  Table  7.  For  Reynolds  numbers 
near  5x10^  per  foot  (comparable  to  high  altitude),  a  laminar  boundary 
layer  will  transition  to  turbulent.  A  suitable  method  for  the  transition 
regime  assumming  an  adiabatic  wall,  is  given  in  Reference  32  and  is  as 
fol lows : 


1.  Determine  the  transition  Reynolds  number 

„  . 22M  6 

ReTRA  '  e  x  10 


2.  Compute  laminar  Reynolds  number 

Reia»,  '  ReTRA  IX0.67/A)1  -^R^-375] 

where,  A  =  0.036  -  .00128M  -  .00072M2 

3.  Compute  turbulent  Reynolds  number 


where  Rg  is  Reynolds  number  based  on  the  characteristic  length  at 
freestream  conditions. 
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4.  Compute  total  friction  of  component 


This  technique  can  be  used  for  each  component  of  the  configuration. 
Heat  transfer  effects  can  be  included  by  using  the  relationship 


1/2 


where  (T„/T  )  is  a  function  of  Mach  number  and  given  in  either  Dunn 

and  Lin,  Reference  20,  or  Van  Driest  and  Boison,  Reference  23.  Note  that 
the  Eckstrom  and  Eckert  methods  handle  transition;  these  methods  are 
recommended  alternates. 

Eaton,  Reference  33,  uses  the  empirical  results  of  Chapman  and  Kester, 
Reference  34.  These  data  are  limited  to  low  supersonic  Mach  numbers  and 
were,  therefore,  excluded. 

Skin  Friction  Compressibility  Correction  -  The  Hoerner  method  is 
recommended.  A  correction  due  to  compressibility  at  Mach  number  is 
commonly  used  to  utilize  incompressible  values  of  skin  friction  at  com¬ 
pressible  speeds.  For  methods  without  heat  transfer,  Hoerner  (Reference 
36)  presents  the  commonly  used  formulas,  which  are  as  follows: 

Laminar  Body  Layer  - 


=  (1 .  +  0.045M2)-4 


Turbulent  Boundary  Layer  - 
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The  turbulent  boundary  layer  relationship  of  Datcom,  Figure  4.2.3.1-68, 
is  also  recommended. 

Surface  Roughness  Correction  -  The  Clutter  method  is  recommended. 

Clutter,  Reference  18,  presents  charts  of  skin  friction  coefficient  as  a 
function  of  Mach  number,  Reynolds  number  and  surface  roughness  (sample  given 
in  Figure  9)  using  a  modified  form  of  the  Prandtl  and  Schlichting  Method 
(Reference  31).  Datcom  presents  a  table,  shown  in  Table  8,  which  equates 
surface  type  to  an  equivalent  sand  roughness.  From  this  admissable 
roughness  height,  a  critical  Reynolds  number  is  determined  from  Datcom 
Figure  4.1.5.1-27.  This  Reynolds  number  is  then  used  to  compute  skin 
friction  using  data  for  a  smooth,  turbulent  flat  plate. 

Friction  is  a  more  significant  effect  at  subsonic  speeds  since  it  is 
a  greater  percentage  of  the  vehicle  drag.  However,  it  should  not  be 
ignored  at  supersonic  speeds.  It  is  recommended  that  the  Clutter  friction 
method  be  used.  It  is  essentially  an  extension  of  the  Van  Driest  II  technique. 
Since  this  method  requires  iteration,  it  is  more  amenable  to  automation.  For 
a  handbook,  it  is  recommended  that  the  Datcom  technique  be  employed  because 
of  its  simplicity. 

Subsonic  Nose  Pressure  Drag  -  The  Hoerner  "form  factor"  correlation  is 
recommended.  At  subsonic  and  transonic  speeds,  the  zero-lift  drag  terms  of 
pressure  and  friction  drag  are  normally  presented  as  a  multiplier  to  the 
friction  coefficient  (Cf),  somstimes  referred  to  as  a  "form  factor". 
Experimental  observation  quantifies  pressure  drag  as  a  direct  function  of 
Cf  and  thickness  to  chord  ratio  (t/c)  as  follows: 


FI  ow 


Pressure  Drag 


2-Dimensional 


3-Dimensi onal 


60  (i)3  cf 

7  (t/c)3  Cf 
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In  addition,  a  correction  is  applied  to  since  the  local  flow  over  the 
surface  is  at  a  higher  speed  than  free-stream  (from  the  Bernoulli  principle). 
Theoretically,  it  can  be  shown  that  the  increase  in  is  as  shown  below: 

FLOW  INCREASE  IN  Cf  DUE  TO  q 

2-Dimensional  2  (t/c)  for  (£)max  =0.3 

1.2  (t/c)  for  (£)max  =  0.5 


3-Dimensional 


1.5  (t/c)1,5 


Hence,  by  summation,  the  following  formulas  are  applicable  at  subsonic 
speeds: 

2- Dimensional : 

=  [1+1 .2  (^-) 1  ’ 5  +  7  (£)3] 

Cf  L  L  ^ref 

3- Dimensional:  ^ 

Although  Datcom  presents  what  appears  to  be  a  radically  different  formula 
for  3-D  flow 


=  [1+0.0025  (^)  +|r7dp] 


Swet 

Sref 


closer  inspection  however  shows  the  results  from  both  methods  are  nearly  the 
same  for  body  fineness  ratios  greater  than  six. 

Transonic  Nose  Pressure/Wave  Drag  -  Chaussee  and  empirical  techniques  are 
recommended.  Since  knowledge  of  the  pressure  distribution  over  a  body  is 
required  to  compute  pressure  drag,  it  is  difficult  to  obtain  such  information 
for  general  body  shapes  without  the  use  of  test  results  or  complex  computer 
programs.  The  "form  factor"  techniques  previously  described  easily  solve  this 
difficul  ty. 

Wave  drag  can  be  handled  theoretically  using  the  variety  of  techniques 
shown  in  Table  9.  Transonic  wave/pressure  drag  is  classically  evaluated 
through  use  of  empirical  correlations  of  cones,  tangent  ogives  and  hemisphere 
shapes  such  as  those  given  in  References  37-42.  Since  nose  bluntness  can  be 
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accounted  for  through  addition  of  the  blunt  nose  cap  drag  to  that  of  a  sharp 
nose,  less  the  increment  of  the  unblunted  cap,  a  series  of  design  charts 
are  available  (MDAC-HB  Reference  32)  for  a  wide  variety  of  cone  and  ogive  shapes 
with  varying  degrees  of  nose  bluntness,  form  Mach  0.5  to  5.0,  and  nose  fineness 
ratios  of  1.5  to  5.0  (Figure  10).  Chaussee,  Reference  43,  performed  the  same 
task  through  numerical  solution  of  the  unsteady  Euler  equations  in  the  transonic 
speed  regime  (Figure  11).  The  results  are  given  in  Table  10  as  a  function  of 
nose  bluntness  ratio,  fineness  ratio  and  Mach  number.  Use  of  the  Chaussee  data 
with  the  empirical  results  for  tangent  ogives  should  accurately  define  the 
transonic  wave/pressure  drag  for  this  common  nose  shape. 

Supersonic  Nose  Wave  Drag  -  The  second-order  shock-expansion  technique 
is  recommended.  At  supersonic  speeds  it  is  shown  in  Reference  14,  and 
Figure  12,  that  modified  Newtonian  theory  applied  to  hemispheres  is  accurate 
to  low  supersonic  speeds.  Since  nose  shape  generality  is  desirable,  the  use 
of  modified  Newtonian  theory  with  the  second-order  shock  expansion  method 
(Reference  44)  appears  to  be  an  excellent  choice.  As  will  be  shown  later, 
second-order  shock  expansion  is  an  excellent  technique  compared  to  other 
theories  for  predicting  normal  force  and  center  of  pressure.  Hence,  this 
one  method  can  be  used  to  evaluate  three  separate  aerodynamic  parameters, 

Cn»  X^p,  C[)y,  that  are  related  by  the  pressure  distribution. 

Other  methods  available  for  ogive  and  cone  nose  shapes  are  shown  in 
Figures  13  and  14.  The  data  presented  in  the  Martin-Marietta  CAMS  program 
(Reference  45)  shown  in  Figure  15  is  very  complete;  wave/pressure  drag  from 
NACA  A52H28  (Reference  46)  has  been  correlated  for  the  3/4-power,  L-D 
Haack,  conical  and  tangent  ogive  shapes.  These  methods  are  selected  as 
alternates  to  the  theoretical  second-order  shock  method  recommended. 

Note  that  these  methods  only  address  axisymmetric  geometries.  There 
is  no  technique  available  to  estimate  wave/pressure  drag  of  arbitrary-shaped 
configurations.  The  only  non-circular  methods  available  are  for  elliptic 
shaped  configurations.  Datcom  provides  a  correction  to  wave  drag  due  to  body 
ellipticity  and  Jorgensen  (Reference  47)  has  derived  an  expression  from 
slender-body  theory  for  constant  elliptic  cross  section  cones.  These  methods 
are  given  in  Figure  16.  Arbi trari ly-shaped  configurations  must  be  analyzed 
through  use  of  a  "paneling"  code  such  as  Supersoni c/Hypersonic  Arbitrary  Body 
Program  (S/HABP),  Reference  48,  whose  inclusion  is  beyond  the  scope  and  intent 
of  Missile  Datcom. 
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Boattail  wave/pressure  drag  -  The  empirical  results  of  Payne  are 

recommended.  The  effect  of  boattail s  or  flares  is  perhaps  the  greatest 

gap  in  methodology.  In  spite  of  the  volume  of  test  data  available  on  the 

subject,  and  reported  by  P.  R.  Payne  (Reference  53),  most  of  the  published 

results  treat  specialized  configurations.  Data  was  available  to  permit 

construction  of  parametric  design  charts  of  conical  and  circular  arc  boat- 

tails,  from  Mach  0.6  to  Mach  1.3  for  the  range  of  parameters  in  Figure  17. 

Since  many  missile  configurations  employ  boattailing,  and  as  much  as  30% 

of  a  configuration  zero-lift  drag  can  be  attributable  to  the  afterbody,  an 

accurate  means  are  required  to  estimate  boattail  pressure  drag  (and  the  effect 

on  base  pressure  drag.  The  Payne  results  are  available  for  pressure/wave 

drag  and  base  drag,  and  their  use  is  recommended  in  this  speed  regime.  Though 

2 

the  range  of  applicability  is  for  (d^/d^  greater  than  0.5,  this  does  cover 
the  great  percentage  of  boattail  designs  considered  for  missile  configurations. 
Since  results  are  not  available  at  subsonic  speeds,  it  seems  appropriate  to 
use  the  Datcom  approach  of  determining  a  separation  point  on  the  afterbody 
and  assume  full  base  drag  over  the  remaining  boattail  area. 

Moore,  Reference  51,  suggests  the  use  of  Wu  and  Aoyoma  theoretical 
technique  (Reference  54),  where 

1/2 

cp  (X)  =  - 1  C*i  -  £)  f  _i  ixi  -Q2  ~  l  -  tji  :  _  /dR\2 

5  /(Y+l)  M»2/3  5-  25  (y+1)  M.2/3  (y+1)  m2  J  Vdx* 

Here,  x,  is  measured  from  the  shoulder  of  the  boattail  and 

1  2 


The  pressure  distribution  so  determined  is  integrated  over  the  boattail  to  deter¬ 
mine  pressure  drag.  This  method  is  limited  to  the  range  from  Mach  0.8  to  1.2. 

To  assure  method  continuity,  the  Payne  results  are  preferred  over  this  theoretical 
technique. 

Although  the  results  presented  by  John  Jack,  Reference  55,  and  shown 
in  Figure  18,  were  performed  using'theoretical  techniques  in  supersonic  flow. 
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the  second-order  shock  expansion  method  has  been  shown  to  given  reasonably 
good  accuracy  and  can  handle  a  wide  variety  of  configurations. 

The  design  chart  methods  for  wave  drag  at  supersonic  speeds  are  not 
recommended  for  automation,  but  they  are  extremely  useful  for  handbook 
computations.  Use  of  the  nose  wave  drag  and  boattail  wave  drag  results 
must  however  be  corrected  for  afterbody  length-,  the  afterbody  pressure 
distribution  is  a  function  of  nose  length.  The  Datcom  design  chart,  given 
in  Figure  19,  handles  this  difficulty. 

Base  Drag  -  Empirical  correlations  are  recommended.  The  methods 
available  to  evaluate  base  drag  are  shown  in  Figure  20.  At  subsonic  speeds, 
the  empirically  derived  Hoerner  result 

,  3  , _ 

C[)b=0 . 029  (Ik)  /VCQf 


is  usee  by  many  because  of  its  simplicity.  Moore  proposed  an  approach  which 
incorporates  an  angle  of  attack  effect,  shown  in  Figure  21.  This  method 
has  the  advantage  of  being  simple  to  apply  and  is  quantitatively  verified 
at  zero  angle  of  attack.  The  Love  (Reference  57)  results  shown  in  Figure 
22  serve  as  an  excellent  data  source.  Aiello,  Reference  58,  presents 
empirically  derived  tables  of  base  drag  due  to  angle  of  attack,  but  are 
limited  in  the  range  1.5  <  M  _<  2.75  and  0  _<  a  <_  30  degrees. 

These  methods  are  for  non-boattai led  bodies.  It  is  recommended  that 
the  angle  of  attack  influence  given  by  Moore  and  Aiello  be  used  to  develop 
a  unified  method  across  the  Mach  regime.  Search  of  the  literature  has  deter¬ 
mined  that  there  are  sufficient  results  available  to  perform  the  following 
method  development  tasks: 

1.  Effect  on  base  drag  due  to  Reynolds  number  for  laminar  boundary 
layers;  if  the  flow  ahead  of  the  base  is  turbulent,  no  effect  on 
base  drag  is  shown, 

2.  Effect  on  base  pressure  due  to  body  fineness  ratio;  for  1/ d 
greater  than  5,  in  supersonic  flow,  the  forebody  effect  can  be 
neglected,  and 

3.  Effect  of  base  pressure  due  to  boattailing  at  supersonic  speeds. 
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Some  methods  to  determine  these  effects  have  been  developed  in  Reference  32. 

It  is  also  recommended  that  the  Datcom  methodology  of  Section  4.6  be  included 
to  account  for  jet  effects. 

Boattail  Wave  Drag  Due  To  Jet  Effects  -  The  method  of  Payne  is  recommended. 
A  body  which  has  a  jet  exhausting  from  the  base  will  experience  an  entrainment 
effect  of  the  jet  exhaust  over  the  aft  body.  The  body  pressure  distribution 
shown  schematically  in  Figure  23  is  altered  depending  upon  the  nozzle  area 
fraction  of  the  base.  The  qualitative  variation  of  boattail  pressure  drag, 
given  in  Figure  24,  illustrates  a  typical  variation  due  to  jet  velocity.  At 
u./u  of  unity,  the  boattail  pressure  distribution  is  unaltered,  hence  this 
point  should  correspond  to  the  true  pressure  drag.  As  jet  velocity  is 
increased,  the  jet  flow  acts  like  a  flow  sink,  thereby  increasing  the  local 
flow  speed  over  the  afterbody  and  increasing  the  boattail  drag.  Since  little 
results  are  available  in  the  region  Uj/uQ=1  to  Mj=1.0,  it  is  assumed  that 
boattail  drag  varies  linearily;  limited  experimental  results  do  reflect  a 
linear  trend.  As  the  jet  velocity  exceeds  sonic  flow,  a  number  of  second- 
order  phenomena  take  place  which  are  configuration  dependent.  Generally, 
exhaust  pluming  will  become  the  dominant  characteristic  as  jet  velocity  is 
further  increased,  acting  like  a  source  and  reducing  boattail  flow  velocity, 
hence,  reducing  boattail  drag.  In  the  extreme,  this  pluming  will  completely 
cancel  the  entrainment  effect  and  cause  a  negative  boattail  drag. 

A  theoretical  method  for  evaluating  this  effect  was  presented  by  Payne, 
Reference  53.  For  a  conical  boattail,  the  change  in  pressure  drag  to  entrain- 


ment  is  ACDgE  =  (K/2)  (ujQ/u0-l)  $ 


where  K  =  entrainment  ratio,  0.0415  at  low  speeds 


UjQ  =  jet  "core"  velocity 
u  =  free-stream  velocity 


o 


and 
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where 


6  = Vl-M2 


e  =  boattail  angle 


K 

r 


=  r  /r 
max  max  s 


Protuberance  drag  -  Recommend  Hoerner  compilation.  Probably  the  most 
comprehensive  compilation  of  drag  data  and  methods  assembled  is  that  by  Sighard 
Hoerner  in  his  book  "Fluid-Dynamic  Drag"  (Reference  36).  The  method  pre¬ 
sented  in  Section  V,  "Drag  of  Surface  Imprefections" ,  is  comprehensive  enough 
to  allow  the  missile  designer  to  quickly  and  accurately  estimate  the  drag 
effect  of  joints,  steps  or  surface  waviness.  A  significant  protuberance 
of  interest  to  the  missile  designer  is  that  due  to  launch  lugs,  or  shoes, 
which  attach  the  missile  to  the  launcher.  Limited  test  results  for  standard 
lugs  were  presented  in  AIAA  Paper  72-969  (Reference  59).  It  is  recommended 
that  selected  data  compilations  from  Hoerner  be  used  to  allow  the  designer 
to  perform  drag  analysis  of  surface  imperfections  or  protuberances  suitable 
for  missile  design  purposes. 

Axial  Force  At  Angle  of  Attack  -  The  methods  for  evaluating  (or 
Cq)  at  angle  of  attack  range  from  an  Allen  and  Perkins  theoretical  equation 
to  empirical  curve  fits.  The  Jorgensen  technique  (Reference  47)  is  simple  to 
apply,  but  is  only  approximate  because  it  assumes  that  varies  as  the  axial 
component  of  dynamic  pressure.  The  CAMS  prediction  code  (Reference  45)  uses 
the  induced  drag  equation 


a 


Cp  =sin  2asin  2  +  nC^ 


REF 


c 


from  Allen  and  Perkins  and  Datcorn  at  subsonic  and  transonic  speeds.  This 
method  can  be  reduced,  through  inspection  to 


Coi  =  CL  sin  a 


Limited  results  from  the  report  "Analysis  of  Datcorn  Methods  as  Applied  to 
Modern  Configurations",  Reference  60,  shows  this  equation  to  be  superior 
to  Jorgensen's  result  at  moderate  angles  of  attack.  Some  authors  have 
instead  chosen  the  relation 
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Cp  =  CL  tan  a 

The  CAMS  technique  at  supersonic  speeds  has  modified  this  method  empirically 
to 

CD  =  CL  tan  ct [1 +K  (0.566  +  0.111M)  (1.1  5-0.075  FA/FN^ 

where  K  is  a  modification  factor  which  is  a  function  of  angle  of  attack. 

Another  empirical  technique  is  also  attributable  to  Martin-Marietta 
(Reference  58).  The  non-linear  axial  force,  from  Mach  0.6  to  1.3,  is  an 
empirical  function  of  Mach  number  and  angle  of  attack. 

CA  =  CA  +  f  (M,  a) 

where  f(M,a)  is  given  in  Figure  25.  The  few  comparisons  with  test  results 
are  observed  to  be  quite  good.  A  good  approximation  is  obtained  at  higher 
supersonic  speeds  by  specifing  CA  to  be  invariant  with  angle  of  attack.  This 
characteristics  is  quite  common  for  bodies. 

The  techniques  for  estimating  body  axial  force  at  angle  of  attack  are 
primarily  empirically  based.  It  is  recommended  that  empirical  results  be 
utilized  at  angles  of  attack  greater  than  approximately  20  degrees.  Use  of 
the  Allen  and  Perkins  and  Jorgensen  results  are  preferred  at  the  lower  angles 
of  attack.  The  method  recomnendations  are  given  in  Table  11. 

3.3  BODY  NORMAL  FORCE  AND  PITCHING  MOMENT 

In  most  methods  the  normal  force  on  a  missile  body  is  assumed  to  be 
composed  of  potential  and  viscous  components.  This  technique  is  used 
extensively,  is  simple  to  apply,  and  models  the  vortex  separation  phenomena 
at  high  angles  of  attack. 

Subsonic  Inviscid  Lift  and  Pitching  Moment  -  Empirical  and  the  Allen 
and  Perkins  methods  are  recommended.  In  1924,  Max  Munk,  Reference  61, 
derived  the  potential  flow  normal  force  of  pointed  slender  bodies  of 
revolution  (based  on  maximum  cross-sectional  area) 

CNp  -  (K2-K-|)  sin  (2a) 

where  (Kg-K, )  is  the  virtual  mass  coefficient  difference  between  transverse 
and  axial  motion  for  ellipsoids  of  revolution  computed  from  Lamb,  Reference 
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62.  Ward,  Reference  63,  examined  Munk's  hypothesis  and  determined  that 
the  resulting  force  vector  should  be  inclined  down-stream  by  an  angle 
a/2  to  the  vertical.  Hence,  normal  force  in  potential  flow  became 

CN  =  (K2"Kl)  sin  (2a)  cos  (a/2) 

This  work  was  extended  by  Allen  and  Perkins,  Reference  64,  65,  and  66, 
when  they  derived  the  potential  normal  force  equation  for  blunt-based 
bodi es , 

CNp  =  (K2-K.|)  ^b  sin  (2a)  cos  (|-) 

The  potential  pitching  moment  was  also  presented,  being 

Cmp  =  (K2"K1 J  [  "~'S^$d~XM)]  Sin  (2a)  C°S  (a/2) 

Using  a  trigonometric  identity,  the  normal  force  becomes 

CNp  =  2  (K2-K.|)  sinacosacos  (a/2) 

Nielsen  (Reference  12)  has  shown  that  for  a  slender  body  with  a  blunt  base, 
slender-body  potential  theory  derives  CNa  as  2.00  per  radian,  based  on  base 
area.  If  CNa  is  known  for  a  body,  then  the  potential  normal  force  and 
pitching  moment  become  sinacosacos  (a/2) 

r  s  r»i  T  ^-ShU-X^)]  sinacosacos  (a/2) 

Hence,  one  is  able  to  approximate  the  potential  normal  force  and  pitching 
moment  of  a  body  through  simple  modification  of  the  slender-body  result. 

Many  empirical  techniques  exist  that  are  curve-fits  of  specific 
data  bases.  The  well-known  empirical  techiques  are  those  of  Martin -Marietta 
Company  (MMC)  and  by  Baker  of  AEDC.  The  MMC  results,  given  in  Reference  58, 
and  the  Baker  methods  of  Reference  74  do  provide  for  excellent  prediction 
as  long  as  the  configuration  being  analyzed  lies  within  the  data  base.  These 
data  may  be  useful  for  analysis  of  selected  configurations  at  transonic 
speeds,  but  the  methods  are  not  recommended  because  of  their  questionable 
extrapolation  capability. 

Since  the  Allen  and  Perkins  potential  normal  force  and  pitching 
moment  results  were  derived  using  the  simplifying  assumption  of  long- 
pointed  forebodies,  many  authors  have  strived  to  obtain  much  better 
results  for  "realistic"  configurations. 
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At  subsonic  speeds,  Moore  (Reference  51)  correlated  data  from  pro¬ 
jectiles  to  evaluate  the  body  normal  force  slope  by  summation  of  the  nose, 
afterbody  and  boattail  components, 

where  (C^)^  is  evaluated  from  empirical  data,  (CNa)^  1s  evaluated  from  the 
theoretical  calculations  of  Wu  and  Aoyoma  (Reference  54),  and  (Cm  )dx  is  an 
empirical  correlation.  One  must  resort  to  empirical  results  at  subsonic  and 
low  transonic  speeds  since  there  are  no  easily  applied  methods  to  compute 
potential  flow  at  these  speeds  other  than  the  Allen  result.  The  well-known 
Neumann  Potential  Flow  (References  90  and  91)  and  Woodward  (Reference  92) 
computer  codes  supply  excellent  results  at  low  angles  of  attack,  but  are 
large,  costly  and  mathematically  complex.  Their  inclusion  is  considered 
outside  the  scope  of  Missile  Datcom. 

Transonic  inviscid  lift  and  pitching  moment.  Klopfer  and  Chaussee 
and  empirical  are  reconmended.  The  method  of  Wu  and  Aoyoma  in  determining 
transonic  theoretical  normal  force  slope  and  aerodynamic  center  is  based  on 
small  perturbation  theory.  Application  to  blunter  nose  shapes  raises 
doubts  as  to  the  methods  applicability.  Chaussee  and  Klopfer  (Reference  93) 
have  numerically  solved  the  three-dimensional  flow  about  axi symmetric 
bodies  and  the  computed  pressure  profile  nearly  reproduced  test  results. 

Since  this  procedure  requires  an  iteration  scheme,  the  technique  is 
extremely  costly  (3900  sec.  computer  time  on  a  CDC  7600)  and  not  desirable 
as  a  module  in  an  automated  Missile  Datcom.  Parametric  results  are  available 
which  provide  theoretical  solutions  in  the  transonic  Mach  regime. 

The  range  of  applicability  of  the  Klopfer  results  are  shown  in  Figure 

26.  The  results  were  curve-fit  with  quadratic  polynomials,  similar  to  the 

one  shown  in  Figure  27.  The  terms  of  the  polynomials  are  Mach  number 

dependent  and  are  presented  in  Tables  12  and  13  for  normal  force  slope  and 

pitching  moment  slope,  respectively.  This  scheme  was  used  because  of  the 

) 

extreme  cost  in  generating  points  for  the  interpolation  tables.  Correlations 
with  test  data  have  been  shown  to  be  quite  good.  The  body  fineness  ratio 
range  investigated  is  smaller  than  desirable,  and  the  range  of  nose  bluntness 
ratios,  from  0.025  to  0.5,  are  limited.  However,  it  is  felt  that  these 
results  are  important  enough  to  warrant  inclusion  since  they  do  cover 
the  lower  body  fineness  ratio  range. 
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As  shown  in  Figure  28,  these  results  do  supplement  the  higher  afterbody 
fineness  ratio  empirical  results  of  Aiello  and  Bateman  (Reference  58)  and 
Krieger  (Reference  49).  Extrapolation  outside  of  the  nose  fineness  ratio 
and  afterbody  length  ranges  covered  should  be  determined.  The  effect  of  nose 
bluntness  is  limited.  It  is  anticipated  that  incremental  nose  bluntness 
effects  can  be  summed  to  sharp-nosed  profiles,  as  was  described  earlier  for 
wave  drag,  to  obtain  the  effects  of  nose  bluntness.  This  was  performed  in 
the  Krieger  compilation  with  success. 

Moore  suggests  computing  pitching  moment  slope  usiing  the  relationship 
CMa  =  -E(CNa)N  (XCP)N+(CNa)A(XCP)A+(%)BT(Xcp)BT 
where  (X^pJ^  and  (X^p)gj  are  evaluated  using  slender-body  theory, 

(YCP)  =  o  -  (VOL)^ 

^CP;N  *N 

ref 

(xCP)BT  -  *  - 

itR^  - 
ref 

and  (X^p)^  is  evaluated  using  the  theoretical  result  from  Wu  and  Aoyoma.  No 
serious  fault  can  be  seen  using  this  technique  at  subsonic  speeds.  At  tran¬ 
sonic  and  supersonic  speeds  more  accurate  results  are  required.  Use  of 
empirical  results,  with  Klopfer  and  Chaussee  theoretical  solutions,  are 
recommended  at  transonic  speeds. 

Inviscid  Lift  and  Pitching  Moment  -  Supersonic  -  Second-order  shock- 
expansion  and  modified  Newtonian  are  recommended.  The  range  of  applicability 
of  the  various  supersonic  theories  are  shown  in  Figure  29.  Since  several 
techniques  are  available,  method  selection  was  based  primarily  on  method 
accuracy  and  the  fundamental  theory  used  in  its  development. 

At  supersonic  speeds,  Van  Dyke,  Reference  94  and  95,  derived  a  second- 
order  solution  for  axial  flow  and  cross-flow.  He  found  that  solutions  for 
the  cross-flow  equations  could  only  be  obtained  for  conical  flow,  so  a 
refinement  of  the  first-order  approximations  were  necessary.  He  proposed 
the  Hybrid  Theory  which  combines  first-order  cross  flow  developed  by  Tsien 
with  second-order  axial  flow.  Other  attempts  to  solve  the  potential  equations 
resulted  in  the  Method  of  Characteristics,  Tangent-Cone  and  the  Shock-Expansion 
Techniques.  The  Method  of  Characteristics  is  the  most  exact  technique  but 
requires  computer  solutions.  Taylor  and  Maccoll,  Reference  70,  formulated 
a  numerical  solution  to  the  shock  wave  equations  proposed  by  Rankine,  where 
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the  position  of  the  Mach  wave  and  the  pressure  over  a  cone  in  supersonic 
flow  were  determined.  Kopal,  Reference  97,  and  the  Ames  Research  Staff 
extended  this  work  in  NACA  Report  1135  (Reference  98)  .  The  Generalized 
Shock  Expansion  method  was  extended  by  Syvertson  and  Dennis,  Reference  44, 
to  cover  a  much  larger  Mach  regime.  Fenter,  Reference  99,  developed  the 
Modified  Second-Order  Shock  Expansion  technique  for  use  with  ogive  cylinder 
configurations.  Newtonian  Impact  Theory  was  developed  by  assuming  that  the 
shock  wave  lies  on  the  body  surface,  attainable  at  very  high  (hypersonic) 

Mach  numbers;  this  method  has  been  modified  to  include  a  pressure  relief 
due  to  flow  centrifugal  force  over  a  curved  surface.  Some  quantitative 

comparisons  with  test  results  are  shown  in  Figure  30. 

In  general,  the  Second-Order  Shock  Expansion  method  performs  better  and 
is  far  superior  to  the  Allen  and  Perkins  result,  since  it  can  adequately 
handle  more  arbitrary  shaped  surfaces.  It  is  available  in  automated  form 
from  many  sources.  Datcom  Figures  4.2.1.1-21  and  4.2.2.1-23,  Figures  31 
and  32,  summarize  these  results  for  both  cones  and  tangent  ogives;  the  test 
results  for  Figure  30  have  been  reported  in  NACA  1328  by  Syvertson  and 
Dennis,  Reference  44.  An  extensive  summary  chart  is  not  readily  available 
for  the  Hybrid  Theory,  because  of  its  limitation  that  the  body  surface  slope 
be  less  than  the  Mach  wave  angle  (see  Figure  33). 

The  second-order  shock  expansion  method  requires  an  attached 
shock,  and  therefore  cannot  handle  the  effects  of  nose  bluntness.  The 
method  proposed  by  Moore,  Reference  131,  is  recommended.  It  uses 
Newtonian  theory  over  the  blunted  cap  and  a  pressure  matching  criteria 
from  perturbation  theory. 

The  effect  of  boattails  or  flares  can  be  handled  using  the  Second-Order 
Shock  Expansion  Theory  for  supersonic  flow,  through  its  accuracy  has  not  been 
thoroughly  examined.  Empirical  boattail  results  at  supersonic  speeds  are 
also  given  in  Datcom  (Figure  4.2 . 1 . 1 -22a) ,  shown  in  Figures  34  and  35; 
their  use  in  a  handbook  is  ideally  suited  with  data  of  Figures  31  and  32. 

At  subsonic  and  transonic  speeds  one  must  resort  to  empirically  derived 
equations  or  charts.  Moore,  Reference  51,  and  Krieger,  Reference  49, 
correlated  a  large  collection  of  test  data  and  then  applied  the  correlation 

<cN«>ur  <cNc,>'BTn-&2]  f  (1/degree) 

cL 
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where  (CNa)'gy  ls  the  boattail  increment  obtained  as  a  polynominal  function  . 
These  results  should  be  adequate  for  missile  design. 

Viscous  Lift  and  Pitching  Moment  -  Allen  and  Perkins  is  recommended. 
Allen  and  Perkins  surmised  that  an  inclined  body  of  revolution  experiences 
a  cross  flow  which  is  the  result  of  viscous  flow  about  the  body. 

The  viscous  cross  flow  acting  on  an  incremental  section  of  the  body  is 
assumed  to  be  produced  from  the  cross-section  drag  coefficient,  Ccjc,  which 
is  a  function  of  Reynolds  Number.  This  incremental  normal  force  is 

AC.  =  Ch  d  sin2  a  AX  /S. 

"vis  c  d 


Integration  of  this  equation  along  the  body,  assuming  Cdc  is  constant, 
results  in 

CNvIS  *  Cdc  ^  sin2c 

and 

cMvis  ■  Cdc  ^  (V^C) 
b  D 

Allowance  must  be  made  for  finite  body  lengths.  Goldstein  (Reference 
73)  determined  a  reduction  factor,  n,  shown  in  Figure  35  for  this  purpose. 
Baker  (Reference  74)  and  Aiello  (Reference  58)  later  modified  n  at  transonic 
speeds.  Hence,  the  Allen  and  Perkins  results  became 


C 

C 


N 

M 


=  (IC,-!^)  Ik  sin  (2a)  cos  (a/2)+r|C(jc  sin2a 


=  (K2-K1)rV'SbU'XM)  ]  Sin  (2a)  cos  (a/2) 


SO 

A 

S 


+  n  cd  [  XM"XC  ]  sin2  a 

c  c  D 


p. 

I  {aA 


Several  methods  are  available  to  determine  n  at  transonic  speeds. 

All  have  been  derived  by  assuming  a  Mach  number  variation  in  Cdc,  several 
of  which  are  shown  in  Figure  37.  A  particular  n  and  Cdc  combination  can¬ 
not  be  recommended  at  transonic  speeds,  but  must  be  determined  quantitatively. 
All  applicable  methods  must  be  assessed  for  accuracy  before  a  choice  is  made. 
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Hill,  Reference  75,  pointed  out  that  the  flow  external  to  the  boundary 
layer  is  potential  in  nature,  so  the  body  model  must  include  the  exterior 
of  the  boundary  layer.  Kelly,  Reference  76,  called  the  use  of  a  constant 
Cdc  alcng  the  body  length  inappropriate;  he  took  the  results  by  Schwabe, 
Reference  77,  and  derived  the  Impulsive  Flow  Analogy,  where  the  develop¬ 
ment  of  cross  flow  on  a  body  is  analogous  to  the  time  dependent  development 
of  cross  flow  force  of  a  cylinder  set  in  motion  from  rest.  The  flowfield 
is  assumed  to  be  developing  along  the  body  length,  being  analogous  to  time. 

In  1966,  Sarpkaya  (Reference  78)  modified  the  data  of  Schwabe  to  remove  the 
inertia  effects.  The  results  are  shown  in  Figure  38.  Perkins  and  Jorgensen 
(Reference  79)  and  Mello  (Reference  80)  did  extensive  studies  of  the  pressure 
and  normal  force  distributions  over  bodies  and  found  that  the  cross  flow  drag 
rose  steadily  and  then  declined  to  a  steady  state  value,  a  trend  similar  to 
the  Schwabe  results.  The  result  of  Kelly's  refinements  became 

/tail 

d  sin^a  dX 

nose 

This  method  should  be  quantitatively  assessed  as  an  alternate  technique. 

In  NOL  TR  73-225  (Reference  89),  Darling  proposed  modifying  the  cross 
flow  drag  along  a  missile  body,  as  shown  in  Figures  39  and  40,  by 
accounting  for  upstream  influence  of  the  base,  the  effect  of  nose  axial 
pressure  gradient,  and  base  influence  at  transonic  speeds.  This  method 
is  also  worth  further  review  because  of  its  more  rigorous  approach. 

Datcom  (Section  4. 2. 1.1)  presents  a  method  by  which  the  viscous 
effect  is  applicable  over  only  that  portion  of  the  body  aft  of  XQ,  the 
position  of  maximum  negative  change  in  cross-section  area.  This  technique, 
though  simple,  is  too  elementary  compared  to  other  available  methods.  The 
accuracy  is  dependent  upon  the  choice  of  XQ,  and  is  subject  to  wide  inter¬ 
pretation.  This  method  was  shown  to  be  less  accurate  than  other  available 
techniques  (Reference  60),  and  is  not  reconmended. 

Since  the  viscous  cross  flow  methods  derived  from  Allen  and  Perkins 
use  the  flow  past  an  infinite  cylinder  to  model  the  effect  of  viscous  normal 
force,  it  is  expected  that  such  models  will  be  most  accurate  in  the  higher 
angle  of  attack  range.  The  time  dependency  noted  by  Kelly  should  serve  as 
an  excellent  analogy  to  the  formation  of  the  vortex  patterns  shown  in 
Fi  gure  41 . 
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Jorgensen,  References  47,  86,  and  87,  has  used  the  method  of  Allen  and 
Perkins  to  develop  a  technique  valid  for  slender  bodies  through  180  degrees 
angle  of  attack.  This  method  adjusts  the  inviscid  and  viscous  components 
by  a  factor  (Cn/CnQ),  which  is  the  local  cross-flow  force  ratioed  to  that  for 
an  infinite  cylinder.  Even  though  most  cross-sectional  shape  experimental 
investigations  have  been  conducted  at  subsonic  speeds,  this  method  is  useful 
at  speeds  where  (Cn/CnQ)  has  been  developed  from  theory.  This  has  been  done 
for  elliptical  cross-section  and  can  be  performed  for  arbi trari ly-shaped 
sections  using  Newtonian  theory.  The  integral  form  of  Jorgensen's  equations 


allow  for  variation  in  body  shape  along  its  length,  making  it  a  far  more  use¬ 
ful  tool  in  missile  design.  Other  methods  for  handling  more  arbitrarily 
shaped  configurations  are  summarized  below. 

Baker  (Reference  74)  has  empirically  modified  the  Jorgensen  pitching 
moment  result  to  bring  the  method  in  line  with  results  of  his  data  base. 

He  has  used  the  equation 
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Figure  42  presents  those  parameters  used  by  Baker  in  his  analysis.  This 
method  is  also  recommended  for  further  consideration  in  the  transonic 
Mach  regime. 

Arbitrary  Shaped  Bodies  -  No  specific  method  is  recommended.  As  can  be 
expected,  there  are  no  simple  techniques  available  to  predict  normal  force 
or  pitching  moment  of  general  body  configurations  at  any  Mach  number.  The 
most  comprehensive  method  available  is  proposed  by  Jorgensen  in  NASA  TN-D- 
6996  and  NASA  TN-D-7228  (References  47  and  86),  where  the  Allen  and  Perkins 
method  is  extended  through  use  of  correction  factors  to  the  potential 
and  viscous  terms  of  and  C  .  The  equations  relevant  to  Jorgensen's 
technique  are  presented  in  Figure  43  as  given  in  Datcom.  Note  that  the 
terms  (C^/Cnc1- r )sb  anc*  (Cjsj/Cnc^ r )^j  correct  the  potential  and  viscous  terms 
due  to  body  ell ipti city.  These  equations  assume  that  the  body  cross-section 
is  uniformly  elliptic  from  nose  to  tail.  Other  results  available  for  (C^/ 

Cl^ci  r )  are  Presentec*  Figure  44  from  Jorgensen.  A  logical  extension  of 
this  method  is  to  permit  the  body  cross-section  shape  to  vary  along  the 
body  length.  The  C|\j  and  Cm  equations  transformed  in  this  manner  were  pre¬ 
viously  shown.  Jorgensen  provided  a  thorough  discussion  of  this  concept 
in  NASA  TR  R-474,  Reference  87,  and  provided  comparisons  with  test  results; 
the  samples  of  which  are  shown  in  Figure  45.  As  anticipated,  fair  agreement 
with  test  was  obtained  at  the  higher  angles  of  attack.  This  is  expected  since 
the  viscous  contribution  is  derived  from  cross-flow  drag  results. 

Another  means  of  using  the  Jorgensen  method  has  been  to  substitute  the 
experimental  cross-flow  drag  coefficient-for  the  particular  shape.  However, 
the  shapes  available  (presented  in  Figure  46)  have  only  been  experimentally 
investigated  at  subsonic  speeds.  Examples  of  other  configurations  tested 
are  shown  in  Figure  47,  indicating  the  probability  of  a  limited  data  base 
for  arbitrary  shaped  slender  bodies.  Since  no  comprehensive  summary  of 
results  have  been  collected,  the  Jorgensen  method  is  constrained  to  circular 
and  elliptical  configurations  at  the  higher  Mach  numbers  by  virtue  of  data 
availability. 

A  second  subsonic  method  used  for  bodies  of  revolution  and  elliptical 
cross-section  bodies  is  similar  to  that  derived  using  the  concept  of  vortex 
lift  of  thin  wings  by  Polhamus,  commonly  referred  to  as  the  Polhamus 
Suction  Analogy  (Reference  102),  and  empirically  extended  in  Datcom.  The 


48 


method,  outlined  in  Figure  48,  has  been  observed  to  fairly  represent  the 
lift  of  elliptical  shaped  bodies  with  a  power  law  planform  shape.  However, 
a  more  unified  approach,  such  as  Jorgensen's  method  is  preferred. 

Williams,  Reference  108,  presented  a  method  based  upon  theory  and  experi¬ 
ment  ot  predict  the  aerodynamics  of  elliptical  lifting  body  geometries, 
similar  to  that  shown  in  Figure  49.  Application  of  the  technique  to  other 
similar  configurations  is  not  known.  The  method  requires  knowledge  of  the 
configuration  pressure  distribution  to  derive  pressure  drag.  The  method 
of  Jorgensen  should  perform  equally  well  yet  have  a  wider  range  of  configura¬ 
tion  applicability. 

A  number  of  lifting  bodies  at  subsonic  speeds  have  been  included  in  the 
Datcom.  The  types  of  configurations  shown  in  Figure  50  have  enabled  the  Datcom 
authors  to  develop  empirically  based  methodology.  These  configurations  are 
applicable  to  missile  design,  and  inclusion  of  the  Datcom  techniques  are  recommended. 
Unfortunately,  methodology  at  supersonic  speeds  is  not  available,  although  programs 
similar  to  S/HABP  are  ideally  suited  to  such  designs. 

In  the  design  environment,  the  S/HABP  code  is  perhaps  the  most  in¬ 
expensive  technique  presently  available.  It  does  however,  require  experience 
in  use  and  methods  choice  to  obtain  good  results,  and  often  requires  "calibra¬ 
tion"  to  a  known,  similar  configuration.  Figure  51  presents  a  pressure  method 
selection  rationale  (Reference  110)  to  reduce  this  "calibration"  time. 

A  promising  technique  still  in  development  is  that  being  performed  by 
Purvis  at  NSWC,  White  Oak  Laboratory.  It  is  designed  to  allow  static  and 
dynamic  aerodynamic  prediction  of  non-axi symmetric  geometries.  This  computer 
program  will  allow  the  aerodynamic  estimation  of  general  body  shapes,  but 
having  the  advantage  of  allowing  the  user  to  build  a  complex  geometry  using 
a  minimal  number  of  inputs.  This  multi-Mach  program  is  considered  a 
significant  advance  in  prediction  capability.  None  the  less,  it  is  not 
recommended  that  such  complex  codes  be  inserted  in  an  automated  Missile 
Datcom.  There  are  many  other  codes  available,  such  as  APAS  from  Rockwell, 

Woodward  from  NASA,  and  PANAIR  (Reference  109)  from  Boeing,  and  referencing  the 
available  methodology  and  data  comparisons  such  as  that  in  Reference  111,  is  a 
significant  benefit  of  Missile  Datcom  which  will  enhance  and  not  detract  from 
its  usefulness. 
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Asymnetric  Forces  -  No  method  is  recommended.  The  cross  flow 
methods  all  strive  to  model  the  flowfield  of  a  body  at  angle  of  attack. 

As  illustrated  in  Figure  41,  experimental  observations  show  two  symmetric 
vortices  forming  at  the  lower  angles  of  attack.  As  the  angle  of  attack 
is  increased  the  pattern  changes  to  a  Bernard- Von  Karman  vortex  street. 

At  much  higher  angles  of  attack  a  new  phenomena  develops,  originally  termed 
"phantom  yaw",  in  which  the  flowfield  shows  unsteadiness  and  the  vortex 
patterns  switch  between  sides  of  the  body  and  exhibit  large  forces  per¬ 
pendicular  to  the  plane  of  the  velocity  vector.  In  theory,  the  time 
integral  of  these  forces  should  be  zero,  but  characteristics  such  as  wind 
tunnel  flow  angularity,  model  surface  imperfections  or  alignment,  cause 
the  vortices  to  favor  one  side  of  the  configuration.  Formation  of  the 
leeside  flow  unsteadiness  disappears  when  the  cross  flow  Mach  number 
exceeds  approximately  0.5  and  its  occurance  is  a  function  of  nose  bluntness, 
nose  cone  semi -apex  angle  and  presence  of  vortex  producers  such  as  strakes. 
Since  this  phenomena  is  due  to  the  boundary  layer  separating  on  the  leeside 
of  the  body,  boundary  layer  blowing  has  also  been  used  to  alleviate  the 
effect.  The  work  of  Wardlaw  (References  81  and  82),  Fleeman  (Reference  83), 
Reding  (Reference  84)  and  Dahlem  (Reference  85)  have  presented  semi- 
empirical  and  empirical  techniques  to  model  this  phenomena. 

For  practical  use  in  the  missile  design  community,  the  following  are 
recommended  for  inclusion: 

o  A  summary  for  the  current  state-of-the-art  in  prediction  methodology 
o  A  summary  of  available  experimental  results 
o  A  concise  description  of  the  phenomena 

o  A  selected  method  which  enables  an  approximation  to  the  magnitude 
of  the  forces/moments  that  are  possible. 

It  is  recommended  that  one  of  the  methods  outlined  above  be  evaluated  to  deter¬ 
mine  its  suitability  for  inclusion.  No  specific  method  is  recommended. 
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TABLE  5  RECOMMENDED  BODY  METHODOLOGY 
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EME1RICAL 


=  Cn  +  Cn  +  Cn  +  Cr 


BASE 


Figure  6. 


Variation  of  Zero-Lift  Drag  Components  due  to  Mach  Number 
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TABLE  6  SUMMARY  OF  SKIN  FRICTION  METHODS 
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HACH  NUNKR  HACH  NUMBER 


CHARACTERISTIC  LENGTH  •  1  FT 


CHARACTERISTIC  LENGTH  -  10  FT. 
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Figure  7.  Applicability  of  Friction  Methods 
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Figure  8.  Applicability  of  Incompressible  Friction  Methods 


TABLE  7  RECOMMENDED  METHODS  FOR  SKIN  FRICTION 


Van  Driest  II  Theory 

(Requires  iteration  for  solution) 

A  -  [1/2  {7-  1)M,/{TW/T00)] 1/3 
B  -  [1  +  .5  (7-l)Mi]/{Tw/T0.|-1 


Schoenher  Theory 
CF  -  (4.13  log  ReCF)*» 


a  -  (2A3  -8)  (B3  +4A 3)‘,/3 
P  -  B(B3  +  4A3)*,/3 

CF  -  [{.242  (sinMa+sin'‘ /?)/[A{Tw/TJ,/3 

Tw 

log  Re  CF  *  .5  (1  +  2w)log  —)  ]3 
Cf  -  [.558/A($in‘*  a  +  sin*  0)CF]/[.558/A 
(sin*1  a  +  sin'1  p)  +  2v^(Tw/T.j‘/a] 


Blasius  Laminar  Theory 

CF  -  1.328  Re',/3 

Cf  -  .664  Re*,/3 
\n 

*/t  -  S.78Re* 

S*/£  -  1.73  R,*,/3 
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Figure  9.  Sample  Presentation  of  Clutter  Results 


TABLE  8  DATCOM  EQUIVALENT  SAND  ROUGHNESS 


Type  of  Surface 

Equivalent  Sand  Roughness 
k  (inches) 

Aerodynamicaiiy  amooth 

0 

Polished  metal  or  wood 

0.02  -  0.08  x  10'3 

Natural  stoat  maul 

0.16  x  10*3 

Smooth  man*  paint,  carafuily  appiiad 

0.25  x  10  3 

Standard  camouflage  paint,  average  application 

0.40  x  10* 3 

Camouflage  paint,  mass- product  ion  x>rty 

1.20  x  10‘3 

Oip^geivenized  matal  aurfaca 

6  x  10’3 

Natural  aurfaca  of  cast  iron 

10  x  10‘3 
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TABLE  9  RECOMMENDED  WAVE/PRESSURE  DRAG  METHODS 
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Figure  10  Sample  Presentation  of  MDAC-HB  Results 


NOSE  FINENESS  RATIO 


TANGENT  OGIVES 


Figure  11.  Transonic  Nose  Wave  Drag  Tangent  Ogives 


PRIORITY  2 


PRIORITY  1 
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TABLE  10  CHAUSSEE  PRESSURE/WAVE  DRAG  AT  TRANSONIC  SPEEDS 
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Figure  12.  Comparison  of  Newtonian  Theory  with  Test  Results 
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Figure  13.  Cone  Pressure/Wave  Drag  Methods  Summary 
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Figure  14.  Tangent  Ogive  Pressure/Wave  Drag  Methods  Summary 
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Figure  15. 


Sample  Results  of  CAMS  Have  Drag  Correlation 
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a)  Datcom  increment  applied  to  circular  sections 


CA  =  a&(2X+l)  +  02ei[3a6X2  +  |  (a2  +  62)  (a  ~b)2  +  ^  a  b 
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8  a  b  ^  +  2  a2  -  b2  log 
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b)  Jorgensen  theoretical  method 


Figure  16.  Wave  Drag  of  Elliptical  Sections 
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Figure  17.  Applicability  of  Payne  Data  Correlations 
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Figure  18.  Jack  Theoretical  Wave  Drag  of  Boattails 
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Figure  19.  Datcom  Charts  for  Determining  Interference  Drag 
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Figure  20.  Base  Drag  Methods  Summary 
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Figure  21.  Moore  Base  Drag  Method 


Figure  22.  Compilation  by  Love,  NACA  3819 
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Figure  23.  Influence  of  Nozzle  Diameter  on  Aft  Body  Pressure  Distributions 
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Figure  24.  Effect  of  Jet  on  Afterbody  Wave  Drag 
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Variation  Of  CA  With  Mach  lfuabir 
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Figure  25. 


Aiello  Empirical  Axial  Force  at  Angle  of  Attack 
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TABLE  11  METHOD  RECOMMENDATIONS  FOR  AXIAL  FORCE  AT  ANGLE  OF  ATTACK 
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Figure  26.  Limitations  of  Klopfer  and  Chaussee  Results 
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Figure  27.  Polynominal  Representation  of  Klopfer  and  Chaussee  Results 


TABLE  12  KLOPFER  AND  CHAUSSEE  COEFFICIENTS,  CN 


COEFFICIENTS  FOR  THE  QUDRATIC  INTERPOLATION 
FORMULA  FOR  THE  NORMAL-FORCE-CURVE  SLOPES 
FOR  THE  ENTIRE  BODY  WITH  A  ONE 
CALIBER  BOATTAIL 


Free  Stream  Mach  Number  Mb 

0.75 

0.90 

0.95 

1.20 

<CNo5o 

1.573 

1.454 

1.534 

2.177 

l  <3W3Vo 

-.625 

-.445 

-.575 

.361 

*  (3CNa/3EN,o 

.296 

.017 

.039 

-.462 

S  (3W3I:a)o 

..107 

.008 

.338 

-.649 

3  (v3Vo 

-1.920 

-1.945 

-1.897 

-.827 

'j2WSKt.2|o 

-.331 

.319 

.476 

-.701 

-.685 

-.821 

1.032 

-.090 

,  "’W’Wo 

.041 

-.814 

-.290 

-.126 

I  '^W’VVo 

-.778 

-1.332 

-.889 

-.712 

j  Ii2Sa/ii:»2'o 

-.278 

-.816 

-.832 

.128 

S  ',2W,E»’E»>o 

-.703 

.240 

.439 

.773 

1  <>J<W»Wo 

-.451 

1.084 

.767 

.562 

-.089 

.297 

-.296 

1.376 

(>2W)rAS?B>o 

.632 

.930 

.038 

.535 

<j2W,5b'o 

.083 

1.465 

1.663 

.104 

vh“'  W<\W  VWw  -  V»-aW 


78 


TABLE  13  KLOPFER  AND  CHAUSSEE  COEFFICIENTS,  C 

m 


COEFFICIENTS  FOR  THE  QUADRATIC  INTERPOLATION 
FORMULA  FOR  THE  PITCHING-MOKENT-CURVE  SLOPE 
FOR  THE  ENTIRE  BODY  WITH  A -ONE 
CALIBER  BOATTAIL 


Free  Stream  Mach  Nunbe: 

* M. 

0.75 

0.90 

0.95 

1.20 

(CW  o 

•  515 

1.591 

1.420 

-1.319 

I 

2.856. 

2.804 

3.639 

.160 

f  (3CMa/3Vo 

-1.280 

-.714 

-1.721 

-7.349 

«  I3W3Ia>o 

2.875 

3.514 

2.122 

4.583 

3  |SWS'b'o 

8.482 

9.478 

9.951 

2.927 

1.544 

-2.282 

-1.431 

4.682 

<s2W3^3r»>o 

11.199 

12.480 

6.177 

-1.214 

«  2  —  — 

-.549 

2.947 

2.823 

-15.080 

E  <3  W3Wo 

3.903 

6.767 

8.013 

10.845 

*  <32W3r»2’o 

-.269 

4.265 

5.151 

7.352 

“  I32W3E»3!:a>o 

6.172 

1.643 

-.611 

-3.682 

1 

14.263 

4.686 

1.562 

-9.186 

°  '^W’V'o 

-.006 

-2.660 

-.135 

-4.087 

7.463 

4.803 

7.536 

-5.228 

l32<W3V>o 

.489 

-7.698 

-9.178 

8.201 

where  ^  =  V(Vmax;  EN  "  Wmax;  LA  =  V(LA>max 
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NOSE  FINENESS  RATIO 


Figure  28.  Transonic  Capability-Ogive-Cylinder  Bodies 
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Configuration 

Neumann  Potential 
Flow 

90,  91 

Subsonic 

Linear  to  Low  a 

Arbitrary  Slender 
and  Blunt  Bodies 

Van  Dyke  Hybrid 

94 

Transonic  — 
Supersonic 
(1.2  <  Mm<  4) 

Linear 

Slender  Bodies 

Second-Order 
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44 
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Method  of 
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100, 
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0° 
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48 
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Figure  29.  Qualitative  Capability  of  Supersonic  Theories 
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Figure  30.  Comparison  of  Supersonic  Methods 
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O  EXPERIMENT 

—  SECONO  OROER  SHOCK  EXPANSION  METWOO 
-—GENEPAH2EO  SHOCK  EXPANSION  METHOO 

—  —  HYBRID  potential  THEORY 

—  -  NEWTONIAN  IMPACT  THEORY 
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Figure  30.  (continued) 


83 


Ctoler  of  pet tsu/f ,  7Af 


0  EXPERIMENT 


1 

L,  J 

Mhb  I  — 

\ 

! 

nr1 

J 

~~  d 

ks 

j  DIMENSIONS  IN  CALIBERS 

1  REFERENCE  LENGTH  * 4 

«'■  SECONOOHOER  SHOCK  EXPANSION  METHCO 
— —  GENE  RAUZEO  SHOCK  EXPANSION  METHOD 

—  — HYanio  potential  theory 

—  -  NEWTONIAN  IMPACT  THEORY 


/ 

/ 

/ 

* 

t 

/ 

/  ^ 

*** 

v 

•*30 

0 

1  1 

/ 

✓ 

✓ 

* 

/ 

*.‘4 

24 

41^/4.61 

1  1 

/ 

/ 

J 

► 

Ho*  5 

05 

*</(.' 72 

1  1 

✓ 

✓ 

/ 

£ 

p 

1 

/ 

A 

<o^6. 

2i 

«n 

'FL. 

02  468  10  02460  10 

Af**rtxxty  f*oen*«  note,  f0 


(b)  Com-,/. -5 


(r)  Cone,  /.  — 3 


(d)  0|nvr,/.-7 


0  2  4  6  8  10 

(f)  Offiv.-,  /.-.l 


Figure  30.  (continued) 
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Figure  31.  Datcom  Design  Charts  Based  on  Second-Order  Shock  Expansion-C^ 
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CONE 


Figure  32.  Datcom  Design  Chart  for  Center  of  Pressure, 
Based  on  Second-Order  Shock  Expansion 
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Figure  33.  Limitations  of  Van  Dyke  Hybrid  Theory 
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Figure  34.  Datcom  Supersonic  Boattail  Correlation 

a 


Figure  35.  Datcom  Center  of  Pressure  of  Boattails 
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Drag  Proportionality  Factor 


Mn  =  M00  sin  a 


o  Baker's  Transonic  Fairing 

n=n*+[(l-n*)/2]  [1+tanh  {  (M-l )  (15/M4)  }] 
n*  -  From  Goldstein's  result 
o  Martin  Marietta  (MMC)  Transonic  Fairing 
n=B0n*+B1 

BQ=  -9.0741  +  31.1111  Mc  -  30.5556  Mc2  +  9.2593  M,.3 

B1  =  10.0741  -  31.1111  Mr  +  30.5556  M  2  -  9.2593  M  3 

^  L  C 

Figure  36.  Cross-Flow  Drag  Proportionality  Factor  Methods 
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USAF  DATCOM  Ref.  (1)  6.  KRIEGER  (MDAC)  (49) 
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Figure  38.  Sarpkaya  Reduction  of  Schwabe's  Results 
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Figure  39.  Variation  of  Cross-Flow  Drag  Used  by  Darling 
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Figure  40.  Correction  Terms  Cross-Flow  Drag  Presented  by  Darling 
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Figure  40.  (continued) 
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Path  of  Symmetrical 
Vortex  Pair - 
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Asymmetric  vortex  formation  and  breakaway 
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Figure  42.  Baker's  Correction  Factor  to  Jorgensen's  Pitching  Moment  Relation 
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When  the  major  axis  (a)  is  perpendicular  to  the  cross-flow  velocity. 


When  the  minor  axis  (b)  is  perpendicular  to  the  cross-flow  velocity. 


0  is  the  angle  of  bank  of  the  body  about  its  longitudinal  axis;  0  =  0 
with  the  major  axis  horizontal,  and  0  =  90°  with  the  minor  axis 
horizontal. 


Figure  43.  Jorgensen  Method  for  Elliptical  Bodies,  as  Given  in  Datcom 
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Figure  44.  Some  Newtonian  Theory  Results  of  Jorgensen 
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Figure  45.  Jorgensen  Integral  Form  Method  Comparison  to  Data 
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Figure  45.  (continued) 
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Figure  45.  (continued) 
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Figure  45.  (continued) 


r  >  kw  \  k  *  .080,  .245, . 


102 


Figure  46.  Other  Shapes  Experimentally  Investigated  at  Subsonic  Speeds 


Votes  B»  ie  the  forebody  of  Bi 
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Figure  47.  Example  Results  for  General  Shaped  Bodies 
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Figure  47.  (continued) 
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Figure  47.  (continued) 
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Figure  48.  Datcom  Version  of  Polhamus  Suction  Analogy  for  Bodies 
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Figure  49.  Geometry  Investigated  by  Williams 
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Figure  50.  Datcom  Subsonic  Lifting  Body  Configurations 
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Figure  51.  Pressure  Method  Selection  Criteria  for  S/HABP 
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SECTION  4 

LIFTING  PANEL  ALONE  METHODOLOGY 


4.1  INTRODUCTION 

Tail  surfaces  on  a  missile  serve  the  primary  purpose  of  stability  and 
control.  In  some  instances,  the  design  requirements  specify  altitude  and/or 
maneuverability  levels  which  cannot  be  met  with  only  a  body-tail  configuration; 
wings  then  are  required  to  increase  the  configuration  normal  force.  As 
described  in  Section  2  of  this  report,  a  significant  number  of  missile  designs 
have  much  lower  aspect  ratio  surfaces  than  those  commonly  used  in  aircraft. 
Missile  aspect  ratios  normally  range  from  one-half  to  four.  Additionally,  the 
most  common  surface  types  are  trapezoidal  and  triangular  panels.  The  following 
paragraphs  describe  the  methodologies  available  for  analysis  for  missile-class 
fin  panels.  Table  14  summarizes  the  recommendations.  Since  Datcom  is 
specifically  oriented  to  aircraft,  and  lifting  surfaces  dominate  the  aerodyna¬ 
mic  characteristics,  it  is  not  surprising  that  the  majority  of  the  lifting 
surface  methods  recommended  are  those  of  Datcom.  The  methods  presented  in 
this  section  assume  a  panel  referenced  axis  system,  analogous  to  the  body- 
axis  system,  where  the  longitudinal  axis  is  oriented  along  the  root  chord. 
Hence,  deflection  or  incidence  of  the  panel  is  analogous  to  angle  of  attack. 

Total  fin  axial  force  is  assumed  to  be  the  sum  of  zero  angle  of  attack 
drag  and  drag  at  angle  of  attack.  The  zero  angle  of  attack  drag  (C/\g)  is 
computed  as  the  sum  of  friction  drag,  pressure  drag  for  sharp  edged  fins, 
leading  edge  drag  due  to  rounding  or  bluntness,  and  trailing  edge  drag.  For 
friction  drag,  methods  which  more  accurately  model  the  two-dimensional  flow 
characteristics  of  fins  are  suggested.  These  are  the  equivalent  skin 
friction  chord  approach  by  Vondrasek  and  the  Barkhem  method  that  approximates 
strip-integration  for  straight-tapered  fins. 

A  good  fin  pressure  drag  method  is  that  of  Moore.  It  may  be 
selected  because  of  its  versatility  as  a  theoretical  method.  However,  this 
technique  requires  a  computer  program  which  will  allow  fin  "paneling"  and 
integration  procedures.  It  may  provide  results  similar  to  that  of  the  Datcom 
which  is  easiest  to  apply. 

The  Datcom  method  is  selected  for  sharp-nosed  pressure  drag  and  leading 
edge  bluntness  drag. 


in 


For  fins  at  angle  of  attack,  the  Datcom  method  is  recommended  due  to 
its  completeness  in  addressing  various  types  of  fins  from  subsonic  to 
supersonic  speeds. 


4.2  AXIAL  FORCE 

Determination  of  fin  drag  for  missile  configurations  has  developed  along 
lines  unique  to  missile  requirements.  Missile  fins  are  typically  designed 
for  efficient  operation  at  high  speeds,  and  have  section  characteristics 
designed  to  have  low  drag,  provide  adequate  control,  and  meet  the  require¬ 
ments  of  manufacturing  relatively  small  panels  in  a  cost  effective  manner. 

Drag  reduction  is  usually  more  important  thart  lift  efficiency.  This  situa¬ 
tion  has  led  to  fin  design  considerably  different  than  typical  airplane  wings 
which  emphasize  lift  performance  in  the  subsonic  and  transonic  regions. 

Missile  fin  planform  characteristics  tend  to  be  low  aspect  ratio  and  simple 
in  shape,  such  as  triangles,  rectangles  and  trapezoids.  Section  characteris¬ 
tics  such  as  twist  and  camber  are  very  seldom  used.  Airfoil  sections  usually 
are  flat  plate,  wedges,  or  biconvex  in  shape  and  have  relatively  sharp 
leading  and  trailing  edges.  Therefore,  the  development  of  drag  methods  for 
typical  missile  fins  has  been  characterized  by  these  traits  and  have  concentra¬ 
ted  on  the  supersonic  speed  regime. 

A  component  build-up  approach  to  fin  drag  is  typically  used  which 
accounts  for  the  following  parameters: 


'FIN 


*^A  +  ^A 
hF  HWave 


+  C, 


+  C, 


+  C, 


Ye 


TE 


where  C.  -  fin  friction  drag 
AF 


C.  -  fin  trail ing-edge  drag 
ATE 


Cfl  -  fin  wave/pressure  drag  C.  -  fin  drag  due  to  angle  of  attack 
MWave  Ra 

Cfl  -  fin  leading-edge  drag 
ALE 

These  component  methodologies  will  be  discussed  for  subsonic,  transonic,  and 
supersonic  speed  regimes. 

Fin  Skin  Friction  -  Recommend  the  use  of  Van  Driest  II.  The  approach 
used  to  calculate  fin  skin  friction  is  similar  to  that  of  bodies.  Flat  plate 
skin  friction  coefficients  are  used  and  corrected  for  three-dimensional 
flow  and  compressibility.  The  Van  Driest  Method  II  is  the  most  comprehensive 
method  in  relation  to  the  requirements,  but  is  an  iterative  technique  and 
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therefore  requires  computer  solution.  The  Schultz-Grunow  and  Schoenher 
methods  are  recommended  for  handbook  application.  For  mixed  laminar  and 
turbulent  flow,  the  transition-criteria  presented  in  Section  3  is  suggested 
when  applicable  to  a  fin. 

Flow  over  a  wing  surface  is  more  two-dimensional  than  over  a  body.  This 
results  in  lower  local  Mach  numbers  and,  therefore,  lower  skin  friction  drag 
relative  to  bodies  at  the  same  conditions.  The  two-dimensional  characteris¬ 
tics  of  fin  flow  has  led  to  the  development  of  methods  to  determine  average 
panel  skin  friction  coefficients.  Table  15  summarizes  three  of  these  methods 
for  fully  turbulent  flow.  The  Butler  method  divides  the  fin  into  strips  using 
the  local  skin  friction  coefficient  based  on  a  mean  local  chord  Reynolds  num¬ 
ber.  Eaton  reduces  the  flat  plate  skin  friction  coefficient  by  an  empirical 
Mach  number  dependent  term  and  a  Reynolds  number  based  on  the  wing  mean  geo¬ 
metric  chord.  Barkhem  proposes  a  method  for  triangular  or  tapered  fins  that 
improves  upon  the  equivalent  chord  approaches.  His  formula  was  developed  by 
correlating  taper  ratio  (x)  and  Reynolds  number  to  values  calculated  by  com¬ 
plete  integration  using  the  Prandtl  flat-plate  turbulent  skin  friction 

5  9 

formula.  For  the  Reynolds  number  range  of  10  to  10  the  formula  agrees 
within  0.2%  with  the  stri p -integrated  values. 

Another  attempt  to  improve  accuracy  is  presented  by  Vondrasek  in  Ref¬ 
erence  118.  She  observed  that  for  a  limited  Reynolds  number  range,  1  og-j oCf -j 
(incompressible)  was  approximately  linear  with  log^Re.  By  simplifying  the 
incompressible  friction  coefficient  equation  to  the  expression 
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she  was  able  to  integrate  to  find  the  average  skin  friction  coefficient  of  a 
straight- tapered  fin  in  terms  of  exposed  root  chord  (Cre)  and  exposed  taper 
ratio  (xE).  The  equivalent  skin  friction  chord  (c*)  was  determined  to  be 
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which  was  shown  to  apply  at  all  Reynolds  numbers  in  subsequent  investigations. 
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The  above  skin  friction  coefficients  apply  to  flat  plates  with  no 
pressure  gradient  and  uniform  temperature.  When  applied  to  fins  of  finite 
thickness,  regions  of  nonuniform  pressure  and  temperature  exist  and  change 
the  coefficients.  A  "farm  factor"  approach  similar  to  that  described  in 
Section  3  is  also  used  for  fins.  Hoerner,  Reference  36,  related  subsonic 
pressure  and  friction  drag  for  slender  airfoil  sections  by  combining  the 
friction  and  pressure  contributions  due  to  thickness.  The  change  in  velocity 
because  of  flow  displacement  results  in  increased  dynamic  pressure  over  the 
fin.  The  differential  is  approximately  2(t/c)  for  maximum  thickness  located 
at  30%  chord.  The  pressure  drag  originates  along  the  afterbody  of  the  fin 
and  was  found  experimentally  to  vary  as  60(t/c)4.  Hoerner's  relationship 
for  a  position  of  maximum  thickness  at  30%  chord  is 

CDo  =1+2  t/c  +  60  (t/c)4 


Examination  of  experimental  data  for  a  maximum  thickness  at  50%  chord 
resulted  in  a  viscous  term  of  1.2  (t/c).  This  has  been  combined  with  the 
above  equation  into 

§ 

CDo/Cf  =  n  +  (2  +  K)  (t/c)  +  60  (t/c)4]  ^ 

ref 


where  K  =  4  [0.3  -  (x/c)m=v]  and  assumed  linearly  valid  between  (x/c)m3v  of 

maX  lilaX 

0.3  to  0.5. 

Datcom  refined  Hoerner  with  the  following  equation. 
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The  L  term  is  identical  to  the  (2+K)  term  from  Hoerner.  A  lifting-surface 
correlation  factor  which  accounts  for  increased  Reynolds  number  length  due  to 
spanwise  flow  has  been  added  from  Reference  119.  The  correction  factor,  R^  ^  , 
is  shown  in  Figure  52.  This  factor  was  empirically  derived  from  wing  data 
having  round-nosed  airfoil  sections.  The  solid  lines  are  used  for  conventional 
straight- tapered  wings,  and  for  outboard  panels  of  non-straight-tapered  wings; 
the  dashed  lines  are  for  use  on  the  inboard  section  of  a  cranked  wing.  For 
transonic  and  supersonic  flow  Datcom  defines  the  friction  drag  as 


Trans oni c 


3  ref 


Supersonic 
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The  wave  drag  term  is  computed  separately.  The  Reynolds  number  charac¬ 
teristic  length  is  the  mean  aerodynamic  chord  of  the  fin.  Compressibility 
functions  for  laminar  and  turbulent  flow  (C^/C^^)  are  the  same  as  described 
for  the  body. 

These  methods  are  reasonably  accurate  techniques  for  determining  fin 
friction  drag.  The  methods  developed  to  increase  accuracy  are  recommended. 

These  are  the  equivalent  skin  friction  chord  (c*)  technique  of  Vondrasek 
and  the  Barkhem  method  that  approximates  a  strip-integration  for  straight- 
tapered  fins.  The  methods  are  valid  for  low  to  moderate  angles  of  attack. 

The  lack  of  methods  at  high  angle  of  attack  is  not  important  because  at  high 
angles  of  attack  skin  friction  is  a  small  contribution  to  drag. 

Fin  Pressure  CAq  -  Subsonic  -  Pressure  drag  arises  from  the  inability  to 
obtain  full  pressure  recovery  of  flow  over  a  fin  due  to  shocks  or  boundary 
layer  displacement.  At  subsonic  speeds  this  drag  is  usually  small  compared 
to  friction  drag,  however,  it  becomes  the  dominant  contributor  at  transonic 
and  supersonic  speeds.  The  empirical  methods  described  for  subsonic  fin 
skin  friction  were  formulated  by  relating  friction  to  pressure  drag.  Datcom 
presents  a  procedure  for  analyzing  non-straight-tapered  planforms  by  treating 
each  component  indivually  then  summing  the  results. 

Fin  Pressure  or  Wave  Caq  -  Transonic  -  The  transonic  regime  is  characterized 
by  regions  of  mixed  subsonic-supersonic  flow  over  the  airfoil.  This  is  shown 
in  Figure  53  and  results  in  a  significant  drag  rise.  The  pressure  drag  contri¬ 
bution  results  from  losses  through  the  shock  on  the  wings.  The  mixed  flow 
pattern  is  not  readily  amenable  to  theoretical  treatment. 

Hoerner  describes  the  application  of  similarity  parameters  to  tran¬ 
sonic  fin  analysis.  Interactions  from  a  variety  of  fin  parameters  such  as 
sweep,  fin  thickness,  maximum  thickness  location,  taper  ratio,  aspect  ratio 
and  leading  edge  characteristics  affect  fin  drag.  Thus,  even  similarity  para¬ 
meters  were  found  difficult  to  correlate  with  data.  Datcom  uses  von  Karman 
similarity  laws  to  show  fin  drag  proportional  to  (t/c)5/3.  Empirical  data 
has  been  correlated  in  Figure  54  as  a  function  of  aspect  ratio  and  thickness- 
to-chord  ratio  for  unswept  wings  and  round-nosed  airfoils.  A  swept  fin  drag 
curve  is  constructed  by  modifying  the  unswept  values  for  peak  drag,  peak  drag 
Mach  number,  and  drag  divergence  Mach  number  (where  3 Cq/3 M=0 .10)  by  the  following 
relationships: 
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Moore  also  uses  the  Datcom  approach  for  transonic  fin  wave  drag.  Other 
techniques  that  are  applicable  to  the  transonic  regime  are  basically  super¬ 
sonic  methods  and  will  be  discussed  as  such. 

Fin  Pressure  or  Wave  C/\q  -  Supersonic  -  Several  methods  exist  for  the 
computation  of  fin  supersonic  pressure  drag.  Most  of  these  are  theoretical 
and  have  been  derived  from  linear  supersonic  theory  as  applied  to  straight- 
tapered  planforms.  The  approximations  of  linear  theory  and  the  relative 
simplicity  of  the  assumed  geometric  shapes  have  allowed  exact  solutions  to 
be  obtained.  The  basic  approach  is  presented  in  Figure  55  which  shows  the 
common  types  of  airfoil  sections  and  the  two  dimensional  and  conical  flow 
regions  established  on  typical  planforms  in  supersonic  flow.  Several  zero 
lift  drag  solutions  have  been  formulated  with  these  models  using  linear 
theory  and  are  summarized  in  Figure  56.  It  should  be  noted  that  these 
results  assume  complete  supersonic  flow  over  the  fin.  In  reality,  the 
Mach  number  variation  on  the  fin  results  in  certain  regions  of  sonic  flow 
in  the  leading  and  trailing  edge  portions  of  the  fin.  When  experimental 
values  are  compared  to  a  typical  theoretical  curve  the  theory  will  over¬ 
predict  drag  in  these  areas  because  the  bow  shock  actually  detaches  resulting 
in  sonic  flow. 

Some  of  the  finite  wing  solutions  are  summarized  in  Figure. 57  through 
59.  The  angle  of  attack  dependency  (a  terms)  are  also  included.  The  equations 
of  Figure  58  are  derived  from  second-order  theory  assuming  an  attached  shock. 

The  pressure  drag  varies  as  the  square  of  the  thickness  ratio  for  a  given  cross- 
sectional  shape.  The  R.A.S.  Data  Sheets  also  provide  drag  equations  for 
three  types  of  airfoils,  Figure  58.  Carafoli  in  Reference  128  used  oblique 
shock  wave  and  expansion  theory  to  compute  the  drag  on  various  supersonic 
profiles;  Figure  gives  the  expressions  (including  angle  of  attack  terms) 
for  several  profiles. 
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The  equations  from  linear  theory  are  adequate  for  determining  trends  of 
fin  wave  drag.  As  already  mentioned,  sonic  flow  regions  over  portions  of  the 
fin  do  result  in  experimental  trends  that  do  not  match  theoretical  "peaks". 
The  Datcom  method  is  obtained  from  linear  theory  for  the  two-dimensional  case. 
For  sharp-nosed  airfoils  the  expressions  for  wave  drag  are, 
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for  subsonic  fin  leading  edge,  (ScotA^^  <  1) 

The  variations  in  fin  thickness  ratio  and  planform  are  accounted  for 
by  defining  the  effective  thickness  ratio.  This  effective  thickness  ratio 
is  defined  by 


(i 


which  is  solved  by  numerical  integration.  The  constant  factor  K  is  defined 
for  specific  sharp-nosed  airfoil  sections  shown  in  Figure  60.  For  fins 
with  variable  thickness  ratios  the  K  factor  is  defined  by  the  average  chord. 
Linear  theory  shows  the  drag  due  to  camber  to  be  exactly  equal  to  the  drag 
due  to  thickness  at  zero  angle  of  attack.  Thus,  fins  with  flat  bottoms  have 
additional  drag  due  to  camber  over  symmetrical  fins  of  the  same  thickness. 

Eaton's  approach  is  similar  to  Datcom.  He  assumes  the  fins  are  con¬ 
stant  thickness  to  chord,  thin  (t/c  <_0.04),  symmetrical,  and  have  a  maxi¬ 
mum  thickness  at  50  percent  chord.  With  these  assumptions  he  uses  an  ex¬ 
pression  similar  to  that  of  Datcom, 

(°0»  •  X  “» 

The  factor  KWw  is  defined  as 


“  (  Ow  (e)  for  X<0.15 

V(C0w  (A\)  for  A <0.1 5 
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The  values  of  (Cj,,),,  are  given  in  Figures  61  and  62,  respectively, 
w  w 

The  R.A.S.  Data  Sheets,  Reference  129,  present  a  series  of  theoretical 
fin  wave  drag  charts.  Figure  63  is  an  example  for  a  double-wedge  airfoil 
with  zero  taper.  The  charts  present  variations  in  maximum  thickness  loca¬ 
tion  and  taper  ratios  for  double-wedge  and  parabolic  airfoil  sections.  As 
previously  discussed,  these  values  will  not  be  accurate  near  the  cusped 
peaks  due  to  sonic  flow  over  portions  of  the  fin. 

Moore,  Reference  131,  suggests  two  techniques  supersonically.  From 
Mach  1.2  to  2.5  he  suggests  using  linear  theory  combined  with  Modified 
Newtonian  to  handle  blunt  leading  edge  flow  regions.  From  Mach  2.5  to  8.0 
he  has  suggested  the  use  of  a  tangent  wedge  strip  theory  that  is  currently 
in  development.  For  the  linear  theory  method,  he  assumes  a  thin  airfoil 
which  is  symmetric,  has  no  camber,  and  is  either  a  biconvex  or  a  modified 
double  wedge  section.  The  procedure  requires  integration  by  numerical 
quadrature  of  the  equation 

8  cb/ 2  rw 
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For  each  flow  region  on  the  panel  the  perturbation  velocity  is  defined  by 
the  individual  surface  slope  of  that  local  element.  This  simplifies  the 
integration  into  a  form  that  is  easily  performed  for  simple  wing  planforms. 

The  airfoil  surface  slopes,  although  assumed  constant  within  an  elemental 
area,  can  vary  between  elements  in  either  the  spanwise  or  chordwise  direction. 

Leadinq/trailinq  edge  bluntness  effects  -  The  methods  described  for  fin 
pressure/wave  drag  applies  to  fins  with  sharp  leading  and  trailing  edges. 

When  fin  edges  are  blunt,  additional  methodology  must  be  applied  in  the  nose 
vicinity  because  the  assumptions  of  perturbation  theory  are  not  valid.  The 
change  in  drag  due  to  bluntness  is  modeled  by  computing  a  leading  edge  drag 

term  that  is  added  to  the  total  drag. 

At  supersonic  speeds  Moore  applies  Modified  Newtonian  Theory  to  blunt 

leading  edges  and  has  derived  the  following  expression  for  leading  edge  axial 
force , 
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The  CpQ  terms  is  the  stagnation  pressure  behind  a  normal  shock  such  that 
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The  term  e  is  the  angle  where  Newtonian  theory  and  perturbation  theory 
match  as  shown  in  Figure  64  and  assumes  that  first  order  perturbation  pres¬ 
sure  coefficients  are  computed  at  e=15°  and  the  flow  is  allowed  to  expand 
to  the  matching  point.  The  leading  edge  radii  are  also  allowed  to  vary 
with  span. 

Datcom  presents  design  charts  for  a  cylindrical  leading-edge  pressure 
correlation  from  Crosthwait,  Reference  132,  which  is  shown  in  Figure  65, 
and  valid  at  subsonic  and  supersonic  speeds.  In  equation  form,  the 
relationship  is 
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Datcom  suggests  using  the  leading  edge  radius  at  the  average  chord  point  for 
variable  radii  fins.  This  approach  has  been  substantiated  over  a  Mach  number 
range  from  0.5  to  8.0  and  for  sweep  angles  from  0  to  75°. 

The  Datcom  or  Moore  methods  are  suggested  for  use  due  to  their  com¬ 
pleteness  and  demonstrated  applicability  in  conjunction  with  sharp-nosed 
pressure  methods. 

Moore  provides  methodology  for  fin  trailing  edge  separation  drag. 

When  the  trailing  edge  is  either  sufficiently  blunt  or  the  surface  slope  suf¬ 
ficiently  large,  the  boundary  layer  will  separate  near  the  trailing  edge. 

This  situation  results  in  increased  drag  due  to  an  equivalent  rear-facing 
step  that  exhibits  two-dimensional  base  drag  characteristics.  Moore  has 
taken  the  experimental  supersonic  blunt  wing  data  of  Chapman,  Reference  134, 
as  a  function  of  Mach  number  and  applied  it  to  fins.  This  incremental  drag 
is  presented  in  Figure  66.  The  curve  was  extrapolated  from  Mach  0.0  to  1.1 
based  on  three-dimensional  base  pressure  trends. 
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Drag  due  to  lift  -  Additional  axial  force  components  can  be  generated 
when  a  fin  is  at  angle  of  attack.  Several  methods  exist  to  determine  drag 
contribution  due  to  lift.  Saffell,  Reference  135,  suggests  that  the  change 
in  drag  due  to  angle  of  attack  be  approximated  by  resolving  a  flat  plate 
normal  to  the  flow  by  the  following  expression: 


The  flat  plate  drag  coefficient  normal  to  flow  is  presented  in  Figure  67. 
This  technique  is  limited  to  0  <_  BA  £  10  and  is  quoted  as  valid  to  180° 
angle  of  attack. 

Datcom  expresses  the  drag  increment  in  terms  of  drag  due  to  lift. 

The  subsonic  drag  due  to  lift  expression  is, 


The  first  term  defines  drag  due  to  lift  resulting  from  induced  drag  and 
viscous  drag.  Induced  drag  depends  on  spanwise  fin  loading  distribution 
and  results  from  the  rearward  rotation  of  the  lift  vector  produced  by  the 
fin  trail ing-vortex  system.  The  viscous  drag  due  to  lift  results  from 
boundary-layer  changes  over  the  fin  resulting  in  an  effective  increase  in 
profile  area.  The  term  is  fin  lift  coefficient  and  the  span-efficiency 
factor,  e,  is 


R^/Aj  +  (1  -  R)ir 


The  factor  R  is  the  leading-edge-suction  parameter,  defined  as  the  ratio  of 
actuel  leading-edge  suction  to  theoretical.  This  factor,  is  shown  in  Figure 
68  as  a  function  of  Reynolds  number  based  on  leading-edge-radius,  Mach 
number,  aspect  ratio  and  leading  edge  sweep  angle.  For  fins  with  sharp 
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leading  edges,  R=0.  These  factors  were  derived  for  fins  with  relatively 
high  aspect  ratios  (greater  than  2.0),  taper  ratios  from  0.0  to  0.713,  and 

leading  edge  sweep  angles  between  19.1°  and  63.4°. 

For  cranked  fins,  represents  the  drag  due  to  lift  resulting  from 

a  breakdown  in  leading-edge  suction  (rounded  edges  are  assumed)  when  sepa¬ 
ration  occurs.  This  term  is  an  empirical  correlation.  The  span  efficiency 
factor,  e,  for  the  cranked  fin  is  determined  by  an  effective  leading-edge- 
suction  parameter  (R1 )  for  the  inner  and  outer  panels  such  that 


R'  =  Rj(rjB)  +  R0(l-V 


Fin  twist  is  accounted  for  in  the  f(e)  term  for  which  charts  are  used.  The 
empirical  factors  v  and  w  are  related  to  induced  drag  such  that 

C  2 

S.  %v+(®%)2  w 

The  w  dependent  term  is  for  a  zero-lift  drag  increment  due  to  twist. 

At  transonic  speeds  Datcom  determines  drag  due  to  lift  of  conventional 
trapezoidal  planforms  of  symmetrical  section  using  transonic  similarity 
parameter  charts,  such  as  Figure  69.  This  approach  is  used  due  to  the  lack 
of  data  at  transonic  speeds.  The  method  covers  the  following  range  of 
parameters,  0<_  AR  tan  £  3.0 

0.5  <  AR  (t/c)1/3  <  2.0 

-4  £  e2/(t/c)2/3  £  2 

At  supersonic  speed  Datcom  suggests  that  fins  be  classified  as  a 
function  of  whether  the  Mach  number  component  normal  to  the  leading  edge  is 
subsonic  or  supersonic.  If  the  fin  has  a  supersonic  leading  edge,  the 
spanwise  pressure  loading  will  be  constant  due  to  two-dimensional  flow. 
Subsonic  leading  edges  can  vary  from  no  suction  to  full  suction  depending 
on  the  subsonic  component  normal  ot  the  leading  edge.  Two-dimensional  flow 
regions  can  be  modeled  by  linear  or  shock-expansion  theory  to  determine  drag 
due  to  lift.  Three-dimensional  flow  regions  have  varying  span  loading  and 
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thus  also  have  a  vortex  drag  contribution.  Datcom  gives  the  two-dimensional 
value  of  drag  due  to  lift  from  linear  theory  as 


The  supersonic  fin  drag  due  to  lift  method  is  given  by  the  relationship 


The  term  in  the  first  bracket  is  shown  in  Figure  70.  The  term  p  is  a  wing 
geometry  parameter  which  is  basically  a  wing  area  fraction.  Comparisons 
with  data  were  made  over  the  following  ranges  of  fin  parameters, 
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The  variation  in  surface  pressure  forces  over  a  fin  at  angle  of  attack 
or  deflection  has  been  determined  from  linear  theory.  The  techniques  shown 
in  Figure  57  include  terms  that  represent  the  influence  of  angle  of  attack 
on  fin  pressure  drag.  It  is  suggested  that  the  results  of  linear  theory  at 
angle  of  attack  be  evaluated  in  conjunction  with  experimental  data  to  deter¬ 
mine  the  accuracy  of  the  method.  Otherwise,  the  Datcom  approach  to  induced 
drag  is  recommended  due  to  its  completeness.  The  Datcom  technique  should 
be  evaluated  quantitatively  for  the  range  of  typical  missile  fin  parameters 
to  determine  the  limits  of  applicability  to  missile  fins  at  angle  of  attack. 

Sufficient  methodologies  exist  to  compute  fin  axial  force  characteris¬ 
tics  for  typical  missile  planforms  and  airfoil  sections.  The  method  des¬ 
cribed  do  not  include  interference  effects.  Such  interference  effects  on 
drag  are  usually  considered  small  and  are  often  neglected. 

4.3  NORMAL  FORCE  AND  PITCHING  MOMENT 

This  section  describes  the  methods  available  to  evaluate  the  linear 
and  nonlinear  normal  force  and  pitching  moment  of  isolated  lifting  surfaces. 
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Many  of  the  techniques  use  design  charts  to  facilitate  calculation;  where 
appropriate  the  relevant  equations  are  presented.  Calculation  of  fin 
alone  or  is  often  suitable  for  hand  calculation. 

Subsonic  Normal  Force  -  At  subsonic  speeds,  the  most  comprehensive  and 
experimentally  verified  straight-tapered  method  used  is  the  lifting  line  theory 
derived  by  Lowry  and  Polhamus,  Reference  136,  and  shown  in  Figure  71.  The 
Lowry-Polhamus  formula  is  considered  valid  over  the  following  ranges; 

M<_  0.8 
t/c<  10% 

An  input  requirement  to  this  method  is  the  determination  of  the  section 
lift-curve-slope.  For  NACA  airfoil  sections,  tables  and  experimental 
results  are  available  in  Datcom  or  "Theory  of  Wing  Sections,"  Reference 
137.  Using  the  Kutta-Joukowski  hypothesis  of  finite  velocity  at  the 
wing  trailing  edge  yields  the  theoretical  method  in  Figure  71,  and  plotted 
in  Figure  72  (bottom)  with  eTE  0f  20  deg  (the  upper  limit  of  the  method). 

Since  this  method  over-predicts  due  to  viscous  effects,  experimental 
results  are  used  to  correct  the  theoretical  result.  The  revised  section 
lift-curve-slope  equation  is  shown  in  Figure  71,  and  the  correction  factor 
is  given  in  Figure  72. 

Although  the  Multhopp  (Reference  138)  lifting  surface  theory  is  highly 
accurate,  it  is  also  difficult  to  apply.  It  has,  however,  been  automated 
for  straight-tapered  surfaces  by  Moore,  Reference  139,  in  the  U.S.  Navy 
Approximate  Aerodynamic  Prediction  Code,  and  is  most  useful  for  low  aspect 
ratio  panels  typical  to  missiles.  Assuming  the  computing  cost  of  this 
technique  is  reasonably  low,  and  the  method  routine  size  is  kept  within 
the  required  computer  field  length  limitations,  its  inclusion  is  suggested  as 
an  excellent  choice  since  it  can  also  be  used  for  non-striaght  tapered  surfaces. 

The  method  more  suitable  to  handbook  application,  and  found  to  be 
fairly  accurate  to  low  aspect  ratio  panels  is  the  Lowry-Polhamus  method. 

It  has  been  extended  by  Spencer,  Reference  140,  to  include  double-delta 
panels  (A<3)  through  calculation  of  an  effective  half-chord  sweep  angle. 

This  area  weighted  sweep  angle  is  used  in  the  Lowry-Polhamus  formula. 

Transonic  Normal  Force  -  Only  empirical  curve-fits 
are  available.  The  most  comprehensive  technique  available  transonical ly  is 
the  empirically  derived  method  outlined  in  Datcom,  Section  4. 1.3. 2.  Datcom 


123 


classifies  wing  panels  as  types  "A"  or  "B".  Thick,  unswept  panels  show  the 
variation  given  in  Figure  73  for  a  type  "A"  wing,  whereas  thin,  low  aspect 
ratio  panels  common  to  missiles  follow  the  trend  for  type  "B".  This  method 
is  a  fairing  of  results  obtained  at  five  distinct  Mach  numbers:  subsonic 
(Mach  0.6),  the  force  break  Mach  number  (Mfb),  Figures  74  and  75  ,  two 
determined  from  Figure  76  (Ma,Mb);  and  supersonically  at  Mach  1.4,  where 

Ma  =  Mfb  +  °'07 
Mb  «  Mfb  +  0.14 

The  lift-curve  slope  at  Mf^  is  determined  from  Figure  71,  and 


(C 


(C 


La) 


=  (1  -  |)  c 
-  (]  -  c 


L«fb 

L«fb 


The  user  is  then  required  to  fair  the  results  according  to  wing  type.  The 
Royal  Aeronautical  Society  data  sheets,  S. 01. 03. 04,  present  transonic  fairings, 
shown  in  Figure  77.  These  charts  present  inviscid  theoretical  flat  plate 
values  of  lift-curve  slope  at  subsonic  and  supersonic  speeds  for  taper  ratios 
from  zero  to  unity.  The  accuracy  of  the  supersonic  values  has  been  shown  to 
be  within  +10%;  accuracy  at  subsonic  speeds  are  not  available.  The  transonic 
fairings  are  suggested  and  not  experimentally  verified.  Aside  from  empirical 
data  results  available  from  Aiello  (Reference  58),  Baker  (Reference  75), 

Nielsen  (Reference  144)  and  Stallings  (Reference  58)  ard  Lamb  (Reference  142) 
for  selected  fin  designs,  no  complete  method  is  available. 

There  are  no  generally  accepted  techniques  for  non-straight  tapered  panels, 
though  Datcom  presents  a  correlation  at  Mach  1.0  for  lift-curve-slope  to  be  used 
as  a  fairing  guideline. 

Supersonic  Normal  Force  -  There  are  more  supersonic  techniques  available 
because  it  is  possible  to  simplify  the  flow  model.  Hoewver,  they  all  reduce 
to  a  form  of  linearized  lifting  surface  theory.  Moore  in  Reference  142  has 
described  in  detail  a  technique  for  straight  tapered,  thin,  uncambered  sur¬ 
faces  for  several  flow  conditions.  General  forms  of  these  relations,  from 
NACA  TN  2114,  Reference  143,  are  given  in  Figure  78.  Five  conditions  for 
which  the  airfoil  pressure  distribution  are  computed  are  shown.  The  applicable 
angle  of  attack  range  is  unknown.  However,  linear  thoery  is  often  valid 
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to  moderate  angles  of  attack.  Datcom  presents  these  theoretical  results  in 
chart  form,  Figure  79.  Regions  where  theoretical  solutions  have  been  obtained 
are  shown  in  Figure  80  and  are  as  follows: 


Region  of  Supersonic  L.  E.  and  T.  E. 

Region  of  Subsonic  L.  E.  and 

Supersonic  T.  E. 

Region  of  tip-root  and  tip- tip 
interactions 

Region  of  Subsonic  L.  E.  and  T.  E. 

Region  of  Ae<0.25  and 

a  <  1.0 

where  a-  0.25  [A(l  +  X )  tanA|_E3 
Region  of  Ag<.  0.25  and  a>_  1.0 


NACA  TN  2114 
NACA  TR  970 


Douglas  SM-13480 


NACA  TR  1050 
NACA  TR  835 


ARC  R&M  2888  and  NACA 
TN  3105 


Datcom  has  combined  this  vast  amount  of  theoretical  work  into  the  easily 
used  design  charts  of  Datcom  Section  4. 1.3. 2,  Figure  4.1.3.2-56.  The  qual¬ 
itative  range  of  applicability  for  the  Datcom  figure  is  shown  in  Figure  8T. 

Note  that  the  R.A.S.  data  sheets  (Figure  77)  will  have  the  same  upper  limit 
of  validity  shown  in  Figure  81,  but  the  ordinate  will  be  half-chord  sweep 
anlge;  the  lowest  limit  is  A^2=0  ^or  aspect  ratios.  Figure  82  graphically 
illustrates  the  conversion  from  half-chord  to  leading  edge  sweep  angles.  It 
is  evident  that  greater  coverage  is  obtained  at  the  lower  aspect  ratios  of 
interest  through  usage  of  the  Datcom  design  charts.  Hence,  the  Datcom  de¬ 
sign  charts  are  recommended  for  handbook  use. 

Thin  airfoil  theory  has  been  assumed  in  development  of  the  Datcom  charts. 
When  the  leading  edge  is  nearly  sonic,  thickness  effects  cause  the  shock  to 
detach  with  a  resultant  loss  of  normal  force.  The  correction  factor  to 
account  for  this  condition  is  presented  in  Figure  83.  One  of  the  unique 
appl i cations  of  the  Datcom  design  charts  occurs  through  use  of  the  reversibi¬ 
lity  theorem  in  supersonic  flow.  Reference  150.  This  theorem  can  be 
summarized  as--"the  normal -force-slope  of  a  panel  in  forward  flight  is  the 
same  as  the  normal -force-slope  of  the  same  panel  in  reverse  flight  at  the 
same  Mach  number."  Hence,  it  is  implied  that  swept  forward  panels  can  be 
handled  using  the  theoretical  results  presented. 
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Figures  84  through  88  present  Datcom  methods  for  non-straight- tapered 
panels  of  interest.  The  results  of  Squire,  Reference  151  and  Figure  88, 
define  a  technique  for  the  little-used  ogee  or  gothic  shapes. 

Normal  Force  at  Angle  of  Attack  -  There  are  four  approaches  used  to 
determine  wing  normal  force  at  angle  of  attack:  (1)  fairing  from  the  linear 
range  to  90  degrees  angle  of  attack,  (2)  computing  non  linear  lift  from  the 
cross-flow  concept,  (3)  variations  of  Polhamus  suction  analogy,  and  (4) 
empirical  curvefits  of  test  results.  Techniques  1,  2  and  4  use  the  fin 
normal  force  at  90  degrees  angle  of  attack  through  empirical  data  correla¬ 
tion;  the  results  given  in  Datcom,  Eaton  (Reference  33)  and  Aiello  (Reference 
58)  are  shown  in  Figure  89.  It  should  be  noted  that  the  results  of  Aiello 
(Figure  89-C)  correlate  better  with  Datcom  (Figure  89-a  subsonically ,  and 
Datcom  Figure  4.1.3. 3-60a  supersonically).  The  Mach  number  idependency  as¬ 
sumption  by  Eaton  is  suspect  since  it  is  not  observed  to  occur  for  planforms  of 
interest  to  Missile  Datcom. 

Since  normal -force -si ope  is  easily  determined,  and  normal  force  in 
normal  flow  is  fairly  well  described,  the  modeling  an  intermediate  angles 
of  attack  become  the  challenge.  The  Datcom  technique  (Approach  1)  uses  the 
relationship 


C..  =  C..  Sin  2a+  C„  Sin  a  | S i n  a 
N  Na - x-  Nact  l 


for  strai ght- tapered  panels  at  subsonic  and  supersonic  speeds.  At  subsonic 
speeds  Cn  is  a  function  of  Cj_MAX;  at  supersonic  speeds  CNaa  is 
a  function  of  the  Mach  detachment  angle  of  attack.  Some  degree  of  empiricism 
cannot  be  avoided  and  this  method  seems  to  be  a  sound  approach.  The  dis¬ 
advantage  is  the  lack  of  methods  in  the  range  0.6<M<1.4.  The  alternate  cross- 
flow  approach 


CN  "  CN 


„  •  2 

Sin  2  a  +  Cdc  sin  a 


does  not  adequately  handle  such  phenomena  as  fin  stall  and  shock  detachment, 
since  the  data  base  deriving  C(jc  is  based  on  only  normal  flow,  i.e.,  panel 
stalled  or  shock  detached  at  all  speeds.  Empirically  derived  curves  of  CL 
versus  angle  of  attack  are  given  in  Datcom,  and  presented  in  Figures  90 
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and  91,  for  double-delta  and  gothic  or  ogee  planforms  at  subsonic  speeds. 
Results  are  also  available  from  NASA  TN  D-5661 ,  Reference  152;  one  chart 
from  this  report  is  shown  in  Figure  92.  These  results  were  generated  by 
employing  the  modified  Multhopp  approach  on  a  series  of  panels  to  generate 
design  charts  at  subsonic  speeds.  Limitations  on  these  design  charts  are 
as  follows: 


Mach  Number 

Trail ing-Edge  Sweep 

Inboard  Panel  Sweep 

Prandtl -G1 auert  Compressibility  Rule 

Zero 

Outboard  Panel  Sweep  Taper  Ratio 

55°  to  85° 

O 

O 

LO 

0.1 

to 

0.5 

0°  to  85° 

o 

o 

LO 

0.1 

to 

0.5 

0°  to  85° 

O 

C\J 

0.1 

to 

0.5 

This  report  also  includes  results  for  the  following  aerodynamic  parameters: 

CL ’  Xac’  %’  Cmq ’  and  CLq' 

Experimental  results  for  sharp  delta  wings  at  low  angles  of  attack 
show  that  the  flow  separates  from  the  leading  edge  and  rolls  up  into  two 
spiral  vortex  sheets.  Acceleration  of  the  flow  in  this  manner  produces  lift. 
The  results  of  flow  separation  at  the  leading  edge  is  called  potential 
lift,  and  the  effect  of  the  spiral  vortex  is  termed  vortex  lift.  These 
effects  are  shown  in  Figure  93.  It  is  assumed  that  the  leading  edge 
suction  force  produced  in  potential  flow  is  the  same  as  the  pressures 
required  to  maintain  flow  equilibrium  and  attached  flow.  The  resulting 
force  is  then  equal  to  the  theoretical  leading  edge  suction  force  (C^) 
and  is  perpendicular  to  the  leading  edge.  This  Polhamus  Suction  Analogy 
concept  (References  102,  153,  and  154)  has  received  much  attention  since 
it  attempts  the  invsicid  plus  viscous  modeling.  The  potential -flow 
normal  force  and  lift  are 


*-L,p  =  Kp  sin  a  cos2a 
CL,p  =  cN,p  cos  « 


and  the  viscous  (suction  force)  is 

Cl,v  =  Ky  cos  a 
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and  Kp  and  Ky  are  derived  from  the  modified  Multhopp  theory  (Reference  155) 
for  triangular  panels.  This  theory  was  found  to  agree  well  with  test  data 
at  very  low  speeds  to  25  degrees  angle  of  attack  for  aspect  ratio  panels  less 
than  two.  The  method  was  also  extended  by  Lamar  (Reference  156),  among 
others,  to  deduce  leading-edge  and  side-edge  vortex  lift  factors  for  appli¬ 
cation  to  straight  tapered  panels.  The  normal  force  (lift)  of  a  panel  is 
represented  as 

CL  *  *  sin  a  cos2  a+ltvlp  sin2  a  cos  a  +  K  sin2  a  cos  a 
*  vjj-c  v,se 


K  sin  a  cos  a 


ref 


K  sin^  a  — 
v,le  c 


'll 

ref 


+  K  sin^  a  — 
v,  se  c 


ref 


where  Kv,|^  is  that  given  above  for  triangular  panels  and  for  panels  with 
subsonic  leading  edges,  and  sonic  trailing  edges 

m  *  B  cot  A 

V  m  'bk cot  A 

Sref  fct’  +  (1  -  m)/2] 


For  rectangular  panels 


K 


v,  se 


Charts  at  subsonic  speeds  were  given  by  Polhamus  in  Reference  157  for  selected 
fins . 

Lecat  and  Rietschlin  of  Grumman  Aerospace  developed  "Goniometric  Aero¬ 
dynamics,"  Reference  158,  and  achieved  good  correlations  through  modification 
of  the  Polhamus  Kp  and  Kv  formulations.  A  planform  shape  parameter,  p,  given 
in  Figure  94,  simplified  the  formulation  for  general  shaped  surfaces  to 


4t 


tan  $  +  j tan2  +  402 

V  Pm 

Kv-[K,-Kj^At]yi+  tan2  ^ 

JS_ 


Kv  * - cos  v 

cos  <t> 


sin  <t>  *  tan  y  [tan  (y  +  4>)  ] 
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Subsoni c 


2p*  tan  0  =  tan  f 


M 


=  P* 


Transom'  c 


Supersoni c 


_El _ 

tan  0m 

2  tan  0 


tan  =  2py(tan  0  -  tan  0M) 


*  p#  (2  tan  0/tan  0M) 
tan  0n  *  2  tan  ^(1  -  tan  0/tan  0«) 


The  authors  claim  accuracy  with  this  technique,  but  use  of  this  method  for 
a  wide  range  of  shapes  should  be  evaluated  before  it  is  chosen  for  imple¬ 
mentation,  because  it  is  not  well  known  within  the  industry.  Oberkampf, 
Reference  159,  extended  the  vortex  concept  to  60  degrees  angle  of  attack 
using  the  following  formulation 

CN(°0  =  [Kp  sin  «  cos  a  +  Kv  sin2a  ] -^coT^-g- 
for  0<a<as 

CN(«>  *  CCN(»s)  +*§  cN(„_„s, 


where , 


for  as<a  <  2as 

CNC»)  -  (C  *  AA  CN(aj), 

for  2as<a 

£  =  0.7  +  0.3M2  (M  < 1 ) 


a. 


 1  .3 


VTTa  (1  .5  +\) 


Kp 


1  .45  A  -  01 7  A 
Kv  =  tt(1  +  — y  +  a*  )  for  A  <  3 
CN  CN  («.A.X,|M  =  0 

Bradley,  Reference  160  and  Figure  95,  has  devised  a  method  where  the 
leading  edge  suction  force  is  computed  across  the  span.  Hence,  the  vortex 
lift  is 

C.  =  [E(Ctn 


Cos  A  n 


)  +  Cy  -H(CTm  tanAj]  cos  a 
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where  Ct„  and  C 
•  n  v 


are  computed  from  lifting  surface  theory. 


Therefore , 


and 


C.  _  Kp  si n  a  Cos +  (Ky  +  K..  )  sin^a  cos  a 

d  LE  VTiP 

Y  r  dsin^  a  (^Yn  ) 

cos  An 


K  =  5"s?n^a  (C^  -£CTn  tan  An} 

Tip 


This  technique  has  shown  remarkable  accuracy  to  30  degrees  angle  of  attack. 
Like  the  Grumman  method,  these  methods  should  be  verified  for  missile-type 
lifting  surfaces  and  general  accuracy. 

There  is  a  large  body  of  test  results  (References  161  through  165,  for 
example)  which  have  explored  the  use  of  wraparound  fins.  For  those  launchers 
which  limit  the  span  of  fins  attached  to  the  body,  this  type  of  fin  is  ideal. 
Experience  has  shown  that  at  low  angles  of  attack  the  aerodynamic  charac¬ 
teristics  are  the  same  as  that  for  an  equivalent  flat  panel.  There  are  no 
theoretical  methods  available  to  analyze  wraparound  fins  at  angle  of  attack. 
Reference  to  experimental  results  are  required. 

There  is  a  large  body  of  panel  alone  empirical  results  available,  due 
primarily  to  the  efforts  of  Baker  (Reference  74),  Nielsen  (Reference  141), 
Hill  (Reference  166),  and  Stallings  and  Lamb  (Reference  142).  The  Stallings 
and  Lamb  paper  presented  at  the  AIAA  19th  Aerospace  Sciences  Meeting  in 
January  1981,  is  a  good  summary  of  the  available  results.  The  available 
data  cover  the  following  ranges  of  conditions: 

Taper  Rati o  0 .0  to  1 .0 

Aspect  Ratio  0.5  to  2.0 

Mach  Number  0.6  to  3.0 

Angle  of  Attack  0  to  60  deg 

The  data  consists  of  normal  force,  longitudinal  center  of  pressure,  and 
lateral  center  of  pressure.  Although  these  results  do  not  quite  cover  the 
range  of  paneT  aspect  ratios  or  Mach  numbers  required,  they  do  cover  the 
transonic  Mach  regime,  which  has  been  shown  to  be  a  difficult  analysis 
area.  It  is  recommended  that  these  results  be  considered  for  use  in 
Missile  Datcom. 
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Pitching  Moment  -  Pitching  moment  methods  exist  for  each  of  the  Mach  regimes 
specified  for  normal  force  above.  The  Datcom  contains  an  extensive  compendium 
of  generally  accepted  methodology  which  utilizes  the  procedures  and  theories  applied 
for  normal  force.  It  is  appropriate  to  select  the  pitching  moment  technique 
which  has  been  developed  using  the  same  derivation  assumptions  for  normal  force. 

The  method  sources  for  the  recommended  normal  force  methods  apply  for 
pitching  moment  as  well. 

Pitching  Moment  at  Angle  of  Attack  -  The  panel  aerodynamic  center  can  be 
evaluated  theoretically  at  subsonic  and  supersonic  speeds.  Datcom  has  corre¬ 
lated  section  center  of  pressure  for  the  double  wedge,  modified  double  wedge 
and  circular-arc  sections  due  to  angle  of  attack,  as  shown  in  Figure  96. 

In  addition,  Datcom  design  charts  derived  through  integration  of  the  theoreti¬ 
cal  pressure  distribution  on  the  panel  are  shown  in  Figure  97.  An  empirical 
method  is  also  presented  in  Figure  97  at  transonic  speeds  through  correlation 
of  aspect  ration  and  thickness-to-chord-ratio.  The  Datcom  methods  are  of 
sufficient  detail  and  it  is  recommended  that  they  be  employed  for  handbook 
use.  If  selection  of  the  Multhopp  (subsonic)  or  pressure  region  (supersonic) 
theories  are  chosen  for  automation,  it  is  recommended  that  the  results  be 
integrated  to  obtain  pitching  moment  or  center  of  pressure.  Use  of  the 
empirical  results  available  are  also  recommended  for  determining  the  effects 
of  angle  of  attack  in  the  transonic  Mach  regime.  It  is  recommended  that 
results  of  the  Polhamus  Suction  Analogy  method  be  quantified  for  low  aspect 
ratio  panels. 
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TABLE  15  FIN  SKIN  FRICTION  METHODS 
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Figure  52.  LIFTING-SURFACE  CORRELATION  FACTOR  FOR  SUBSONIC  MINIMUM  DRAG 


Figure  53.  MIXED  FLOW  REGIONS  AT  TRANSONIC  SPEEDS 
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Figure  54.  TRANSONIC  ZERO-LIFT  WING  WAVE  DRAG  FOR  UNSWEPT  WINGS  AND 
ROUND-NOSE  AIRFOILS 
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Figure  55.  SUPERSONIC  FLOW  REGIONS  OVER  FINS 
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Figure  56.  SUMMARY  OF  FINITE  WING  SOLUTIONS 
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Figure  57.  Aerodynamic  Characteristics  for  Supersonic  Airfoils  -  Bonny 
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Figure  58.  Supersonic  Airfoil  Section  Data  -  R.A.S.  Data  Sheets 
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Figure  59.  Supersonic  Airfoil  Characteristics  -  Carafoli 
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SHARP-NOSED  AIRFOILS 


Figure  60.  Sharp-Nosed  Airfoil  Thickness  Factor 
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Figure  61.  Wave  Drag  of  Zero  Taper  Rotio  Lifting 
Surfaces  (VVings  or  Tails)  ( Xw  *  0.15) 
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Figure  62.  Wove  Drag  of  Lifting  Surfaces  (Wings  or  Tails) 
Whose  Toper  Ratio  is  Greater  Than  0*15 
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(I)  DOUBLE  WEDCE  (A-0.  m-03) 


Figure  63.  Fin  Theoretical  Wave  Drag  -  R.A.S.  Data  Sheets 
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Fi  gure  64.  Combined  Newtonian  and  perturbation  theory  for  a  blunt  leading  edge. 
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Fi  gure  65.  CORRELATION  OF  CYLINDRICAL  LEADING-EDGE  PRESSURE  DRAG 
COEFFICIENTS 
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Figure  65.  (continued) 


.8  1  2  4 

MACH  NUMBER 


qiErnxrnxxi 


4(b)  LEADING-EDGE  SWEEPS  FROM  55  TO  80  DEGREES; 

ztm 


bw 


(deg) 


- 


70- 


146 


147 


0.0  1.0  2.0  3.0  4.0 


MACH  NUMBER 


Fi  gure  67  .  Drag  Coefficient  for  a  Flat  Plate  Xormai  to  the  Flow 
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Figure  68.  LEADING-EDGE  SUCTION  PARAMETER  AT  SUBSONIC  SPEEDS.  M  s  0.8 
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Fi  gure  69.  TRANSONIC  DRAG  DUE  TO  LIFT 
(»)  A  TAN  *  0 
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FIGURE  4.1.3.2-49  SUBSONIC  WING  LIFT-CURVE  SLOPE 
Theoretical  section  lift  curve  slope: 

«  6.28  +  4.7  t/c  [  1  +  .00375  (per  rad) 

Section  lift-curve  slope  corrected  for  viscous  effects: 
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Figure  71.  Subsonic  Wing  Lift-Curve  Slope 
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TYPE  “A”  WINGS  TYPE  “B”  WINGS 


Figure  73.  Wing  Type  Classification 
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Fi  gure  74 .  TRANSONIC  FORCE-BREAK  MACH  NUMBER  FOR  ZERO  SWEEP 
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Figure  75 .  TRANSONIC  SWEEP  CORRECTION  FOR  FORCE-BREAK  MACH  NUMBER 
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TRANSONIC  SPEEDS 
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CORRECTION  TO  LIFT-CURVE  SLOPE  AT  FORCE-BREAK  MACH  NUMBER 


b)  CHART  FOR  DETERMINING  LIFT-CURVE  SLOPE  AT  M, 
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c )  CHART  FOR  DETERMINING  LIFT-CURVE  SLOPE  AT  Mb 

Figure  76.  Datcom  Transonic  Fairing  Parameter 


156 


157 


Figure  77.  Panel  Li ft-Curve-Slope-R.A.S .  Data  Sheets 
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Figure  78.  NACA  TN  2114  Pressure  Regions 
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Fi  gure  79 .  WING  SUPERSONIC  NORMAL-FORCE-CURVE  SLOPE 


Figure  80.  Datcom  Normal  Force  Slope  Design  Chart  Theoretical  Sources 
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LEADING  EDGE  SWEEP  ANGLE,  DEG 


Figure  81.  Applicability  of  Datcom  Charts 
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Figure  82.  Comparison  of  Half-Chord  and  Leading  Edge  Sweep  Angles  for 
Straight-Tapered  Panels 


160 


FLAT  PLATE 


1.0 

.9 

.8 


N, 


(cnJ 


.7 

.6 

.5 


theory 


0 


TANA 


LE 


TAN  Ale 


Figure  83.  SUPERSONIC  WING  LIFT-CURVE-SLOPE  CORRECTION  FACTOR  FOR 
SONIC-LEADING-EDGE  REGION 
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Figure  4.1 .3.2-63  may  dto  ba  applied  with  equal  facility  to  obtain  normal -tor oa-curve  slop#  of  tbagto*. 
Uing  (CNa/A)f  from  Figure  4.1 .3.2-63 


Figure  84.  Datcom  Method  for  Gloved  Panels 
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Fi  qure  85 .  LIFT-INTERFERENCE  FACTOR  FOR  NORMAL-FORCE-CURVE 
SLOPE  AT  SUPERSONIC  SPEEDS 
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Figure  86.  LEADING-EDGE-EFFECT  FACTORS  (CLE)g  AND 

FOR  NORMAL-FORCE-CURVE  SLOPE  AT  SUPERSONIC  SPEEDS 
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Figure  87.  WING  SUPERSONIC  NORMAL-FORCE-CURVE  SLOPE,  X  *  0 
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Figure  88.  CORRELATION  OF  NORMAL-FORCE-CURVE  SLOPE  AT  SUPERSONIC  SPEEDS 
FOR  GOTHIC  AND  OGEE  PLANFORMS  HAVING  SHARP-NOSED  AIRFOILS 


166 


a)  Subsonic  Lift  Variation  with  Wing  Aspect  Ratio  at  90°  Angle  of 
Attack,  Datcom. 


b)  Eaton  Fin  Cross-Flow  Drag 


Figure  89.  Methods  for  Fin  Cross-Flow  Drag 
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Figure  89  (Continued) 
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Fi  gure  90.  PREDICTION  OF  NONUNEAR  LIFT  OF  DOUBLE-DELTA  PLANFORMS  AT 
SUBSONIC  SPEEDS 
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SUBSONIC  SPEEDS 


Fi  gure  91 .  CORRELATION  OF  LIFT  CURVES  OF  GOTHIC  AND  OGEE  PLANFORMS 
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Figure  92.  Subsonic  Lift-Curve-Slope  for  Double-Delta  Panels 


171 


Contribution  of  various  affects. 


Figure  93.  Polhamus  Suction  Analogy 
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Figure  94.  Panel  Shape  Parameter  for  "Goniometric  Aerodynamics" 
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Figure  95.  Results  of  Method  From  GD/Convair  -  Bradley 
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Fi  gure  96. 


Effect  of  Thickness  and  Angle  of  Attack  on 
Center  of  Pressure  Location 
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SUBSONIC  AND  SUPERSONIC  SPEEDS 
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Figure  97.  Datcom  Charts  for  Panel  Aerodynamic  Center 


SECTION  5 

COMPONENT  INTERFERENCE  METHODOLOGY 


5.1  INTRODUCTION 

The  component  interference  methodology  is  a  critical  concept  in  the 
build-up  approach  selected  for  Missile  Datcom.  These  mutual  interference 
effects  are  normally  described  as  two  separate  phenomena,  carryover  inter¬ 
ference  and  vortex  interference.  The  methods  available  to  evaluate  each  are 
discussed  in  the  following  sections. 

5.2  CARRY-OVER  INTERFERENCE 

This  section  presents  a  summary  of  the  methodologies  for  predicting 
non-vortical  inteference  effects  on  missile  configurations.  The  inter¬ 
ference  prediction  techniques  available  include  body-fin  and  panel-panel 
interference  for  deflected  and  undeflected  straight-tapered  fins.  Method 
development  has  emphasized  interference  effects  on  lift  and  center-of -pressure. 
However,  some  empirical  methods  model  the  effects  of  fin/afterbody  combinations 
on  axial  force. 

The  development  of  prediction  techniques  for  the  interference  among 
components  of  high  speed  missile  designs  began  in  the  1950's  due  to  use  of 
low  aspect  ratio  fins  on  relatively  large  diameter  bodies.  The  work  of 
Pitts,  Nielsen  and  Kaattari  (PNK)  published  in  1957,  Reference  167,  presented 
an  approach  for  calculating  body-fin  interference  based  on  slender-body  theory. 
The  approach  adopted  for  lift  interference  relates  fin  alone  lift  to  combined 
body-fin  lift  by  the  factor  such  that 

KC  =  (ClciW(Clc^W 

The  interference  factor  is  handled  by  separately  considering  the  interactions 
of  the  body  and  fin  as  a  function  of  angle  of  attack  and  fin  deflection.  The 
individual  interference  effects  are  thus  defined  as 

K^I(B)  -  fin  interference  due  to  body,  variable  angle  of  attack 

Kg(w)  -  body  interference  due  to  fin,  variable  angle  of  attack 

^W(B)  "  ^in  1nter^erence  due  to  body,  variable  fi-n  incidence 

kg(w)  -  body  interference  due  to  wing,  variable  fin  incidence 

Figure  98  illustrates  the  influence  regions  used  by  PNK  for  a  typical  missile 
concept. 

Body/Mi nq  Interference  (no  deflection)  -  The  angle  of  attack  interference 
of  the  body  on  the  wing  is  shown  by  slender  body  theory  to  be 


177 


-Kirc#»  •" 


HT 


tin 


This  factor  is  applicable  at  subsonic  and  supersonic  speeds,  small  angle  of 
attack  (linear  regime),  and  low  aspect  ratio  fins  of  straight-tapered  planform 
with  non-swept  trailing  edges.  To  compute  the  fin  influence  on  the  body  due 
to  angle  of  attack,  PNK  derived  a  slender-body  theory  value  of 
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The  interference  factors  K^g^  and  Kg^  are  shown  in  Figure  99  as  determined 
by  slender-body  theory.  It  is  interesting  to  note  that  whenever  slender-body 
theory  is  used,  the  sum  of  these  two  factors  can  be  expressed  in  terms  of 
geometry  by  the  relationship  (Reference  171), 

’Sf(B)  +  KB(W)  “  (£  + 

When  slender-body  values  were  compared  to  the  results  from  conical  flow 
theory  in  supersonic  flow  for  triangular,  rectangular  and  trapezoidal  plan- 
forms,  differences  were  observed  for  certain  Mach  number/planform  combinations. 
The  conical  lifting  solutions  were  developed  by  assumming  the  planar  model 
shown  in  Figure  100;  the  lift  on  the  body  was  calculated  by  integrating  pres¬ 
sures  due  to  the  half-fin  over  the  influenced  region  of  the  infinite  body. 

The  resulting  interference  factor  equation  for  a  supersonic  leading  edge 
was  found  to  be 
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For  subsonic  leading  edges  the  equation  is 
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where 
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When  the  body  is  cut  at  the  fin  trailing  edge,  the  solution  yields  the  following 
equations ,  K.lw>[(i(Ctm)w](\+l)  (j-l)= 
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Design  charts  developed  for  these  solutions  are  shown  in  Figure  101,  for  the 
full  afterbody  and  the  no  afterbody  conditions. 

Since  the  assumption  of  conical  flow  assumes  that  no  Mach  lines  emanating 
from  the  fin  tip  leading  edge  projects  into  the  body  region  of  fin  interference, 
this  condition  then  establishes  the  following  criteria  for  selecting  Kg^ 
between  slender-body  theory  and  conical  flow  theory  for  triangular  and  non- 
triangular  planforms  in  supersonic  flow: 
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The  body-fin  carryover  interference  for  the  case  of  a  low  aspect  ratio 
delta  wing  mounted  on  a  conical  body  has  been  determined  by  Spreiter, 
Reference  169.  This  slender-body  theory  solution  correlates  well  with  experi 
mental  data  throughout  the  speed  range.  Figure  102  is  reproduced  from  Datcom 
and  presents  the  interference  factor  such  that 

K(WB)3  ^aW^cAl 

The  work  of  PNK  was  extended  by  Vukelich  and  Williams,  Reference  170, 
to  avoid  the  inaccuracy  of  interpolating  between  full  and  no  afterbody 
interference  values.  This  was  done  by  changing  the  conical  flow  theory 
integral  limits  to  match  the  physical  situation.  The  integral  limits  for 
general  afterbodies  are  described  in  Figure  103. 


Body/Fin  Interference  (with  deflection)  -  PNK  used  slender-body  theory 
to  determine  body  interference  on  a  fin  for  variable  fin  incidence  at  a=0 
to  be 
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The  effect  of  fin  incidence  on  body  carryover  lift  is  derived  from  the 
reciprocal  theorem  of  Reference  150  which  states  that  for  cylindrical  body/ 
fin  combinations  the  following  equality  is  valid  under  the  assumptions  of 
slender-body  theory, 

kaon-i —  ktr,at 


The  slender-body  values  of  incidence  interference  are  presented  in  Figure  104 
Corresponding  linear  theory  values  are  shown  in  Figure  105  for  rectangular 
fin-body  combinations.  Slender-body  values  are  used  for  all  fins  except 
rectangular  planforms  where  &A>2  in  supersonic  flow. 

Empirical  values  of  body-fin  interference  due  to  deflection  (k^)  are 
available  in  tabular  form  (Reference  141)  for  the  transonic  Mach  numbers 
(0.8-1. 3)  and  is  assumed  constant  for  all  deflections  and  fin  characteristics 
due  to  lack  of  available  data.  Between  Mach  1.3  and  3.36  the  lack  of  data 
has  led  Nielsen  to  suggest  the  following  relationships  due  to  the  similarity 
rule  (SBT=slender-body  theory): 
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If  *W  =•  <kw>SBT  THEN  kw  -  <kw>SBT 

If  **  <  <kw>SBT  THEN  kw  "  *w 


Above  Mach  3.36,  slender-body  (K^),.g-|-  values  are  always  used. 

Center  of  Pressure  -  The  PNK  theory  assumes  the  fin  center  of  pressure 
is  unaffected  by  carryover  except  for  two  special  cases  for  which  solutions 
have  been  obtained.  One  case  is  for  rectangular  fins  for  which  Figure  106  pre¬ 
sents  the  center  of  pressure  variation  at  supersonic  speeds  as  derived  by 
linear  theory.  Reference  167.  This  method  is  valid  for  4  >  BA  >2. 
outside  this  range  are  considered  equal  to  fin  alone  variations.  Figure  107 
presents  the  case  of  triangular  planform  fins  with  no  trailing  edge  sweep 
as  derived  from  slender-body  theory. 

The  supersonic  center  of  pressure  on  the  body  due  to  the  presence  of  a 
fin  was  derived  by  PNK  from  conical  flow  theory.  For  supersonic  fin  leading 
edges,  the  result  is 
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c  leading  edges  the  relationship  is. 
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Hence,  the  center  of  pressure  for  fin  interference  on  the  body  is  defined  by 


( mg+l )— 


t  ^ri 

Kn.w)LwCt 


The  PMK  derived  body  interference  center  of  pressure  between  full  and  no 
afterbody  cases  were  improved  for  finite  afterbodies  by  the  same  approach 
illustrated  in  Figure  103. 

For  subsonic  body  center  of  pressure  due  to  the  presence  of  a  fin,  Pf!K 
offers  an  approximate  method  based  on  lifting-line  theory.  This  method  pro¬ 
jects  an  image  quarter-chord  line  with  elliptical  loading  onto  the  body  in 
the  cross-flow  plane.  The  resulting  equation  is. 
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and  is  valid  for  eA>4.0  and  r/s<0.5.  Extrapolation  up  to  r/s=0.8  gives  good 
results.  When  bA<4.0,  the  interference  center  of  pressure  is  obtained  by 
interpolating  from  the  slender-body  value  at  bA=0.  The  theoretical  3A=0 
aerodynamic-center  locations  are  shown  in  Figure  108  from  Datcom.  Other 
values  are  obtained  from  Figure  109. 

The  work  of  PNK  is  the  basis  for  methods  to  determine  body-fin  inter¬ 
ference  for  missile  configurations  today.  An  improvement  to  the  method  was 
made  by  Moore  in  Reference  172  when  he  extended  the  approach  to  handle  fins 
with  swept  trailing  edges.  Figure  110  illustrates  the  procedure  for  deter¬ 
mining  interference  lift  using  slender-body  theory  for  swept  trailing  edge 
fins.  The  approach  defines  a  pseudo-panel  which  can  be  analyzed  using  slender- 
body  theory.  It  is  assummed  that  the  interference  lift  is  concentrated  at 
the  wing  root.  Hence,  the  interference  lift  is  directly  proportional  to  the 
actual  chord-to-pseudo-chord  fraction.  Although  this  seems  like  a  rather 
sever  assumption,  there  is  good  agreement  with  experimental  data.  The  technique 
is  described  by  the  following  equations: 

G=  Cril/Cn 

1  ttm  (  [*.<•, 1  l~  1  /)G 

Although  the  PNK  results  are  generally  limited  to  angles  of  attack  less 
than  ten  degrees,  and  deflection  angles  less  than  15  degrees,  it  is  recommended 
that  they  be  used  in  lieu  of  empirical  results.  Adequate  theoretical  methods 
do  not  exist  at  higher  angle  of  attack  or  deflection  angles. 

Effect  of  Angle  of  Attack  -  Empirical  techniques  for  calculating  fin-body 
interference  have  been  developed.  These  approaches  reduce  appropriate  incre¬ 
ments  for  normal  force  and  moment  interference  from  test  data  of  body  and  fin 
alone  at  angle  of  attack.  Examples  of  this  type  approach  are  that  of  Baker 
(AEDC),  Reference  74,  Aiello  (Martin  Marietta),  Reference  58,  and  Nielsen, 
Reference  141.  These  approaches  are  limited  by  the  bounds  of  the  data 
available  which  are  summarized  in  Figure  111.  It  is  important  to  note  that 
these  data  do  cover  angles  of  attack  well  outside  the  linear  range  assumed 
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for  slender-body  and  linear  theory  and  thus  provide  the  non-linear  effect 
of  angle  of  attack.  Typical  angle  of  attack  effects  on  body-fin  interference 
are  shown  in  Figure  112.  The  trends  are  very  nonlinear  at  intermediate  angles 
of  attack,  even  showing  adverse  interference  on  the  wing  at  certain  conditions. 
These  effects  must  be  handled  empirically  since  no  theoretical  methods  are 
available.  The  empirical  data  base  should  be  expanded  further  to  include  a 
greater  range  of  applicability. 

Panel-Panel  Interference  -  The  theoretical  methods  described  to  this 
point  assume  planar  fin  orientation,  and  only  reflect  small  angle  of  attack 
effects.  The  next  step  in  the  interference  analysis  determined  the 
impact  of  cruciform  fins  and  out  of  plane  effects  (e.^O).  An  approach  to 
the  analysis  of  arbitrarily  deflected  cruciform  fins  at  angle  of  attack  and 
roll  has  been  developed  by  Nielsen,  Reference  141  and  is  termed  the  "equivalent 
angle  of  attack"  (aeq)  concept.  The  term  aeq  is  defined  as  that  angle  of  attack 
of  the  fin  alone  for  which  its  normal  force  is  that  of  the  "influenced"  fin 
accounting  for  the  various  interference  effects.  The  differences  between 
planar  and  cruciform  configurations  is  panel-panel  interference  which  results 
from  a  couplinq  of  the  sidewash  velocities  due  to  angle  of  attack  and  sideslip. 

The  panel-panel  interference  due  to  fin  deflection  of  other  cruciform 
fins  has  been  calculated  using  slender-body  theory  as  an  incremental  a  , 
and  is  presented  in  Figure  113.  All  fins  are  assumed  to  be  within  each 
others  region  of  influence  at  subsonic  and  transonic  Mach  numbers.  For 
supersonic  flow,  a  technique  is  provided  that  determines  the  ratio  of  area 
of  influence  of  a  fin  to  the  total  fin  area.  The  application  of  the  equiv¬ 
alent  angle  of  attack  approach  to  data  within  the  data  base  has  shown  reason¬ 
able  ability  to  duplicate  fin  nonlinearities  at  angle  of  attack  and  arbitrary 
roll  as  shown  in  Figure  114. 

An  approach  similar  in  concept  to  aeq  is  that  of  Oberkampf,  Reference  159. 
Figure  115  defines  the  relevant  geometry.  An  effective  leading  edge  sweep  angle 
is  determined  for  each  fin  by  the  relationship, 

«  *12  2  \ 

A^  *  cos*  (cos  4  cos  -  sin  sin  *  sin  A^  ♦  cos  sin  $  cos  A^J 

This  expression,  and  one  for  the  trailing  edge  effective  sweep  angle,  result 

is  an  effective  aspect  ratio  defined  by, 
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♦  tan  Afcc  -  tan 


Oberkampf  applied  the  effective  aspect  ratio  to  a  subsonic  lift  technique 
based  on  the  Polhamus  Suction  Analogy  with  mixed  results.  Figure  116  shows 
that  the  aeq  approach  of  Nielsen  predicted  better  roll  characteristics  than 
the  equivalent  aspect  ratio  concept  of  Oberkampf  in  the  roll  angle  regime  of 
most  interest.  The  aeq  approach  is  preferred. 

Another  multi-fin  interference  approach  has  been  proposed  by  Darling, 
Reference  89.  He  suggests  the  lift  increase  due  to  the  addition  of  fins  to 


a  cruciform  arrangement  should  be, 

Subsonic 

Supersonic 

^NcAfins^Wafins 

1 .37 

1 .50 

^CNct^8FINS^CNa^4FINS 

1 .62 

2.00 

However,  six  or  eight  fin  panels  in 

combination  are  rare. 

Since  the  effect 

of  body  radius  to  span  ratio  can  be  substantial,  it  is  recommended  that  these 
results  be  used  as  a  first  approximation. 

Fin  Gaps  -  Gaps  between  the  fin  root  and  the  body  mold  line  cannot  be 
avoided  with  all  movable  control  surfaces.  Limited  results  from  tests  has 
shown  that  the  decrease  in  body/fin  carryover  is  relatively  small  for  typical 
gap  widths.  No  method  is  recommended  to  account  for  fin/body  gaps,  but  re¬ 
ference  to  typical  experimental  results,  as  a  design  guide,  are  recommended 
for  inclusion  in  Missile  Datcom. 

Body  El  1 i pti city  -  Although  the  PNK  results  have  been  formulated  for 
circular  bodies,  they  are  often  applied  to  other  body  shapes  as  well.  The 
body  width  at  the  fin  panels  is  chosen  as  the  body  diameter  for  the  carryover 
calculations.  Krieger,  Reference  49,  has  used  the  method  of  Jorgensen  to 
determine  the  effect  of  elliptical  bodies;  the  ratio  ( cN^b ^ ELLI P^ ^ CnWB ^CI R 
was  computed  and  has  been  shown  to  correlate  well  with  test  results.  It  is 
recommended  that  this  approach  be  used  in  Missile  Datcom. 
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Interference  Drag  -  A  limited  amount  of  methodology  is  available  for 
boattail-base-fin  interference.  Darling  provides  an  approach  to  determine 
the  effect  of  a  boattail  on  subcaliber  fin  lift  at  subsonic  speeds;  a  sub¬ 
caliber  fin  is  one  whose  span  is  less  than  the  maximum  body  diameter,  such 
as  a  small  fin  on  a  boattail.  The  fin  effectiveness  is  reduced  by  a  factor 
which  depends  on  the  fin  nose  bluntness.  The  fin  lift  is  calculated  using 
the  boattail  diameter  for  carryover  purposes,  then  is  modified  empirically 
by  the  factor  as  shown  in  Figure  117. 

Two  empirical  techniques  are  available  that  modify  base  drag  due  to 
the  presence  of  tail  fins.  Reference  14  suggests  the  change  in  base  drag 
is  a  function  of  fin  thickness-to-chord  ratio  and  Mach  number  such  that 

ACdBASE  =  (0.825/M2-0.05/M)n 

where  n  =  number  of  fins 

in  the  range  1 .4  <  M  <  2.8 

Moore,  Reference  172,  gives  the  change  in  base  drag  to  be 

(A C^)/=  -  f  (A  C„/V/c) )  <r/c) -O.lx/c]-,  l/c>O.Jx/c 

(A t/c<OJ  x/c 


This  approach  is  empirical  and  assumes  the  fins  are  flush  with  the  base. 

Figure  118  shows  the  Mach  number  trend  and  the  means  to  extrapolate  the 
results  for  non-flush  mounted  fins.  This  method  is  recommended. 

Sunmary  -  The  primary  areas  of  concern  for  methods  development  or  improve¬ 
ment  are  in  the  area  of  high  angle  of  attack  carryover  and  panel-panel  inter¬ 
ference  for  generalized  body-fin  combinations.  The  theoretical  methods 
presented  here  are  generally  applicable  to  typical  missile  configurations 
at  angles  of  attack  of  10°  or  less  and  fin  deflection  angles  of  15°  or  less. 

The  results  of  Hill  and  Kaattari ,  Reference  173,  are  presented  in  Figure  119, 
and  shows  the  non-linear  behavior  at  high  angles  of  attack  and  incidence. 

Enough  data  is  available  to  develop  empirical  nonlinear  methods  for  typical 
low  aspect  ratio  fin-body  combinations.  Additional  data  should  be  sought  to 
extend  the  applicability  of  such  a  method  even  further.  Being  able  to  predict 
this  nonlinear  behavior  is  crucial  to  configuration  stability  and  control 
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analysis  at  high  angles  of  attack.  Most  configurations  utilize  the  fully- 
movable  panel  for  control  and  therefore  require  adequate  definition  of  carry¬ 
over  nonlinearities.  In  addition,  only  planar  or  curciform  panels  (panels 
90  degrees  apart)  are  treated.  There  are  few  methods  available  which  describe 
the  influence  of  other  panel  arrangements,  such  as  tri-form  or  arbitrary  di¬ 
hedral  angle  cruciform  panels.  One  may  use  the  inertial  coefficients  derived 
by  Nielsen  in  the  text  "Missile  Aerodynamics",  Reference  12,  and  determine 
the  approximate  effect  from  slender  body  theory.  This  has  been  done  and  the 
results  are  presented  in  Figure  120  in  a  ratioed  form  to  that  for  a  cruciform 
configuration  in  the  "plus"  orientation. 

It  is  recommended  that  a  mix  of  PNK  theoretical  and  other  empirically 
derived  results  for  carryover  be  incorporated  in  Missile  Datcom.  The 
deficiency  in  methods  due  to  high  angle  of  attack  or  panel  incidence  should 
be  corrected  for  maximum  accuracy  and  utility  of  Missile  Datcom. 

Methods  could  not  be  found  which  address  carryover  interference  of 
non-straight  tapered  panels  or  those  panels  which  are  swept-forward.  These 
designs  are  not  at  all  uncommon  in  missile  design  today.  It  is  recommended 
that  these  deficiencies  be  corrected. 

5.3  VORTEX  INTERFERENCE 

In  analyzing  configurations,  the  effect  of  vortices  is  an  important 
aerodynamic  consideration  that  cannot  be  neglected.  Four  specific  tasks 
must  be  performed  to  analyze  the  vortex  effects  as  follows: 

a)  Position  of  vortex  shedding,  either  position  on  the  body  or  span 
location  on  the  wing 

b)  Wing  produced  flowfields;  the  wing  flowfield  can  be  modeled  as  a 
flat  sheet  downwash  field,  as  a  fully-rolled-up  line  vortex,  or  in 
combination. 

c)  Vortex  tracking;  mapping  the  vortex  position  through  the  flowfield 
and  determining  its  proximity  to  the  configuration 

d)  Vortex  strength;  determination  of  the  vortex  circulation  strength 

Figure  121  illustrates  the  number  and  type  of  vortex  interactions  which 

occur  for  a  typical  missile  design.  Vortices  generated  by  the  body  influence 
any  fins  present  (the  effect  on  the  body  aft  of  the  separation  point  is 
accounted  for  through  use  of  viscous  cross  flow),  whereas  fin  produced 
vortices  influence  both  the  body  and  any  aft  lifting  surfaces.  The  system 
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of  vortices  can  get  extremely  complex  through  addition  of  inlets,  launch  lugs, 
conduits  or  even  body  shape.  Axisymnetric  bodies  without  protuberances  have 
been  explored  in  sufficient  detail  to  develop  methodology  to  track  and  deter¬ 
mine  the  strength  of  the  body  nose  vortices.  The  method  in  Datcom  Section 
4. 3. 1.3  and  the  empirical  correlations  of  Nielsen,  briefly  presented  in 
Figure  122,  illustrate  the  methods  in  use.  No  design  methods  are  available 
for  general  shaped  configurations. 

The  vortex  effect  of  forward  lifting  surfaces  has  also  been  explored  in 
detail.  The  first  such  efforts  were  presented  by  Sprieter  and  Sacks, 

Reference  175,  in  1951  and  Decker,  Reference  174,  in  1956.  These  methods 
considered  the  downwash  behind  wings  and  characterized  the  vortex  effect 
as  either  a  sheet  or  a  fully-rolled  up  vortex  core.  The  presence  of  a  vor¬ 
tex  sheet  is  dependent  upon  panel  aspect  ratio;  the  roll-up  into  a  vortex 
is  inversely  proportional  to  aspect  ratio  and  directly  proportional  to  lift. 

For  larger  aspect  ratio  panels,  or  low  lift  coefficients,  the  wing  will  shed 
a  trailing  vortex  sheet  which  slowly  rolls  into  a  fully  developed  vortex  core 
at  angle  of  attack.  Since  missile  panels  are  frequently  low  aspect  ratio,  it 
is  appropriate  to  consider  only  the  fully-rolled  vortex  concept,  which  was 
explored  in  detail  by  Pitts,  Nielsen  and  Kaattari  in  NACA  1307  (Reference  167). 
However,  for  cruise  missiles,  or  other  aircraft-type  designs,  the  "sheet"  con¬ 
cept  is  appropriate  and  the  Datcom  method  should  be  retained  for  that  purpose. 

The  method  of  NACA  1307  is  highly  appropriate  for  classical  missile 
configurations.  Single  vortex  cores  are  shed  from  each  panel  and  trail  aft 
approximately  along  the  free-stream  velocity  vector  and,  hence,  are  dependent 
upon  the  total  angle  of  attack.  Since  the  method  of  NACA  1307  assumes  fin 
panels  on  both  sides  of  the  body,  the  vortex  effect  on  one  panel  at  dihedral 
must  be  handled  with  care.  A  method  available  to  perform  this  task  is  pre¬ 
sented  in  Datcom  and  shown  in  Figure  123.  This  process  assumes  that  super¬ 
position  of  multiple  vortex  effects  can  be  performed.  In  reality,  the  vortices 
interact  with  each  other  and  do  not  necessarily  follow  the  free-stream 
velocity  vector  at  angle  of  attack.  Simulations  of  these  vortex  tracks 
should  provide  more  accurate  results.  It  has  not  been  determined  what 
degree  of  accuracy  improvement  is  realized  compared  to  computing  costs. 

Such  a  determination  is  beyond  the  scope  of  this  study,  but  should  be 
quantified  during  the  development  efforts.  Since  the  vortices  tend  to 
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follow  the  velocity  vector,  track  nearly  parallel  to  the  chord  of  a  lifting 
surface  in  close  proximity,  and  continue  aft,  it  may  be  possible  to  devise  a 
simple  but  suitably  accurate  method  using  such  observations,  particularly 
for  those  configurations  involving  tandem  lifting  surfaces  where  otherwise 
the  mathematics  can  be  extremely  complex. 

As  illustrated  in  Figure  124,  a  body  at  angle  of  attack  sheds  a  number 
of  vortices  which  have  been  observed  to  separate  at  predictable  positions 
along  the  body.  These  body  vortices  can  be  significant  for  tail  orientations 
other  than  cruciform  "plus"  or  at  mixed  angle  of  attack  and  sideslip  angle. 
Figure  125  illustrates  the  panel  local  angle  of  attack  for  both  the  vortex 
free  and  vortex  present  cases.  For  leeside  flow  conditions  the  total  vortex 
effects  results  in  as  much  as  a  12  deg.  change  in  local  angle  of  attack.  The 
magnitude  of  error  induced  in  neglecting  the  afterbody  vortices  is  unknown, 
though  expected  to  be  significant  at  the  higher  angles  of  attack.  A  more 
significant  question  in  this  highly  empirical  regime  is  the  justification 
of  its  incorporation  for  design.  Can  a  suitable  design  be  evolved  without 
its  use?  This  question  must  be  addressed  in  the  Missile  Datcom  development. 
It  is  recommended  that  changes  in  local  dynamic  pressure  at  the  lifting  sur¬ 
faces  due  to  flow  conditions  (compressive  or  expansive  flow)  given  in  Datcom 
be  retained  and  expanded  for  missile  design  purposes. 

It  is  recommended  that  the  wing  vortex  tracking  and  vortex  strength 
methods  of  NACA  1307  and  Datcom  be  used.  It  is  also  recommended  that  the 
empirical  correlations  of  nose  vortex  strength  and  shedding  position  of 
Nielsen  be  adopted.  Until  the  complex  vortex  tracking  methods  can  be 
evaluated  with  respect  to  accuracy  and  cost,  tracking  along  the  velocity 
vector  is  recommended.  Those  recommended  techniques  are  given  in  Table  16. 
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Figure  98.  NACA  1307  Interference  Regimes 
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Figure  99.  lift  ratios  Kw<„  and  k*w>  — slender-body  theory- 
fixed  INCIDENCE— ALL  SPEEDS 
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(b)  PLANAR  MODEL 

Figure  TOO.  Carry-Over  Interference  Model 
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a)  ^b(W)  Design  Chart  for  Full  Afterbody 
(Conical  Flow  Theory) 


b)  Kg(w)  Design  Chart  for  No  Afterbody 
(Conical  Flow  Theory) 

Figure  101.  Body  in  Presence  of  Wing  Carry-Over-Supersonic 
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Figure  102.  Spreiter  cone-delta  wing  interference  factor 
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Figure  103.  Integration  Limits  for  General  Afterbody  Geometries 
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Figure  104.  urr  ratios  kw«  and  k*wi — slender-body  theory 
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Figure  1 05 .  lift  on  body  in  presence  of  wing-variable  incidence 

SUPERSONIC  SPEEDS 
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U/O  *(«)•  °*  WO  n$ it 


Figure  106. Rectangular  Fin  Center  of  Pressure  Variation  due  to  Incidence  and  Body 

Carryover  at  Supersonic  Speeds 
(Linear  Theory) 


Figure  107.  Triangular  Fin  Center  of  Pressure  Variation  due  to  Incidence  and 

Body  Carryover 
(Si ender-Body  Theory) 
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Figure  108.  THEORETICAL  AERODYNAMIC-CENTER  LOCATIONS  FOR  0A*  *  0 
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(b)  NO  AFTERBODY 


Figure  109. 


AERODYNAMIC-CENTER  LOCATIONS  FOR  LIFT  CARRYOVER  OF  WING  ONTO 
BODY  AT  SUPERSONIC  SPEEDS  WHEN  0A«  (1  +  X#)  (\  +  JcOTaT)  >  40 
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Procedure  used  to  calculate  interference  lift  for  wings  with 
sweptback  trailing  edges  when  slender  bod>  theory  is  used:  a)  wing  for 
which  interference  lift  is  desired;  b)  assumed  slender  body  represen¬ 
tation. 


Figure  110.  Swept  Trailing  Edge  Interference  Approach  By  Moore 
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Figure  111.  Empirical  Interference  Methods  Data  Base  Limitations 
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Wing-body  interference  factor 


Figure  112.  Effect  of  Angle  of  Attack  on  Body-Fin  Interference 
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Figure  113.  Panel-panel  Interference  due  to  deflection  -  Cruciform  Fins 

(SI  ender-Body  Theory) 
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Figure  114.  Equivalent  Angle  of  Attack  Concept  Comparison 


Georetry  for  Determining  Effective 
Aspect  Ratio 

Figure  115.  Oberkampf  Geometry 
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Figure  116.  Comparisons  of  Nielsen  vs  Oberkampf  Approaches 
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Figure  117.  SUBCAUBER  FIN  EFFECTIVENESS  FACTOR  VARIATION 
WITH  NOSE  BLUNTNESS 


204 


Base  pressure  coefficient  change  with  fins  located  flush  with  tee. 


j£  *  distance  from  base  to  tail  trailing  edge  where  a  fin  of  given 
thickness  has  no  affect  on  base  pressure. 


Figure  118. Effect  of  Fins  on  Base  Drag 
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Figure  120.  Apparent  Hass  Effect  of  Number  of  Fin  Panels  Versus 
Diameter  to  Span  Ratio 
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NOSE  VORTEX 
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Figure  121.  Vortices  Present  for  a  Typical  Configuration 
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Figure  122. 


Body  Nose  Vortex  Correlations 
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TYPICAL  ORIENTATION  OF  BODY  VORTICES  WITH  RESPECT 
TO  TAIL  PUNE  FOR  CON  KK  IF  RATION  AT  COMBINED  ANGLES 
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Figure  123.  Datcom  Vortex  Effect  Due  to  Panel  Orientation 
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Figure  124.  Afterbody  Vortices  Can  be  Well  Predicted 


212 


Figure  125.  Effect  of  External  Vortices  on  Panel  Angle  of  Attack 
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SECTION  6 

CONFIGURATION  SYNTHESIS 


It  is  recommended  that  Missile  Datcom  utilize  the  component  buildup 
approach  for  missile  configuration  synthesis.  The  equivalent  angle  of  attack 
approach  with  a  total  angle  of  attack,  bank  angle  coordinate  system  is  recom¬ 
mended  to  facilitate  interference  calculations  at  arbitrary  roll  angles.  The 
computer  program  should  be  structured  so  that  user  supplied  experimental  re¬ 
sults  available  for  a  particular  configuration  component  can  be  used.  This 
then  establishes  the  need  for  a  modular  program  in  which  the  aerodynamic  com¬ 
putations  are  isolated  in  individual  routines.  The  program  structure  would 
be  similar  to  the  Missile  Datcom  handbook  outline  presented  in  Section  2. 

Configuration  synthesis  is  the  process  of  estimating  the  aerodynamic 
characteristics  of  a  complete  configuration.  Typically,  synthesis  involves 
the  summation  of  the  aerodynamic  characteristics  of  each,  of  the  configuration 
components  and  then  adding  mutual  interference  effects.  Two  approaches  exist 
for  combining  components  to  obtain  total  configuration  aerodynamics.  Figure 
126  classifies  the  two  methods  for  normal  force  of  a  wing-body- tai 1  configura¬ 
tion.  The  "classic"  approach  is  the  linear  lift  method  presented  in  NACA 
1307  extended  to  non-linear  angles  of  attack.  The  carryover  effects  between 


wing/body  and  tail/body  are  represented 


Kg(T))  to  the  panel  alone  characteristics.  Since  the  body  generates  a  pair 
of  trailing  nose  vortices,  these  vortex  effects,  Ivw(b)  anc*  *VT(b)’  on 
panels  are  separately  computed  and  summed  with  the  wing  vortex  effect,  I Vj( vj)  5 
on  the  tail  surfaces.  Although  wing  and  tail  vortices  are  present,  their 
effect  on  the  body  is  assummed  negligible.  The  body  nose  vortex  effect  on  the 
body  is  modeled  using  the  viscous  cross-flow  procedure  described  by  Allen 
and  Perkins,  and  presented  in  Section  3. 

The  "equivalent  angle  of  attack"  concept  by  Nielsen  (Reference  176) 
identifies  the  same  interference  contributors,  but  combines  the  effects  of 
carryover  and  vortices  on  the  panel  local  angle  of  attack.  The  sources 
of  interference  considered  are  1)  panel-panel  interference,  where  the 
change  in  loading  of  one  panel  affects  another,  2)  nose  vortex-fin  interference, 
3)  wing  vortex-tail  interference,  and  4)  afterbody  vortex-tail  interference. 

All  of  these  effects  are  evaluated  as  an  incremental  effect  on  equivalent 
angle  of  attack.  Once  the  equivalent  angle  of  attack  has  been  evaluated, 
the  panel  normal  force  is  obtained  from  the  empirical  or  theoretical  panel 
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alone  characteristics.  This  technique  is  advantageous  since  at  transonic 
speeds  only  empirical  panel  aerodynamic  characteristics  are  available. 

These  two  methods  are  easier  to  apply  if  aerodynamic  characteristics 
are  estimated  as  a  function  of  total  angle  of  attack  and  roll  angle.  Bank- 
to-tum  and  skid-to-turn  missile  configurations  have  different  angle  of 
attack  and  sideslip  angle  envelopes,  as  shown  in  Figure  127.  The  synthesis 
technique  chosen  must  encompass  the  most  extreme  condition,  such  as  20  degrees 
of  angle  of  attack  and  sideslip.  For  those  cases,  the  "equivalent  angle  of 
attack"  approach  is  the  easiest  to  apply.  The  loads  on  the  individual 
panels  are  directly  determined;  therefore,  the  panel  hinge  moment  and  bending 
moment  can  be  computed.  This  system  is  best  for  the  vortex  calculations, 
since  the  vortex  paths  are  related  to  the  configuration  total  angle  of  attack. 
References  177-185  describe  in  detail  the  synthesis  of  configurations  at 
angle  of  attack  and  bank  angle  and  serve  as  excellent  documentation  for  the 
effects  of  high  body  incidence.  For  arbitrary  shaped  bodies,  an  angle  of 
attack/sideslip  angle  system  is  still  required  because  of  the  non-axial 
symmetry  of  the  body.  It  is  difficult  to  obtain  total  angle  of  attack 
characteristics  of  arbitrary-shaped  bodies  throughout  the  required  envelope. 
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CLASSIC  APPROACH: 

CN  =  CMBODY  +  [Kw(B)  +  Kb(w)]  CnWING  + 

[kT(B)  +  KB(T)]  cntail  + 

’VW(B)  +  IvT(B)  +  IvT(W) 


NIELSEN  EQUIVALENT  ANGLE  OF  ATTACK  APPROACH: 

CN  "  CNB0DY  +  KB(W)  CNwing  +  KB(T)  CNtail 

+ycN  (0£cn)u  +y^cN  («eq)t 
nwing  eq  w  ntail  t 

- - V - - 

INCLUDES  EFFECT  OF  CL, 5,<j)  AND  ALL 
EXTERNAL  VORTICES  (NOSE,  AFTERBODY,  WING) 


Figure  126.  Configuration  Synthesis  Concepts 
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Figure  127.  a  and  g  Envelopes 
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SECTION  7 

PROPULSION  SYSTEM  EFFECTS 


Missiles  generally  use  one  of  three  propulsion  system  types:  (a)  rocket, 

(b)  turbojet,  or  (c)  ramjet.  Propulsion  system  effects  on  aerodynamics  is 
perhaps  the  least  developed  area  in  missile  aerodynamics,  due  to  the  complexity 
and  variety  of  propulsion  systems.  Many  missile  aerodynamic  prediction 
codes  ignore  the  propulsion  effect  on  aerodynamics,  although  it  is  often 
significant.  Section  3  has  described  the  effect  of  jet  plume/exhaust  on 
boattail  wave  drag.  This  section  will  describe  those  other  effects  which  will 
be  required  for  Missile  Datcom. 

PI ume/Ai rframe  C^,  Cm  -  A  major  problem  for  the  aerodynamicist  are  the 
plume/airframe  interaction  effects.  As  the  vehicle  climbs  in  altitude,  the 
exhaust  plume  from  a  rocket  will  expand  producing  the  effect  illustrated  in 
Figure  128.  Normal  force  and  pitching  moment  can  be  significantly  altered. 

The  classic  means  of  solving  this  problem  is  to  assume  that  the  plume  effect 
is  similar  to  that  of  a  transverse-jet  control  device,  but  requires  detailed 
computations.  The  recommended  approach  is  that  presented  by  Aiello  and 
Bateman,  Reference  58,  where  the  incremental  body  normal  force,  body  center 
of  pressure,  and  tail  normal  force,  have  been  empirically  correlated  over 
the  Mach  range  from  0.6  to  2.2  through  180  degrees  angle  of  attack.  The 
extrapolation  of  these  results  to  higher  and  lower  Mach  numbers  should  be 
investigated.  It  will  be  necessary  to  obtain  additional  data  for  further 
method  development. 

For  an  airbreathing  propulsion  system,  the  effect  of  jet  plume/airframe 
interaction  should  be  minimal.  The  captured  air  will  be  exhausted  at  a 
pressure  close  to  that  for  the  free-stream,  and  pluming  will  be  small. 

Operational  airbreathing  propulsion  systems  typically  have  high  pressure 
recovery  coefficients.  For  these  propulsion  systems  the  effect  of  plume/ 
airframe  interaction  can  be  assumed  to  be  negligible  for  preliminary  design 
purposes. 

Plume/Base  Interaction  -  The  base  pressure  methods  described  in  Section 
3  were  for  jet-off  characteristics.  The  jet  effects  on  base  drag  occur  for  all 
propulsion  systems  and  were  theoretically  analyized  by  Addy,  Reference  187, 
and  Korst,  Reference  188.  These  methods  assume  a  strong  plume/freestream  inter¬ 
action  such  as  that  for  an  over-expanded  plume  typical  of  rockets.  However,  they 
do  not  address  the  base  aspiration  effect  typical  to  airbreathing  missiles  where 
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the  plume  may  not  interact  strongly  with  the  free-stream.  Although  therKorst 
technique  is  automated,  Reference  189,  it  is  a  complex  program;  a  method  of 
characteristics  computation  is  performed  to  determine  the  plume  shape.  Various 
simplified  forms  of  the  technique  are  available,  one  such  method  being  that  of 
Reference  190.  This  technique  is  recommended  for  rocket  propulsion  systems.  The 
plume  shape  and  pressure  have  been  correlated  into  simple  relationships,  and  the 
calculation  of  base  pressure  becomes  a  routine  computation. 

For  airbreathing  inlets,  a  simple  approach  is  recommended.  For  turbo¬ 
jets  and  ramjets  the  power-on  base  drag  is  usually  higher  than  that  for  power- 
off,  but  the  annular  base  area  is  usually  smaller  than  the  power-off  area. 
Therefore,  the  inaccuracy  is  using  power-off  base  drag  from  Stoney,  Reference 
35,  is  usually  small.  If  necessary  a  compendium  of  data  and  empirical  tech¬ 
niques  are  available  in  Reference  32  for  more  detailed  estimates. 

Inlet  Effect  -  Airbreathing  inlets  change  the  body  shape  to  non- 
axi symmetric  and  are  usually  treated  empirically.  A  limited  number  of 
methods  are  available  in  Reference  192  using  empirical  correlations. 

Generally,  detailed  analysis  of  the  inlet  components  (e.g.,  cowl,  duct, 
boundary  layer  diverter)  involves  a  number  of  theoretical  and  empirical 
methods  (which  have  been  summarized  in  Reference  191  for  several  inlet 
types)  similar  to  the  component  build-up  approach  for  bodies.  Use  of 
sophisticated  panel  techniques  have  also  been  employed.  The  simple 
inlet/circular  body  build-up  method  is  seen  as  an  ideal  candidate  for 
preliminary  design.  A  compendium  of  data  for  2-D,  axisymmetric  and  chin 
inlet  designs  is  available  to  develop  this  technique.  The  chin  inlet 
methods  of  Reference  58  are  recommended  for  this  task.  A  significant 
amount  of  test  data  has  been  obtained  through  a  cooperative  effort  between 
NASA-Langley  and  the  Naval  Weapons  Center  in  China  Lake,  CA.  These  data 
contain  configuration  build-up  data  for  a  missile  configuration  with  the 
2-D  and  axisymmetric  inlet  types.  The  inlets  were  mounted  on  the  body 
in  various  longitudinal  and  circumferential  positions.  However,  this 
data  has  yet  to  be  put  into  a  form  suitable  for  engineering  design  purposes. 

It  is  recommended  that  an  R&D  task  be  funded  specifically  for  the  purpose 
of  deriving  preliminary  and  conceptual  design  methodology  using  this  data 
base  and  the  methods  of  References  58  and  191.  In  the  interim,  methods  which 
have  been  derived  from  empirical  results  and  presented  in  the  "Ramjet  Design 
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Handbook",  Reference  192  should  serve  as  excellent  design  methodology  when 
data  is  unavailable. 

It  is  recommended  that  the  component  build-up  approach  be  applied  to 
inlet  effects  on  configuration  aerodynamics.  It  is  also  recommended  that 
the  automated  version  of  Missile  Datcom  employ  a  structure  which  will  allow 
substitution  of  test  data  when  it  is  available. 
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Figure  128.  Plume/Airframe  Interactions  can  Influence  Aerodynamics 


SECTION  8 

CONTROL  DEVICE  METHODOLOGY 


The  method  recommendations  for  control  device  methodology  are  as  follows: 

(a)  Panel  deflection,  use  of  the  "equivalent  angle  of  attack" 
concept  by  Nielsen 

(b)  Plain  flap;  Datcom 

(c)  Hypersonic  flaps;  Datcom 

(d)  Transverse-jet  control ;  Datcom 

(e)  Strakes;  Aiello  and  Bateman 

These  methods  are  briefly  summarized  in  this  section. 

Only  the  simplest  of  control  devices  are  used  in  missile  design  because 
of  cost  and  complexity  constraints.  The  most  often  used  control  device  is 
the  all-movable  surface  with  deflection  ranges  typically  of  +30  degrees. 

Other  control  devices  less  frequently  used  include  thrust-vector-control 
(TVC),  jet  interaction/reaction  jet  control  (JI/RJC),  plain  trailing  edge 
flaps  and  body  flaps.  More  sophisticated  control  devices  such  as  split 
flaps,  leading-edge  flaps,  or  jet  flaps  are  extremely  rare  and  their  inclusion 
is  not  recommended. 

Panel  Deflection  -  Section  5  on  component  interference  described  in 
detail  the  effect  of  high  angle  of  attack  and  deflection  on  panel  effective¬ 
ness,  emphasizing  the  need  for  carry-over  interference  results  beyond  that 
obtained  through  slender-body  or  conical  flow  theory.  Through  use  of  the 
empirical  results  available  for  a  variety  of  fin  panel  designs,  sufficiently 
accruate  results  in  the  transonic  and  supersonic  Mach  regimes  can  be  obtained. 
Effective  use  of  the  "equivalent  angle  of  attack"  concept  by  Nielsen  allows 
determination  of  panel  hinge  moment  and  bending  moment  characteristics. 
Correlations  with  test  results  have  shown  acceptable  accuracy  for  preliminary 
desi  gn . 

Plain  Flaps  -  The  Datcom  methods  for  plain  flaps  are  recommended.  Typical 
variation  of  the  flap  characteristics  are  presented  in  Figures  129-132.  The 
method  reference  and  limitations  due  to  Mach  number  are  as  follows: 

(a)  Subsonic  (Reference  136) 

(b)  Transonic  (Approximate  technique.  Reference  1) 

(c)  Supersonic  (Reference  193;  supersonic  leading  edge  and  trailing  edge 

and  flap  surfaces;  control  located  on  surface  tip) 
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Typical  design  charts  for  the  Datcom  method  are  presented  in  Figures  133-141. 

The  method  results  in  Figure  142  show  good  correlation  with  test  data. 

Hypersonic  Flap  -  The  Datcom  method  for  a  plain  flap  at  hypersonic  speeds 
(Hypersonic  Flap)  is  recommended.  This  technique  is  for  Mach  numbers  greater 
than  five  and  for  the  most  part  was  taken  from  Reference  194.  A  qualitative 
variation  of  the  pressure  distribution  is  shown  in  Figure  143. 

Transverse-Jet  Control  -  The  transverse-jet  control  methods  of  References 
195-197  have  been  incorporated  into  the  Datcom  method.  These  methods  are 
also  recommended  for  inclusion  in  Missile  Datcom.  The  methods  presented  in 
Datcom  cover  the  Mach  range  from  2  to  20.  A  representative  pressure  distri¬ 
bution  due  to  jet  interaction  with  the  flow  is  shown  in  Figure  144. 

Strakes  -  Strakes  are  sometimes  included  for  the  purpose  of  enhancing 
control  effectiveness.  As  shown  in  Figure  145,  a  strake  produced  vortex  will 
sweep  across  the  major  lifting  surface  and  result  in  increased  lift  effective¬ 
ness  at  much  higher  angles  of  attack.  This  phenomena  benefits  a  bank-to-turn 
missile  configuration  due  to  increased  pitch  normal  force  as  Figure  146 
illustrates.  A  thorough  discussion  of  these  results  is  presented  in  References 
198  and  199.  There  are  no  theoretical  methods  available  to  predict  the 
favorable  interference  effect  of  wing  strakes  but  the  empirical  (subsonic) 
lift  methods  in  Datcom  are  available  for  cranked  or  double-delta  panels 
(straked  wings)  for  angles  of  attack  to  20  degrees.  Isolated  strakes  (i.e. 
not  a  physical  part  of  a  primary  lifting  surface)  are  also  in  coiranon  use  to 
influence  the  stability  characteristics.  These  isolated  strakes  may  be 
treated  through  use  of  very  low  aspect  ratio  wing  theory  and  wing-body  carry¬ 
over  interference.  An  empirically  derived  technique  is  available  from  Aiello 
and  Bateman  (Reference  58)  and  is  shown  in  Figure  147.  Body-strake  aerodynamic 
effects  are  predicted  as  an  increment  to  the  body-alone  aerodynamics.  The 
criteria  which  defines  those  panels  which  are  classified  as  "wings"  or  "fins" 
and  those  which  are  classed  as  "strakes"  is  lacking.  A  cut-off  aspect  ratio 
or  span-to-diameter  ratio  should  be  defined  to  distinguish  between  the  methods 
to  be  utilized.  It  is  recommended  that  the  Aiello  method  be  compared  with 
low  aspect  ratio  wing  theory  and  carry-over  interference  effects  so  that  a 
method  can  be  chosen.  In  addition,  a  criteria  should  be  defined  which 
specifies  those  surfaces  which  should  be  analyzed  using  wing  theory  or  strake 
methodology.  The  Aiello  and  Bateman  method  is  recommended,  subject  to  a  thorough 
quantiative  analysis. 
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Figure  129.  Panel  Section  Lift  Due  to  Plain  Flap,  from  Datcom 


Figure  130.  Panel  Section  Pitching  Moment  Linear  Range,  from  Datcom 
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Figure  131.  Center  of  Pressure  Across  a  Flapped  Surface, 
from  Datcom 


Figure  132.  Typical  Hinge  Moment  Results  for  A  Plain  Flap, 
from  Datcom 
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Figure  133.  Theoretical  Lift  Effectiveness  of  Plain  Trailing-Edge  Flaps 


Figure  134.  Empirical  Correction  for  Lift  Effectiveness  of  Plain 
Trailing-Edge  Flaps 


227 


0  .2  A  .6  .8  1.0 

FLAP-CHORD  RATIO,  cf/c 

Figure  135.  Empirical  Location  of  Center  of  Pressure  of  Incremental  Load 
due  to  Trail ing-Edge,  Mechanical  Flaps 
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Figure  136.  Effect  of  Trailing  Edge  Flap  Deflection  and  Flap-Chord-to- 
Wing-Chord  Ratio  on  Section  Incremental  Pitching  Moment  Due 
to  PI ai n  Flaps 
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Figure  137.  Approximate  Maximum  Control-Surface  and  Deflections  for 
Linear  Control  Characteristics  of  a  Plain,  Sealed  Flap 
(NACA  0009  Airfoil) 
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Figure  138.  Rate  of  Change  of  Section  Hinge  Moment  Coefficient  With 
Angle  of  Attack  for  A  Plain  Flap 
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Figure  139.  Rate  of  Change  of  Hinge-Moment  Coefficient  with  Control 
Deflection  for  a  Plain  Flap 
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Figure  140.  Thickness  Correction  Factor  for  Hinge-Moment  Derivatives 
for  Symmetric,  Circular-Arc  Airfoils 
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Figure  141.  Flap-Chord  Factor 
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Figure  142.  Correlation  of  Datcom  Plain  Flap  Method  for  Lift  With 
Test  Data 
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Figure  144.  Pressure  Distribution  Due  to  Transverse-Jet  Control,  from  Datcom 
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CONSIDERABLE  BENEFIT  TO  BANK-TO-TURN  MISSILES 

-  SWEEP  BREAK  SEVERS  VORTEX  FEEDING  SHEETS 

-  SHED  VORTICES  INCREASE  LIFT  EFFECTIVENESS 


Figure  145.  Wing  Strakes  in  Aerodynamics 
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Figure  146.  Benefit  of  Wing  Strakes 
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Figure  147. 


Aiello  Strake  Method 
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SECTION  9 

DYNAMIC  DERIVATIVE  METHODOLOGY 


The  most  comprehensive  compendium  of  dynamic  derivative  methodology  is 
that  of  Datcom  which  is  valid  only  at  low  angles  of  attack  with  attached  flow 
conditions.  Few  theoretical  methods  are  available  to  compute  the  dynamic 
derivatives  of  configurations  not  covered  by  Datcom.  The  Datcom  methods  are 
sunmarized  in  Table  17.  These  techniques  are  presented  in  the  stability-axes 
coordinate  system.  Transfer  equations  from  stability  to  body  axes  have  been 
recommended  in  Section  2  of  this  report.  The  Magnus  effect  is  discussed  in 
Section  10. 

Body  and  Body  Fin  Configurations  -  Most  of  the  available  methodology 
for  body  configurations  is  derived  from  slender  body  theory  (see  for  example 
References  12,  200,  and  201),  or  correlations  of  large  amounts  of  test  re¬ 
sults".  Although  rigorous  numerical  techniques  are  available,  such  as  in 
References  92  and  201,  their  complexity  is  considered  beyond  the  suitability 
for  Missile  Datcom.  The  recent  works  of  Ericsson  (Reference  204-206)  have 
extended  the  simple  slender  body  theory  and  Newtonian  theory  methods  to  more 
properly  account  for  Mach  number.  The  method  suggested  by  Erricson  for 
bodies  is  (from  Reference  204) 


2 


0.77  +  0.23  M  »  -  -^)2 


based  on 


d 


This  method  has  shown  acceptable  accuracy  to  Mach  6.0  but  further  substantia¬ 
tion  of  the  method  is  required.  Comparison  of  results  with  experimental  data 
and  the  empirical  results  for  the  SPINNER  code  (Reference  208)  is  shown  in 
Figure  148.  In  Reference  206,  Ericsson  presented  methods  which  give  excellent 
results  for  body-wing  pitch  damping  at  all  speeds.  These  methods  are  recommended 
techniques  for  Missile  Datcom.  Other  theoretical  methods  sunmarized  in  Datcom 
are  given  in  References  211-221. 

Other  Configurations  -  The  methodology  available  can  be  summarized  as 
follows : 
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o  DATCOM  -  THEORY  (SLENDER-BODY ,  LIFTING  SURFACE)  (*  SUBSONIC  ONLY) 

WINGS  -  Cmq»  ^Aq*’^Yp>  p »  ^A^* 

BODY  -  Cfjq,  Cm q »  Cn^» 

WING-BODY  -  C^,  Cmq,  C£q,  C[^a,  C^a 

BODY-WING-TAIL  -  CNq,  C^,  CAq,  C*p*,  Cnp*,  CYp*,  Cy/,  C*r*,  CPr*,  CNa,  Cma, 

cAa*.  cY*.  cn*,  cq 

o  "SPINNER"  PROGRAM  -  EMPIRICAL 
EMPIRICAL  DATA  BASE 

Cmq  +  Cmd 

c£,p 

MAGNUS 


A  suitable  means  of  predicting  the  effects  of  forward  lifting  surfaces 
on  the  dynamic  derivatives  of  aft  lifting  surfaces  are  required  for  these 
configurations  which  employ  multiple  lifting  surfaces.  Direct  summation  of 
the  contributions  of  the  configuration  components  provide  an  adequate  first- 
order  approximation,  but  the  interference  effects  between  components  can  be 
substantial.  The  only  comprehensive  interference  methods-  available  are  those 
presented  in  Datcom.  It  is  recommended  that  they  be  employed. 

Angle  of  Attack  Effect  -  No  methods  are  available.  The  primary  dis¬ 
advantage  of  dynamic  derivative  methodology  is  its  restriction  to  low  angle 
of  attack,  attached  flow  conditions.  No  methodology  was  found  that  specifically 
addresses  the  effect  of  body  incidence. 

Jet  Effects  -  No  methods  are  recommended.  Although  jet  pluming  can  affect 
the  dynamic  stability  of  a  missile,  it  cannot  be  evaluated.  It  is  ignored  in 
most  preliminary  or  conceptual  designs. 

Jet  effects  on  dynamic  derivatives  are  rarely  calculated  in  existing 
prediction  codes.  Not  all  derivatives  are  necessary  for  preliminary  or 
conceptual  design,  but  the  major  contributors  to  dynamic  stability:  Cn^,  Cpr, 
C^,  C^,  and  C$,p  should  be  included  in  Missile  Datcom.  An  important  considera¬ 
tion  in  the  choice  of  necessary  derivatives  is  the  missile  configuration  type 
and  guidance  steering  class.  Langham,  Reference  207,  performed  an  analysis  of 
two  missile  configurations:  the  Interlab  Air-to-Air  Technology  (ILAAT)  and 
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the  Aerodynamic  Data  Correlation  (ADC)  configurations,  to  determine  the 
sensitivity  to  the  dynamic  stability  derivatives.  The  following  parameters 
were  found  to  cause  variations  in  missile  motion:  Cm  ,  Cm-,  Cn  ,  Co  ,  Co  and 
Cnp,  the  derivatives  Cmr,  Ci_q,  Cl^,  Cyp»  and  Cyr  were  deemed  insignificant  to 
motion  for  both  bank-to-turn  and  skid-to-turn  configurations.  It  is  recommended 
that  a  complete  compendium  of  dynamic  stability  methodology  be  included.  Al¬ 
though  the  results  of  Langham  will  prove  valuable  for  some  design  configurations, 
a  thorough  analysis  of  a  specific  configuration  is  required  to  quantify  those 
derivatives  which  are  design  critical. 
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TABLE  17  SUMMARY  OF  DATCOM  DYNAMIC  DERIVATIVE  METHODS 
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TABLE  17.  (CONTINUED) 

METHODS  SUMMARY 
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Figure  148.  Ericsson  Method  Comparison 
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SECTION  TO 

SPECIALIZED  METHODOLOGY 


Section  9  of  the  Missile  Datcom  Outline,  Table  3,  covers  "specialized" 
areas  of  missile  aerodynamic  analysis.  Two  specific  'areas  were  identified 
for  inclusion,  (1)  tumbling  motion,  and  (2)  spinning  motion.  Although  not 
required  for  typical  missile  aerodynamic  design,  these  techniques  will  supply 
methodology  to  support  projectile  design  (Magnus  effect)  and  test  range 
safety  analysis  of  jettisoned  components  (tumbling  motion). 

Tumbling  Motion  -  The  simplest  techniques  found  were  empirically  derived 
equations  and  charts  to  compute  the  average  tumbling  drag  coefficient  of  plates, 
cubes,  spheres  and  cylinders.  This  technique  was  assembled  by  MDAC-HB  in 
Reference  32  and  is  the  recommended  approach.  Other  data  available  is  summarized 
in  Reference  252. 

Magnus  Effect  -  The  effect  of  spinning  motion  has  been  the  subject  of  com¬ 
prehensive  investigations  (References  175,  253-255).  A  computer  program  (SPINNER, 
Reference  208)  is  available  to  numerically  predict  the  Magnus  effect.  A  number 
of  theoretical  methods  include  an  angle  of  attack  dependency.  It  is  recommended 
that  results  from  the  theoretical  and  semi -empirical  techniques  be  compared  with 
the  SPINNER  code  and  test  results  for  accuracy.  A  method  selection  can  be  made 
once  this  quantitative  analysis  is  conducted.  No  specific  recommendation  is  made 
at  this  time. 
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SECTION  11 


CONCLUSIONS  AND  DEVELOPMENT  RECOMMENDATIONS 
11.1  Feasibility  and  Development 

Based  upon  established  realistic  requirements  the  development  of  Missile 
Datcom  is  feasible.  Methodology  is  available  to  perform  preliminary  or  con¬ 
ceptual  missile  aerodynamic  design.  There  is  no  comprehensive  collection  of 
missile  design  methods  available.  The  assembly  and  publication  of  such  methods 
as  a  handbook  is  seen  as  a  major  requirement  to  fill  this  important  void. 

Existing  computer  programs  are  available  but  they  are  often  limited  in 
applicability,  poorly  programmed,  too  complex,  or  are  substantially  undocumented. 
To  overcome  these  problems  Missile  Datcom  should  be  developed  as  follows: 

a)  A  handbook  and  its  companion  computer  program  should  be  developed  in 
parallel.  The  handbook  should  include:  (1)  A  brief,  but  thorough 
discussion  of  the  physical  phenomena  being  modeled,  including  a 
summary  of  the  techniques  available  and  the  reasons  for  choosing 
the  recommended  technique.  (2)  A  description  of  the  selected 
technique  including  all  equations,  tables  and  charts  necessary  for 
its  use.  This  description  is  expected  to  include  the  geometric  and 
flight  condition  range  of  applicability,  data/theory  comparisons, 
and  specific  recommendations  for  method  interpolation  and  extra¬ 
polation.  (3)  A  bibliography  to  allow  more  sophisticated  analysis 
if  required. 

b)  The  computer  program,  developed  concurrently  with  the  handbook, 
should  reflect  as  a  minimum  the  same  capability  as  the  handbook. 

Where  methods  have  been  selected  for  the  program,  that  are  more 
sophisticated,  or  iterative,  a  description  of  the  differences  between 
the  handbook  and  programmed  techniques  should  be  documented.  The 
program  should  have  the  following  characteristics:  (1)  It  should 

be  written  to  conform  with  the  current  American  National  Standards 
Institute  (ANSI)  standard  FORTRAN  computer  code  for  easy  use  on 
many  computer  systems,  (2)  The  code  should  be  structured  using  the 
concepts  of  structured  programming  as  applied  to  the  FORTRAN 
language.  This  will  insure  code  readability.  (3)  The  program 
should  be  well  documented  internally,  with  liberal  but  effective 
use  of  comment  cards.  This  will  aid  in  code  development  and  increase 
its  utility.  (4)  The  program  should  be  modular  with  the  capability 
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to  interchange  subroutines.  Each  subroutine  should  perform  one 
specific  task,  such  as  computing  panel  lift  curve  slope  at  subsonic 
speeds.  A  group  of  routines  which  perform  related  computations,  such 
as  wing  lift  at  angle  of  attack  and  Mach  number,  should  be  assembled 
within  the  same  program  segment.  This  concept  will  allow  easy  sub¬ 
stitution  of  alternate  methodology  or  experimental  results. 

(5)  The  program  should  be  structured  to  minimize  the  number  of 
user  inputs  required.  The  input  scheme  should  be  easy  to  use. 

(6)  The  program  should  be  structured  to  minimize  computer  execution 
time  and  core  requirements,  thus  inexpensive  to  operate.  Setting 
the  maximum  computer  core  at  100,000  octal  words,  as  is  presently 
done  for  Digital  Datcom,  is  a  realistic  goal.  (7)  Finally,  a 
user's  manual  should  be  developed  which  defines  the  inputs  required 
and  the  outputs  available.  It  should  reference  the  handbook  where 
possible.  All  extrapolations  of  the  methods  should  be  clearly 
identified. 

11.2  Recommended  Tasks 

Those  areas  which  require  further  method  development  are  summarized  in 
Table  18.  The  following  tasks  in  decending  order  of  priority  should  be 
performed  to  eliminate  major  deficiences  in  missile  aerodynamic  prediction  at 
the  conceptual  and  preliminary  design  levels. 

A)  Development  of  analytic  methods  for  estimating  the  stability  and 
control  effects  of  airbreathing  inlets.  Present  and  future  missile 
designs  emphasize  the  need  for  longer  cruise  range.  This  requires 
the  use  of  airbreathing  propulsion  systems.  Section  7  details 

the  present  methodology  status. 

B)  Development  of  analytical  methods  to  determine  mutual  interference 
effects  at  all  Mach  numbers  for  the  following  conditions  or  configura¬ 
tions  : 

(1)  angle  of  attack  and  panel  deflection  angles  recommended  in 
Table  2  where  slender-body  or  conical  flow  theory  is  no 
longer  valid.  The  available  theoretical  and  empirical 
methods  are  discussed  in  Section  5. 

(2)  swept  trai ling-edge  straight-tapered  panels;  method  substantia¬ 
tion  required. 

(3)  swept  forward  panels,  such  as  for  oblique  wings  employed  on 
cruise  missile  class  vehicles;  no  methods  available. 
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(4)  non-straight- tape red  panels;  no  methods  available 

(5)  fin  panel  mutual  interference  for  other  than  planar  or  cruciform 
fin  arrangements;  no  methods  available. 

C)  Development  of  analytical  techniques  to  determine  fin  panel  normal 
force  and  center  of  pressure  at  high  angles  of  attack,  at  all  Mach 
numbers.  No  theoretical  methods  are  available.  Empirical  results 
are  limited. 

D)  Development  of  easily  applied  techniques  to  determine  the  stability 
and  control  characteristics  of  arbitrary  shaped  geometries  at  all 
speeds.  Present  methods  are  too  expensive  or  too  difficult  to  use. 

E)  Development  of  analytical  techniques  to  determine  the  wave/pressure 
drag  of  nose  shapes  other  than  axisymmetric  tangent  ogives  or  cones. 
Data  availability  is  limited  at  transonic  Mach  numbers. 

F)  Development  of  analytical  techniques  to  determine  the  effect  of 
angle  of  attack,  panel  deflection  and  jet  exhaust  on  configuration 
dynamic  stability.  Very  few  methods  are  available. 

The  listed  deficiencies  cause  significant  shortcomings  in  missile  aero¬ 
dynamic  prediction  capability.  It  is  important  that  these  deficiencies  be 
corrected  as  quickly  as  possible.  Such  an  effort  requires  the  leadership 
and  support  of  the  United  States  Air  Force. 
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TABLE  18  GAPS  IN  METHODOLOGY 


GO 
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NOMENCLATURE 


SYMBOL 
LOWER  CASE 

a 

b 

c 

d 

e 

g 

h 

i 

P 

q 

r 

s 

t 

u 

v 

w 

X 

y 

z 


DEFINITION 

Body  radius 
Fin  span 
Fin  Chord 

Body  diameter;  also  reference  length 
Oswald's  wing  efficiency  factor 
Gravitation  constant 
Altitude 

Interference  factor;  also  fin  incidence  angle 
Planform  shape  parameter;  also  roll  rate,  pressure 
Pitch  rate;  also  dynamic  pressure 
Yaw  rate;  also  body  radius 
Fin  span 

Airfoil  thickness 
Axial  component  of  velocity 
Lateral  component  of  velocity 
Vertical  component  of  velocity 
Longitudinal  distance 
Lateral  distance 
Vertical  distance 


CAPITALS 

B 

D 

E 

I 

K 

M 

P 

Re 

S 

V 


Mach  similarity  parameter 

Maximum  body  diameter;  also  reference  length 

Elliptic  integral  of  the  second  kind 

Vortex  interference  factor 

Constant  of  proportionality 

Mach  number 

Pressure 

Reynolds  Number 

Reference  Area 

Velocity 
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SYMBOL 

COEFFICIENTS 


DEFINITION 


GREEK 

a 


a 1 

B 

6 

Y 

X 

0 

<f> 


Axial  force  coefficient;  body  axis 
Drag  coefficient;  wind  axis 
Local  skin  friction  coefficient 
Incompressible  skin  friction  coefficient 
Total  (average)  skin  friction  coefficient 
Hinge  moment  Coefficient 
Lift  coefficient;  wind  axis 
Section  lift  coefficient;  wind  axis 
Pitching  moment  coefficient 
Normal  force  coefficient;  body  axis 
Yawing  moment  coefficient;  wind  axis 
Leading-edge  suction  coefficient 
Leading-edge  thrust  coefficient 
Side  force  coefficient;  wind  axis 


Angle  of  attack 

Zero  lift  angle  of  attack 

Angle  of  attack  limit  for  linear  lift 

Total  angle  of  attack 

Sideslip  angle 

Deflection  angle,  or  surface  slope 

Ratio  of  specific  heats  for  a  gas;  1.4  for  air 

Panel  taper  ratio;  tip  chord/root  chord 

Panel  twist  angle 

Roll  angle;  or  bank  angle 
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SYMBOL 

DEFINITION 

GENERAL  SUBSCRIPTS 

A.C. 

Aerodynamic  Center 

B 

Body 

(B) 

In  presence  of  body 

C.P. 

Center  of  pressure 

T 

Tail 

(T) 

In  presence  of  tai  1 

Theo 

Theoreti cal 

W 

Wing 

(W) 

In  presence  of  wing 

WB 

Wing- body  combination 

CO 

At  free-stream  conditions 

NON-DIMENSIONAL 

FACTORS 

All  forces  and  moments  are  non-dimensional i zed  by  the  free-stream 
dynamic  pressure  and  the  reference  area.  The  reference  area  is  the  maximum 
cross-sectional  area  of  the  body.  In  addition,  the  moments  are  non- 
dimensionalized  by  the  body  maximum  diameter.  The  reference  area  and 
length  is  the  same  for  both  longitudinal  and  lateral-directional  aerodynamic 
coefficients. 

The  dynamic  derivatives  are  non-dimensionalized  by  the  free-stream 
dynamic  pressure,  the  reference  area  and  the  reference  length.  In  addition, 
the  rate  derivatives  are  non-dimensionalized  using  the  free-stream  velocity 
and  the  reference  length.  For  example,  the  dynamic  derivative  Cmq  is 
non-dimensionalized  using 
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